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Preface 
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of  this  Investigation.  His  wisdom  In  the  realm  of  scientific  inquiry 
will  not  be  forgotten.  I  am  grateful  to  Dr.  James  E,  Hitchcock  and 
Captain  Wesley  R.  Crow  for  the  expertise  that  they  generously  gave  In 
their  respective  fields  of  heat  transfer  and  metallurgy,  .-iso,  I  wish 
to  express  my  appreciation  for  the  assistance  provided  by  Majo  John  L. 
Kurzenberger,  Major  Keith  Gilbert,  and  Captain  William  T.  Laugh! in  of 
the  sponsoring  Laser  Effects  Branch,  Laser  Division,  Air  Force  Weapons 
Laboratory.  And  finally,  my  family  deserves  a  very  special  thanks  for 
their  patience  end  understanding  afforded  me  during  the  prolonged  work 
on  my  thesis. 
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Abstract 

A  two-dimensional  numerical  model  of  the  surface  recession 
dynamics  and  transient  heating  of  melai  plates  simultaneously  sub¬ 
jected  to  high  intensity  laser  radiation  and  aerodynamic  frictional 
forces  is  developed.  The  model  employs  heat  and  force  balance  calcu¬ 
lations  over  successive  finite  time  increments  on  an  array  of  finite 
elements  constructed  within  the  plate's  cross-sectional  profile  lying 
parallel  to  the  wind  vector  and  bisecting  the  laser  beam  spot.  A  com¬ 
puter  program  of  the  model  in  FORTRAN  Extended  is  provided.  The  model, 
exercised  at  various  mach  numbers,  was  found  to  give  results  consistent 
with  simplified  limiting  cases  of  instantaneous  melt  removal.  Melt- 
through  times  for  0.05,  0.10,  and  0.20  centimeter  thick  sheets  of 
aluminum,  magnesium,  stainless  steel,  and  titanium  absorbing  1  and 
10  kilowatts  cf  power  over  spots  of  1  cm2  and  10  cm2  at  mach  numbers 
ranging  from  0.0  to  4.0  are  presented. 
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A  NUMERICAL  MOOEL  OF  SURFACE  RECESSION  PHENOMENA 
OF  METALS  SUBJECTED  TO  LASER  RADIATION 
IN  AN  AERODYNAMIC  ENVIRONMENT 

I .  Introdv-’  tion 

Light  has  long  been  equated  with  illumination;  but  light  as  a 
means  of  transmitting  man-generated  energy  was  only  science  fiction 
until  the  advent  of  the  laser.  Focused  laser  light  can  now  deposit 
energy  at  intensities  heretofore  never  seriously  considered.  The 
virtual  void  of  experience  has  prompted  a  flourish  of  laser  effect 
studies.  This  thesis  presents  a  study  concerning  the  problem  of 
modeling  laser  effects  on  metals  in  an  airstream. 

Statement  of  the  Problem 

The  objective  of  this  study  is  to  model  the  surface  recession 
and  transient  heating  of  metals  in  aerodynamic  ervironments  under  expo¬ 
sure  to  high  intensity  laser  radiation.  The  proposed  model  shall  be 
programmed  for  computer  solution  and  be  adaptable  to  various  mQtals  and 
alloys  under  varying  initial  and  boundary  conditions. 

Significance  of  the  Problem 

No  explicit  means  is  Known  to  exist  whereby  melt-through  or 
"melt-to"  times  can  be  predicted  for  laser  beams  incident  upon  metals 
over  which  air  is  flowing.  A  model  to  predict  such  times  and  to  gen¬ 
erate  ancillary  data  describing  the  transient  heating  and  surface 
recession  dynamics  will  enable  various  analyses  of  material  response 
to  laser  heating.  Varying  the  initial  parameters  such  as  laser  power, 
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beam  diameter,  exposure  times,  r£tal  thicknesses  .and  properties,  and 
aerodynamic  conditions  (definable  by  mach  number  and  altitude  in  a 
"standard"  atmosphere)  in  successive  simulations  can  provide  information 
necessary  to  an  increased  understanding  of  the  relationships  involved. 
Corresponding  data  through  experiment  is  and  can  be  expected  to  continue 
to  be  minimal-  due  to  nonexistence  of  lasing  devices  In  various  power/ 
spot  size  regimes,  difficulty  in  accurately  and  fully  instrumenting 
actual  experiments,  and  excessive  cost  and  time  required  to  conduct 
sufficient  tests  over  the  range  of  parameters  of  possible  interest. 
Development  of  a  model  Is  thus  considered  to  be  an  appropriate  sup¬ 
plement  to  laser  effects  tests  that  are  possible  and  a  necessary 
substitute  for  tests  which  are  not  yet  feasible. 

Definition  of  Terms 

The  symbols  used  in  this  text,  the  corresponding  mnemonic  names 
used  in  the  resultant  computer-programmed  model,  and  their  definition? 
are  listed  on  p?ge  viii.  Laser  radiation  as  treated  herein  is  electro¬ 
magnetic  radiation  within  the  infrared- visible-ultraviolet  spectrum. 
"High  intensity"  refers  to  a  laser  energy  flux  sufficient  to  create  sig¬ 
nificant  heating  effects,  namely  melting,  when  incident  upon  a  metal 
under  consideration;  a  precise  numerical  definition  is  not  intended. 
Other  terms  will  be  defined  as  they  occur. 

The  Physical  Phenomena  of  Laser 
Interaction  with  Metals' 

Given  a  laser  beam  directed  toward  a  metal,  only  a  fraction  of 
the  photons  constituting  the  initial  beam  energy  emitted  from  the  lasing 
device  will  be  absorbed  by  tiie  metal.  First,  the  intervening  medium 
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inappreciably  attentuate  the  transmitted  energy.  The  degree  of  such 
attentuatiOR  is  dependent  upon  the  following:  the  density  and  composi¬ 
tion  of  the  medium  (for  example,  atmospheric  humidity  can  be  a  strongly 
influential  factor);  the  wavelength  of  the  photon;  the  “f)ux  intensity 
within  the  beam;  and  the  distance  traversed  by  the  beam.  Secondly,  a 
portion  of  the  incident  photons  are  reflected  oy  the  metal.  The  por¬ 
tion  reflected  is  a  function  of  several  factors;  the  photon  wavelength; 
the  incidence  angle;  the  metallic  composition;  the  surface  characteristics 
including  degree  of  roughness,  foreign  coatings  (ranging  for  example 
from  a  simple  inadvertent  fingerprint  to  a  purposeful  protective  paint), 
surface  scale  from  oxidizing  reactions  involving  the  metal;  and  surface 
temperature  which  in  turn  governs  the  physical  state  (solid  or  liquid) 
of  the  metal's  face. 

The  phenomenon  by  which  the  absorbed  photons  impart  heat  to  the 
metal  is  primarily  by  the  successive  mechanisms  of  internal  photo¬ 
electric  effect  and  electron  de-excitation  to  either  optic  or  acoustic 
phonons  (Ref  1:1-2).  The  process  occurs  essentially  instantaneously  a1: 
the  immediate  surface  of  the  metal.  The  accuracy  of  this  approximation 
can  be  evaluated  by  considering  a  "worst  case"-~one  in  which  penetration 
by  the  photon  should  be  maximal --where  highly  energetic  photons  impinge 
upon  a  comparatively  light  metal.  For  the  case  of  magnesium  subjected 

O 

to  far  ultraviolet  photons  of  12.4  angstrom  (A)  wavelength  X,  the 
photons  have  an  energy  E  expressed  by 


where  h  is  Planck's  constant  ( h  -  4.136><10"15  evsec),  and  o  is  the 
speed  of  light  (e  =  3.00x10®  m/sec).  In  this  example  E  equals  1  Kev 
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(one  thousand  electron  volts).  From  Ref  2,  the  total  mass  attenuation 
coefficient  y/p  for  1  Kev  photons  is  993.3  cm2/gm  with  the  photo¬ 
electric  effect  accounting  for  all  but  0.21%  of  the  value  of  the  mass 
attentuation  coefficient  (Ref  2:49).  The  half- thickness  fi/2,  which 
is  the  penetration  depth  by  which  one-half  of  all  the  nonreflected 
photons  are  absorbed  is  (Ref  3:92) 


whe  is  the  linear  attenuation  coefficient  found  by  multiplying  the 
mass  ac  nation  coefficient  previously  noted  by  the  density  p  of  the 
metal  (p  for  magnesium  is  1.740  gm/cm3).  Thus,  for  magnesium  ti/2 
equals  0.0004  cm,  a  very  small  distance  relative  to  even  thin  sheets  of 
metals  used  in  practical  applications. 

Electrons  excited  by  the  photoelectric  absorption  mechanism  in 
turn  de-excite  by  transferring  their  excess  energy  to  phonons,  or  vibra¬ 
tional  energy  within  the  crystalline  lattice  of  the  solid,  and  thereby 
complete  the  transformation  of  light  to  heat.  The  total  time  of  photon 
to  phonon  conversion  is  the  time  to  absorption  (after  passing  the  face 
of  the  metal)  plus  the  electron  de-excitation  time.  The  order  of  mag¬ 
nitude  of  the  time  to  absorption  may  be  calculated  by  dividing  the 
half- thickness  t\/z  by  the  speed  of  light;  and,  in  the  "worst  case" 
considered  above,  the  time  to  absorption  is  about  1 . 33x10“ 1 4  sec.  The 
"worst  case"  for  de-excitation  times,  that  is  where  the  relaxation  takes 
the  longest,  occurs  where  the  incident  photons  have  much  less  energy 
than  the  1  Kev  used  to  calculate  the  half-thickness  ti/Z  above  (Ref  1:21) 
It  is  also  hypothesized  that  transition  metals  have  longer  relaxation 
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times  than  do  simple  metals.  Typical  de-excitation  times  have  been 
calculated  to  be  about  10'7  to  10“6  sec  for  transition  metals  such  as 
stainless  steel  and  titanium  compared  to  IQ”11  to  10~10  sec  for  simple 
metals  such  as  aluminum  and  magnesium  (Ref  l:v). 

The  combined  times  to  absorption  and  to  de-excitation  can  be 

approximated,  due  to  its  dominance,  solely  by  the  time  to  de^excite; 

that  is,  the  "worst  case"  total  time  for  conversion  of  light  to  heat  is 
~6 

about  10  sec.’  Such  times  are  large  when  compared  to  the  duration  of 
a  typical  single  IQ"8  sec  pulse  from  a  (^-switched  laser  (Ref  4:8),  and 
the  assumption  of  instantaneous  light  to  heat  conversion  in  this  case 
is  inappropriate;  but  for  continuous  wave  (cw)  lasers,  for  pulsed  lasers 
where  the  individual  pulse  is  much  longer  than  10"6  sec,  and  possibly  for 
repetitively  pulsed  lasers  where  the  pulse  repetition  rate  is  very  high, 
the  assumption  of  "instantaneous"  conversion  of  light  to  heat  is  approxi¬ 
mate  for  most  analytic  purposes. 

Subsequent  to  the  optical -to- thermal  energy  transformation,  heat 
is  transferred  from  the  laser  beam-metal  interface,  or  "spot,"  by  four 
mechanisms,  three  of  which  are:  conduction  within  the  metal,  radiation 
from  the  surface,  and  convection  to  (or  from)  the  adjoining  air.  If 
the  absorbed  energy  flux  is  sufficiently  intense,  a  finite  element 
within  the  laser  spot  receiving  the  greatest  flux  will  after  some  initial 
interval  of  time  reach  the  melting  temperature;  and  after  accumu  ating 
an  additional  amount  of  energy  (the  latent  heat  of  fusion)  the  element 
transforms  to  the  molten  state.  As  a  liquid,  the  metal  is  subject  to 
mass  flow  due  to  the  surface  shear  forces  created  by  aerodynamic  skin 
friction;  this  is  the  fourth  mechanism  by  which  heat  is  transferred. 
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The  rate  of  the  liquid  transfer  (and  consequential  heat  transfer}  is  a 
function  of  the  magnitude  of  the  drag  forces,  the  viscosity  of  the 
liquid  metal,  and  the  depth  of  the  molten  layer.  Transfer  of  the  melt 
along  the  surface  of  the  metal  to  a  region  of  lower  flux  absorption  an.: 
correspondingly  cooler  surface  temperatures  may  result  in  either 
freezing  of  the  liquid  metal  or  melting  of  the  underlying  solid  layer. 
As  time  progresses,  heat  will  continue  to  be  absorbed  and  conducted 
such  that  interior  temperatures  will  continue  to  rise  and  additional 
incremental  layers  will  continue  to  melt  thereby  causing  the  solid- 
liquid  interface  to  recede  until  melt-through  occurs. 

Additional  phenomena  which  can  further  complicate  tha  total 
modeling  problem  include  vaporization  of  the  'liquid  metal  and  oxidation 
at  the  metal's  surface.  Vaporization  can  occur  when  incident  flux 
exceeds  the  combined  modes  of  heat  dissipation.  If  vaporization  does 
occur,  any  vapor  within  the  laser  beam  may  strongly  attenuate  the 
oncoming  laser  energy  flux  (Ref  1:2).  Oxidizing  reactions  liberate 
heat  of  combustion  at  rates  dependent  upon  the  temperature-time  history 
of  the  metal  and  the  ambient  air  and  wind  conditions.  The  rates  may 
be  insignificant  except  for  certain  metals  such  as  magnesium  and  tita¬ 
nium,  where  at  elevated  temperatures  in  a  high  velocity  airstream 
catastrophic  oxidation  may  occur.  The  conditions  under  which  this 
ignition  begins  are  not  well-defined  (Ref  5:1),  but  when  ignited  the 
liberated  heat  of  combustion  combined  with  the  laser  flux  will  accel¬ 
erate  the  melting  process  to  one  of  burning.  The  rate  of  the  process 
is  highly  dependent  upon  the  oxygen  available  and  the  temperature  of 
the  combustive  metal;  thus  the  combustion  should  propagate  faster  to 
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the  leeward  but  may  not  be  self-sustaining  beyond  the  area  influenced 
by  the  laser  radiation. 

Scobe  of  Study 

This  study  includes  the  two-dimensional,  numerical  modeling  of  , 
the  surface  erosion  dynamics  and  the  transient  heating  within  metal 
plates  given  the  mach  number,  the  altitude,  the  distance  from  the  plate's 
leading  edge,  the  energy  absorption  rate  (flux)  distribution,  and  the 
physical  and  thermodynamic  properties  of  the  metal.  Attenuation  of  the 
laser  energy  from  the  lasing  device  to  th;e  plate  and  reflection  of  the 
light  at  the  surface  of  the  metal  are  not  treated  in  this  work.  Trans¬ 
port  of  molten  metal  d >e  to  gravity  and  energy  release  from  any  surface 
oxidation  or  combustion  were  not  studied  sufficiently  to  permit  inclusion 
of  their  respective  effects  in  the  resultant  model.  Only  the  shear 
stresses  from  aerodynamic  skin  friction  are  incorporated  within  the 
model;  pressure  drag  forces  on  the  molten  metal  are  neglected.  The 
model  assumes  an  axial!"  symmetric  flux  distribution  within  a  norma 'ly 
incident  laser  beam.  Heat  conduction  in  the  axial  ?,id  windward/leeward 
directions  is  included;  heat  conduction  perpendicular  to  the  plane 
formed  by  the  wind  vector  and  the  normal  to  the  plate  surface  is 
disregarded. 

Four  metals  (aluminum,  magnesium,  stainless  steel,  end  titanium) 
were  specifically  modeled,  including  variation  of  their  physical  and 
thermodynamic  properties  as  a  function  of  temperature.  Other  metals 
can  be  accommodated  upon  appropriately  incorporating  their  respective 
properties  within  the  model.  Although  it  was  originally  intended  to 
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model  continuous  wave  (cw)  laser  radiation,  the  resultant  model  appears 
adaptable  to  certain  pulsed  laser  radiation  where  the  duration  and 
repetitive  rate  of  the  pulses  are  consonant  with  the  approximation  of 
instantaneous  photon  to  phonon  conversion  as  discussed  earlier  in  this 
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Complicating  Factors 

The  problem  of  modeling  the  heat  transfer  within  two-phase 
metal  has  defieJ  a  general  solution  (Ref  6:477).  The  primary  obstacle 
has  been  the  complexity  of  the  mathematics  involved  (Ref  7:951).  Com¬ 
plicating  factors  include  the  following:  icnree-dimensional ,  unsteady 
temperature  distributions  within  the  solid  and  liquid  phases;  temperature- 
dependent  thermodynamic  properties;  diverse  initial  and  boundary  condi¬ 
tions  fauna  in  practical  applications;  and  the  presence  of  two  moving 
physical  boundaries. 


II.  Analysis  of  the  Problem 

Although  lasers  have  provided  a  novel  mode  of  transmitting  inten¬ 
sive  energies  by  means  of  photons,  analysis  of  the  surface  recession 
and  transient  heating  of  metals  coincidentally  subjected  to  such  laser 
energy  and  aerodynamic  friction  reduces  to  using  fundamentals  of  heat 
transfer  and  mass  flow— established  fields  of  science.  Justification  for 
such  analytic  methods  lies  in  the  nature  of  light  absorption  by  metals, 
as  described  in  Chapter  I,  whereby  the  absorbed  light  is  approximately 
instantaneously  converted  to  thermal  heat  at  the  immediate  surface  of 
the  metal.  Such  thermal  heat  is  subject  to  all  of  the  conventional 
equations  employed  in  the  field  of  heat  transfer.  Thus  the  total  phe¬ 
nomenon  being  modeled  conveniently  divides  Itself  into  (1)  a  problem  of 
heat  transfer  within  a  static  two-  (or,  including  the  vapor  phase,  a 
three-)  phase  metal,  and  (2)  a  problem  of  mass  flow  within  the  liquid 
metal  resulting  from  shear  forces  generated  by  aerodynamic  friction 
over  the  plate's  surface. 
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The  second  problem,  that  of  determining  the  molten  mass  flow 
rate  and  skin  frictional  heating  due  to  aerodynamic  flow  over  the  metal 
plate  relates. to  boundary  layer  studies  that  have  been  extensively  pur¬ 
sued  since  the  requirement  arose  for  determining  heat  and-mass  transfer 
effects  on  reentry  vehicles.  While  even  though  most  existent  solutions 
to  these  problems  are  approximate  ones,  it  does  appear  that  adequate 
investigative  study  has  been  accomplished  for  now  and  that  the  relations 
obtained  have  been  verified  sufficiently  to  give  reasonable  tools  for 
the  purpose  of  this  study  (Refs  8  and  9).  Factors  that  in  general  can 
complicate  the  calculations  of  drag  forces  and  convective  heat  transfer 
due  to  aerodynamic  friction  are  as  follows:  arbitrary  shapes  of  the 
objects  to  be  modeled,  imprecise  ability  to  predict  the  transitional 
Reynolds  number  between  laminar  and  turbulent  flow,  and  locally  varying 
temperatures  along  the  line  of  flow. 

Previous  Studies  and  Current  Knowledge 

A  study  on  polar  ice  thickness  including  a  mathematical  treatise 
on  the  associated  heat  conduction  with  freezing  was  published  by  Josef 
Stefan  in  1891:  and,  although  it  was  Franz  Neumann  who  first  found  the 
particular  solution  for  freezing  of  a  semi-infinite  slab  with  a  constant 
boundary  temperature,  the  general  problem  of  conduction  with  freezing  or 
melting  has  become  known  as  "Stefan's  Problem."  Muehlbauer  and 
Sunderland  (Ref  7),  the  source  of  the  above  historical  note,  have  pub¬ 
lished  a  comprehensive  literature  review  of  works  covering  "Stefan's 
Problem"  to  1965.  Reference  10  refers  to  a  survey  of  later  works. 

While  Stefan's  Problem  is  relevant  to  numerous  contemporary  pro¬ 
cesses,  perhaps  the  most  studied  of  these  have  been  metal  casting, 
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welding,  and  reentry  vehicle  ablation.  Despite  the  general  similarity 
of  all  these  problems  to  that  of  'laser  heating,  subtle  differences  do 
exist;  and  no  single  solution  to  any  of  these  problems  was  found  to  be 
readily  adaptable  to  modeling  melt-through  by  laser  heating. 

Several  mathematical  approaches  have  been  employed  in  seeking 
solutions  to  Stefan's  Problem.  They  include  the  variational,  heat- 
balance  and  contour  integral,  analogue,  and  numerical  techniques 
(Ref  7:951-6).  Of  these,  numerical  analysis  using  the  finite  differ¬ 
ence  and  finite  element  approaches  was  selected  for  the  purposes  of 
this  study— primarily  on  the  basis  of  the  availability  of  a  large, 
high-speed,  digital  compute/  (such  a  computer  is  nearly  mandatory  for 
the  large  array  of  elements  necessary  for  any  satisfactory  degree  of 
accuracy).  This  approach  in  part  paralleled  the  work  of  Torvik  (Ref 
10)  who  has  recently  modeled  laser  heating  exclusive  of  aerodynamic 
erfects. 

Method  of  Attacking  the  Problem 

The  final  model  of  this  study,  as  found  in  Appendix  A,  evolved 
from  a  rudimentary,  one-dimensional,  constant  flux,  constant  thermo¬ 
dynamic/physical  properties,  instantaneous  melt  removal,  finite 
difference  model  by  progressively  removing  the  enumerated  qualifi¬ 
cations.  The  next  two  chapters  review  the  logic  employed  in  construct¬ 
ing  the  final  mooel.  Chapter  III  deals  with  heat  transfer  by  absorption, 
radHtion,  conduction,  and  vaporization,  and  Chapter  IV  with  heat 
transfer  by  convection  and  melt  flow. 
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III.  Modeling  the  Heat  -Transfer 

Assuming  a  normally  incident  laser  beam  in  still  air  with  axially 
symmetric  energy  deposition  *  no  heat  transfer  occurs  circumferentially 
about  the  beam  axis..  An  array  consisting  of  concentric  annuli,  as 
shown  in  Fig.  1,  can  be  used  to  model  the  resulting  transient  heating 
and  progressive  melting  phenomena  two-dimensionally  (radially  and 
axially)  in  a  cylindrical  coordinate  system  with  no  loss  in  accuracy 
occurring  by  ignoring  the  circumferential  direction  ("'.ef  10:6). 

The  presence  of  wind-driven  melt  flow,  however,  contravenes  the 
two-dimensionality  of  the  total  surface  recession  problem  since  the 
axial  symmetry  of  temperatures  is  lost  as  soon  as  melting  commences. 
Nevertheless,  thermal  and  topographical  symmetry  should  theoretically 
be  maintained  about  the  spot-bisecting,  wind-parallel  line.  Furthermore, 
deepest  penetration  should  occur  within  the  plane  directly  beneath 
this  same  line.  A  two-dimensional  array  to  become  the  basis  of  the 
finite  element,  numerical  approach  to  modeling  the  total  phenomenon 
was  therefore  constructed,  as  shown  in  Fig.  1,  from  the  volume  lying 
directly  beneath  an  infinitesimal  strip  lying  along  the  wind  line  and 
intersecting  the  beam  axis.  The  origin  of  the  array  is  taken  to  be  the 
element  furthest  downwind  at  the  back  surface  of  the  metal  plate.  The 
first  subscript  of  an  ?-,jth  element  refers  to  the  column  and  the  second 
to  the  row  in  the  array  where  the  array  is  divided  into  m  columns  and 
n  tiers.  With  an  implement  for  the  numerical,  finite  element  analysis 
thus  established,  distribution  of  the  flux  absorption  is  next  modeled. 
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Fig.  1.  Top  and  Cross-Sectional  Views  of  a  Metal  Plate  Impinged  by  a 
Gaussian  Laser  Beam  of  Radius  2a,  The  radius  thus  defined  is 
the  distance  at  which  the  local  flux  density  is  0.135  of  the 
intensity  at  the  beam  axis  as  is  illustrated  by  flux  distri¬ 
bution  along  the  cellularly  arrayed  cross  section  of  plate 
underlying  the  beam-bisecting,  wind-parallel  strip. 
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Modeling  the  Flux  Absorption  Distribution 

The  Incident  laser  energy  was  assumed  -to  be  absorbed  as  a 
Gaussian  distribution;  that  is,  the  intensity  at  any  point  on  the  sur¬ 
face  of  the  metal  is  - 


?£&*)  =  <7q*‘ 


Cl/2) (r2/o2) 


(3) 


where  is  the  absorbed  flux  intensity  at  the  beam  axis  which,  in  turn, 
is 


<?o  = 


2Q 

7i(2a)2 


(4) 


Q  is  *he  total  absorbed  flux  from  the  lasing  device,  and  2a  is  defined 

as  the  radius  of  the  laser  beam  spot.  The  spot  radius  defined  in  this 

manner  is  the  distance  from  the  beam  axis  at  which  the  locally  absorbed 
-2 

flux  intensity  is  e  ,  or  13.5%,  of  its  value  at  the  axis.  The  spot  so 
defined  absorbs  approximately  86.5%  of  the  total  absorbed  incident 
laser  energy. 

The  above  representation  is  reputed  to  approximate  the  spatial 
distribution  for  the  flux  of  cw  gas  lasers  in  general  (Ref  4:13).  Ex¬ 
perience  has  shown  that  in  actuality  hot  spots  can  occur  within  such 
laser  beams;  but,  because  of  the  unpredictability  of  the  intensity, 
size,  and  location  within  the  beam  of  such  hot  spots,  the  Gaussian  dis¬ 
tribution  was  deemed  the  most  appropriate.  If  more  precise  information 
is  available  on  the  spatial  flux  distribution  from  a  particular  lasing 
device,  the  affected  subroutine  (FLUXX)  within  the  computer-programmed 
model  (ree  Appendix  A)  may  be  modified  accordingly. 
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Modeling  the  Internal  Heat  Transfer 

The  transient  heat  conduction  within  the  metal,  which  transfers 
energy  both  axially "and  radially  from  the  area  of  deposition  (the  laser 
spot),  is  mathematically  described  in  Cartesian  coordinates  for  a 
volume  element  by  the  following  parabolic  partial  differential  energy 
balance  equation  (Ref  11:31): 


3  T 

<2  P  - 

P  H 


(5) 


where 

T  =  temperature  ("Kelvin), 

^  =  specific  heat  (cal/gm/°K), 
p  =  density  (gm/cm3), 
k  =  conductivity  (ca1/sec/cm/°K) , 
t  =  time  (sec), 

a  =  rate  of  energy  gain/loss  due  to  laser  flux  absorption, 
surface  radiation,  or  convection  ( cal /sec/cm3 ) ,  and 

x,y,z  -  orthogonal  directions. 

The  restriction  in  dimensionality  to  the  analysis,  as  imposed 
by  the  two-dimensional  finite  element  array,  in  essence  ignores  the 
heat  conduction  in  the  z  direction  when  x  and  y  are  taken  as  the  wind¬ 
ward  and  axial  directions  respectively  (or,  equivalently,  the  restric¬ 
tion  assumes  an  invariant  flux  q  in  the  z  direction).  Some  heat  is 
indeed  transferred  in  this  direction  since,  even  though  the  crosswind 
temperature  gradient  at  the  plane  of  symmetry  is  theoretically  zero, 
its  derivative  is  not.  The  consequences  of  ignoring  conduction  in  the 
z  direction  are  slightly  faster  melt-through  times  and  somewhat  elevated 
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temperatures a  particularly  close-in  to  the  beam  axis  where  32T/3z2  is 
larger  in  magnitude  than  at  outlying  radii.  The  two-dimensional  model¬ 
ing  thus  becomes  more  nearly  exact  with  increasing  distance  from  ;the 
beam  axis  and  with  increasing  spot  size.  The  approximation  is,  however, 
mere  precise  than  a  strictly  one-dimensional  model  which  in  turn,  as  is 
indicated  in  Ref  12,  gives  close  approximations  to  three-dimensional 
conduction  for  large  total  absorbed  powers  in  thin  metal  sheets.  From 
these  considerations,  a  two-dimensional  finite  element  array,  as  con¬ 
structed,  should  not  in  itself  severely  restrict  the  resultant  model's 
applicability. 

Using  the  numerical  technique  of  finite  elements,  F.q  (5)  can  be 
reduced  to  a  more  tractable,  algebraic  equation  for  an  i.jth  element  as 
follows: 

AT.  .o  p  Ax  Ay  As  =  A t*HNET .  .  (6) 

IjJ  p  t,J 

where  AT.  .  is  the  temperature  change  of  the  i,j th  element  over  the 
time  period  At;  Ax  and  A u  are  the  fixed  dimensions  as  chosen  for  the 
finite  elements  (A z  is  taken  to  be  a  unit  distance  in  this  and  subse¬ 
quent  equations  so  as  to  maintain  dimensional  integrity);  and  hhet.  . 

I'}  J 

is  the  net  heat  transfer  rate  into  the  i.jth  element. 

This  numerical  equivalent  to  the  parabolic  partial  differential 
transient  heat  conduction  equation  (Eq  (5))  requires  an  iterative  solu¬ 
tion  over  finite  time  differences  At.  For  stability  and  convergence  of 
the  numerical  solution  it  is  necessary  that  the  forward  time  increment 
At  in  each  iteration  be  limited  by 


/ _ 

■■  (3 /Ax2) 


1 


+ 


(1/Ay2) 


At  < 


(7) 
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lie  net:  heat  transfer  rate  esst^  .-  is  the  summation  of  the  heat 
flow  rates  through  the  top,  bpttom, right,  and  left  of  the  element  as 
illustrated  for  the  middle  column  of  a  3x3  array  in  Fig.  2.  The  heat 
flows  through  the  right  face  HXFRE  and  through  the  top  eXFRT  for  an  i,jth 
embedded  element  are 


radiation 

convection 


1 


HXl’RB  -  0 


Heat  Transfer  Rates  From  the  Top,  Bottom,  Left,  and  Right 
Constitute  the  Net  Heat  Transfer  Rate  for  Each  Finite 
Element  in  the  Analytic  Cellular  Array 
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The  corresponding  equations  for  the  heat  transfer  through  the  Isfl  face 
HXFRL .  .  and  bottom  face  HXFRB,.  •  are  identical  to  Eqs  (8)  and  (9)  with 
the  exception  that  the  plus  signs  in  the  temperature  subscripting  are 
replaced  with  minus  signs.  Adiabatic  conditions,  which  are  assumed  at 
the  back  surface  of  the  metal  .plate,  require  that  all  HXFRB .  .  values 

Vy  X 

equal  zero.  Left  or  right  faces  of  elements  which  are  coincident  with 
the  respective  ends  of  the  analytic  array  or  which  form  a  portion  of 
the  boundary  of  the  melted  hole  or  lip  are  likewise  assumed  to  be 
adiabatic. 

The  conductivity  k  in  Eqs  (8)  and  (9)  is  temperature-dependent; 
but,  since  two  finite  elements  and  therefore  two  temperatures  are 
associated  with  each  of  the  heat  transfer  equations,  calculation  of  an 
average  conductivity  between  the  centers  of  adjacent  elements  is  called 
for.  Such  an  expression,  derivable  from  the  concept  of  thermal  resis¬ 
tance,  is  analogous  to  the  conductance  of  two  series  resistors  in  an 
electrical  circuit  (Ref  11:10),  and  for  two  radially  adjacent  elements 


k 


WK. .  a/k.tlj 


(10) 


1,0 


The  average  conductivity  for  two  axially  adjacent  elements  is  identical 
except  for  appropriate  changes  of  subscripts. 


Modeling  the  Heat  Transfer  at  the  Surface 

The  non-embedded  elements  at  the  front  surface  in  the  finite 
element  array,  in  addition  to  transferring  heat  by  conduction  to  adja¬ 
cent  elements,  absorb  laser  energy,  radiate  heat,  and  transfer  by 
convection  an  additional  amount  of  heat.  The  rates  of  transferring 
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energy  by  these  additional  means  were  lumped  Into  the  term  "q"  In 
Eq  (5) 

Energy  deposition  from  the  lasing  device  was  modeled  for  an 
assumed  Gaussian  distribution  in  an  earlier  subsection.  The  absorption 
rate  for  a  surface  element,  using  Eq  (3),  is 

w  •  hz 

«  — 2—  e-20'/m  te  (11) 

it  (2a)  ’ 

where  r  is  the  distance  from  the  beam  axis  to  the  center  of  the  element 
and  the  other  terms  are  as  previously  defined. 

The  total  hemispherical  radiative  heat  transfer  rate  from  any 
surface  element  is 

q  A.  .  .  =  e oTu  Am  As  (12) 

‘radiation 

where  e  is  the  emissivity  of  the  metal  (nondimensional),  T  is  the  sur- 

“*  8 

face  temperature,  and  a  is  the  Stephan-Boltzmann  constant  (=  5.67*10 
w/m2/0^).  The  emissivity  e  is  a  function  of  the  metallic  and  alloying 
elements,  the  surface  temperature,  and  the  surface  condition.  Emissivity 
values  in  general  were  found  to  be  both  limited  and  inconsistent  {see 
Appendix  C). 

Modeling  the  convective  heat  transfer  rate  a  is  closely 

^convection 

allied  to  determining  surface  shear  stresses  from  aerodynamic  friction. 

The  derivation  of  the  expression  for  q  is  therefore  deferred 

^convection 

(see  Eq  (23))  until  Chapter  IV. 

Combining  these  three  modes  of  heat  transfer,  the  totcl  energy 
transfer  rate  for  Hie  exposed  face  of  a  surface  element  is 
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HXPRT  =  q„  +  <?  ,.  .  .  +  4  :  .  .  .  (13 

^flux  ^radiation  ^convection 

It  should  be  noted  that  a  _  may  be  either  positive  or  negative 

•convection 

depending  upon  whether  aerodynamic  heating  or  cooling  takes  place. 

Modeling  the  Solid-Liquid  Phase  Change 

The  transition  from  the  solid  to  liquid  phase  occurs  differently 
for  alloys  than  for  pure  metals  in  that  while  the  temperature  of  a  pure 
metal  remains  at  a  constant  melting  temperature  throughout  a  steady 
influx  of  the  latent  heat  of  fusion,  the  temperature  of  an  alloy  rises, 
quite  slowly  initially  but  later  at  an  increased  rate,  from  a  solidus 
to  a  liquidus  temperature.  This  phenomenon  was  modeled  by  means  of  a 
quasi -phase,  coined  the  "slush"  phase,  whereby  the  heat  of  fusion  was 
accounted  for  by  incorporating  it  with  the  specific  heat  of  the  metal 
to  give  a  "pseudo-specific  heat"  cpslush*  defined  as 


'Pslush  2 


=  2  (°P@T 


+  aP@T  -e 

solidus  liquidus 


Heat  of  Fusion 


(T  --  T  ) 

liquidus  solidus 


(14) 


To  better  approximate  the  increasing  rate  of  the  temperature  rise,  a  more 
elaborate  formula  could  be  devised;  but,  more  simply,  the  solidus  tem¬ 
perature  may  be  manipulated  upward  so  as  to  in  effect  account  for  the 

entire  temperature  rise  over  an  upper  portion  of  the  T  , 

solidus-iiquidus 

temperature  range. 


Modeling  the  Vaporization 

A  sufficiently  intense  laser  flux  will  vaporize  a  portion  of  the 
subjected  metal  during  the  melt-through  process.  The  lower  limit  of  the 
flux  intensity  necessary  to  cause  vaporization  is  dependent  pon  the  air 


20 


GAW/MC/72-17 


flow  as  well  as  the  thermophvsical  and  thermodynamic  properties  of  the 
metal.  Reduced  air  flow  and  a  poorly  conductive,  highly  vi.cous  metal 
having  a  small  latent  heat  of  vaporization  and  a  low  specific  heat  in 
the  liquid  phase  enhances  vaporization  at  lower  flux  intensities. 

Using  the  postulation  of  instantaneous  vapor  removal,  no  vapor 
phase  need  be  provided  for  in  modeling  the  melt-through  process.  Instead, 
the  mass  and  thermal  energy  removed  from  finite  elements  undergoing 
vaporization  can  be  fully  accounted  for  by  allowing  the  temperatures  of 
such  elements  to  "overshoot"  the  vaporization  temperature.  Each  "super¬ 
heated"  element  will  possess  energy  in  excess  to  the  heat  content  of 
the  same  element  at  the  vaporization  temperature.  The  excess  energy 
divided  by  the  latent  heat  of  vaporization  gives  the  mass  removed  by 
vaporization.  A  subroutine  "VAPORIZ"  (see  Appendix  A),  within  the  com¬ 
puter  programmed  model  determines  if  there  are  any  "superheated"  elements; 
how  much  mass,  or  equivalently  thickness,  is  removed;  appropriately 
reduces  the  thickness  of  the  affected  elements;  and  sets  the  temperature 
of  the  remaining  mass  to  the  vaporization  temperature. 
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IV.  Modeling  the  MeTt  Transfer 

The  slush  phase,  as  modeled  in  the  previous  chapter,  is  considered 
to  consist  of  a  solid  matrix  interspersed  with  liquid  metart.  Thus  the 
slush  phase  is  assumed  to  be  unaffected  by  surface  shear  forces  generated 
by  aerodynamic  friction.  Once  at  the  liquidus  temperature,  however,  the 
molten  metal  is  subject  to  the  influence  of  aerodynamic  drag  forces; 
and  at  this  point  several  postulations  will  facilitate  the  ensuing 
analysis  of  the  melt  flow  process. 

Modeling  the  Melt  Transfer— Assumptions 

First,  for  the  purposes  of  this  study,  the  airflow  is  assumed  to 
be  parallel  to  the  profile  of  the  hole  as  the  hole  deepens  and  enlarges 
with  continued  melting  (see  Fig.  3).  This  approximation  will  overestimate 
drag  forces  on  the  windward  edge  of  the  hole  and  the  leeward  side  of 
the  lip  (see  Fig.  3).  Similarly,  pressure  forces  on  the  lip's  windward 
side  are  neglected.  As  the  diameter-to-depth  ratio  of  the  hole  decreases, 
this  assumption  becomes  less  plausible.  While  the  qualitative  effect 
of  this  approximation  on  the  hole's  cross-sectional  profile  is  easily 
surmised  (i.e.,  the  modeled  hole  and  lip  will  appear  deeper  and  steeper, 
respectively,  and  the  nadir  of  the  modeled  hole  will  lie  slightly  more 
windward  than  should  actually  occur),  the  quantitative  effect  of  the 
parallel  air  flow  approximation  on  the  melt-through  time  is  not  so 
apparent.  The  insulative  effect  of  the  thicker  modeled  lip  should 
counteract  to  some  extent  the  heat  flow  enhanced  by  the  thinner-than- 
actual  liquid  layer  within  the  windward  side  of  the  modeled  hole.  A 
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Fig.  3.  The  Aerodynamic  Boundary  Layer  Creates  Shear  Forces  Over  the 
Plate's  Surface  Causing  the  Molten  Metal  (Cross-Hatched)  Flow 


slightly  faster- than-actual  melt-through  time  is  probably  the  overall 
effect  of  the  assumed  parallel -to-the-surface  airflow. 

Secondly,  it  is  postulated  that  the  molten  metal  adheres  to  the 
solid  (or  slush)  surface  of  the  plate  rather  than  entering  the  airstream 
as  droplets.  This  is  believed  to  be  a  good  approximation  for  flow  from 
the  leading  edge  of  the  hole  to  as  far  as  the  lip.  At  that  point 
pressure  forces  acting  <n  the  leading  edge  of  the  lip,  combined  with  the 
frictional  drag  forces,  conceivably  could  exceed  the  opposing  surface 
tension  and  internal  shear  forces;  if  so,  consequent  "blow-off"  would 
result.  The  effect  of  blow-off  would  be  a  smaller-than-predicted  lip 
size  and  faster,  but  probably  insignificantly  so,  melt-through  times. 

The  argument  for  this  latter  effect  is  that  a  thinner  lip  would  allow 


/ 
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slightly  faster  conductive  heat  transfer  to  the  underlying  and  inwardly 
situated  material. 

A  third  assumption  made  to  facilitate  the  analysis  is  that  of 
instantaneous  melt  flow  acceleration.  Thjat  is,  it  is  assumed  that  the 
flow  of  molten  metal  reaches  steady  state  in  a  time  that  is  small  com¬ 
pared  to  the  total  time  of  the  flow.  While  this  approximation  tends  to 
overpredict  mass  removal  rates,  it  eliminates  a  formidable  task  of  main¬ 
taining  ah  account  of  momentum  and  changes  thereto  throughout  the  flow. 

Determining  the  Surface  Shear  Stresses 

The  surface  shear  forces,  in  conjunction  with  the  viscosity  of 
the  liquid  metal,  determine  the  mass  removal  rate  of  the  molten  mecal. 
The  shear  forces  t  at  the  atmospheric-metallic  interface,  or  "wall"  as 
indicated  by  the  subscript  "w,“  are  determined  by  the  following  equa¬ 
tion  (Ref  8:30}: 

<'5) 

where  p*  is  the  air  density  at  a  reference  enthalpy  to  be  defined  later 
(see  Eq  (19));  Vm  is  the  free  stream  wind  velocity;  and  cf  is  the  local 
skin  friction  coefficient  at  a  distance  x  from  the  leading  edge  of  the 
plate.  The  form  of  the  solution  for  of  is  dependent  upon  whether  the 
aerodynamic  boundary  layer  is  laminar  or  turbulent.  One  set  of  solu¬ 
tions  which  has  been  found  to  give  good  results  (Ref  8:30,35)  is 

cfx  =  0.664A /(4)  (16) 

laminar 

ofx  =  0 . 370/ (logi o  fffl)2,58"  {17) 

‘  turbulent 
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where  the  Reynolds  number  J 

R*  =  xp*Voo  /ii*  (T8) 

and  p*  and  p*  are  the  density  and  dynamic  viscosity,  respectively,  of 
the  air  at  a  reference  enthalpy  i*  (the  asterisk  indicating  the  ref¬ 
erence  condition),  defined  by 

*  - 2 — ~  +  °'22ri  (19) 

The  subscripts  °°  and  w  refer  respectively  to  the  conditions  in  the  free 
stream  beyond  the  aerodynamic  boundary  layer  and  at  the  Wall;  M  is  the 
mach  number,  or  the  velocity  of  the  wind  divided  by  the  speed  of  sound, 
y  is  the  ratio  of  the  specific  heats  at  constant  pressure  c*  and  at 
constant  volume  o*  (for  air  cp/cv  «  1.4);  and  is  a  recovery  factor 
defined  for  laminar  and  turbulent  boundary  flow  conditions  as 


^lair.inar  =v/(P^ 

(20) 

turbulent  =^Pn'> 

(21) 

where  the  Prandtl  number  at  the  reference  enthalpy  conditions  is 

P*  -  op*/k*  (22) 

and  k  is  the  conductivity  of  the  air  at  i*. 

Turbulence  along  flat  plates  generally  occurs  at  a  Reynolds 
number  in  the  range  of  3.5x10s  to  106  (Ref  14:435).  An  arbitrary  value 
of  4><105  for  the  transition  Reynolds  number  was  tentatively  incorporated 
within  the  numerical  model.  Note  tnat  the  Prandtl  number,  used  to  com¬ 
pute  the  recovery  factor  arid  in  turn  the  reference  enthalpy  (Fq  (19)), 
must  also  be  computed  at  the  reference  enthalpy  conditions  (Eq  (22)). 

A  two-step  iteration  ,<as  employed  within  the  model  to  accomplish  this. 
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Equations  (16)  to  (22)  were  in  large  part  empirically  derived. 
Equation  (16),  with  Rfi  replaced  by  ,  is  Blasius'  1908  exact  solution 
for  a  laminar  boundary  layer  assuming  constant  properties,  independent 
of  pressure  and  temperature,  throughout  the  thickness  of  the  boundary 
layer  (Ref  14:128).  Equation  (17),  less  the  asterisk,  is  by  Schulz- 
Grunow  (Ref  8:35).  Modification  of  the  equations  by  computing  various 
properties  at  a  reference  enthalpy  i*  so  as  to  make  Eqs  (16)  and  (17) 
valid  through  hypersonic  mach  numbers  was  devised  by  Eckert  (Refs  8 
and  9).  These  modifications  correct  for  the  variation  of  the  air 
density  and  viscosity  with  changes  in  temperature  as  actually  occurs 
in  boundary  layers  associated  with  high  speed  flow. 


Approximating  the  Convective  Heat 
Transfer  Rate 

Solution  of  the  skin  friction  coefficient  cf  according  to  'iqs 
(16)  or  (17),  above,  facilitates  determination  of  the  convective  neat 
transfer  rate  at  the  same  local  area.  The  appropriate  formula  for  a 
flat  plate  with  a  constant  surface  temperature  (Ref  8:30-31 )  is 


<7 


convection 


%  p*1'- 


(iy  -  Ty) 


(23) 


where  iw  is  the  enthalpy  of  the  air  at  the  wall  temperature  and  iY  is 
the  enthalpy  at  a  recovery  temperature  T  defined  by 


r.lK 


+  T„ 


2a’ 


(24) 


Although  Eq  (23)  has  been  exclusively  used  to  approximate  the  local 
convective  heat  transfer  rate  in  the  model  presented  in  Appendix  A, 
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further  elaboration  has  since  appeared  necessary  to  fully  account  for 
the  effects  of  temperature  variation  along  the  plate  as  would  be  en¬ 
countered  with  localized  laser  heating.  Treatment  of  such  surface 
temperature  variation  may  be  divided  into  that  of  an  essentially 
unheated  starting  length  and  the  temperature  variation  encountered  sub¬ 
sequent  to  that  distance. 

Corrective  factors  for  unheated  starting  lengths  x0  (i.e.,  the 
distance  between  the  plate’s  leading  edge  and  the  point  where  signifi¬ 
cant  temperature  changes  commence)  are  treated  in  Ref  15.  The  factors 
C.F.  given  for  laminar  and  turbulent  flow  (Ref  15:217,242)  to  be  applied 
to  Eq  (23)  are 


C.F. .  . 

laminar 

[1  -  (*o/*)mr1/3 

(25) 

C.F..  V  ,  » 

turbulent 

[1  -  (xo/x)wr1/9 

(26) 

For  a  laser  beam  having  a  10  cm2  cross  section,  directed  at  a  point 
five  feet  from  the  leading  edge  (this  same  distance  was  used  in  all  test 
cases)  and  assuming  a  radial  distance  of  one  inch  from  the  laser  beam 
axis  at  which  the  temperature  rise  begins,  the  turbulent  flow  corrective 
factor  is  computed  to  be 


C.F. 


turbulent 


[1  -  (59/60) 0,9  j  ~ 1/9 


=  1.6 


(27) 


That  is,  the  convective  heat  transfer  rate  should,  under  the  given  situa 
tion,  be  1.6  times  that  approximated  by  Eq  (23).  However,  this  correc¬ 
tion  applies  for  unheated  starting  lengths  which  would  occur  only  if 
the  initial  temperature  r0  of  the  plate  was  at  the  recovery  temperature 
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(see  Eq  (24)).  For  other  initial  temperatures,  Eq  (25)  may  not  be 
entirely  appropriate;  and  additional  investigation  is  considered  neces¬ 
sary  before  adopting  the  factor  as  is  for  cases  in  general. 

To  account  for  the  effect  of  significant  variation  in  surface 
temperature  subsequent  to  modest  variations  along  an  "unheated"  starting 
length  on  the  convective  heat  transfer  rate,  Eckert  (Ref  9:16-21)  offers 
a  numerical  procedure  which  could  prove  adaptable  in  refining  Eq  (23). 
The  compatibility  of  this  correction  procedure  in  conjunction  with  the 
previous  correction  for  unheated  starting  lengths  remains  to  be  deter¬ 
mined  however.  The  change  on  the  convective  heat  transfer  rates  for 
turbulent  flow  due  to  the  local  temperature  variation  alone  (neglecting 
unheated  starting  lengths)  should  not  exceed  ±40%  (Ref  9:18,19).  More 
pronounced  effects  should  occur  with  laminar  flow.  For  the  sake  of 
expedience,  pending  further  development  and  evaluation  of  a  composite 
correction  factor  to  account  for  the  combined  effects  of  relativtly 
unheated  starting  lengths  followed  by  intense  temperature  variations, 

Eq  (23)  was  provisionally  incorporated  within  the  numerical  model.  The 
consequence  of  this  oversimolification  should  not  be  overly  adverse  to 
most  of  the  modeling  results  for  turbulent  flow  (no  laminar  flow  results 
are  presented  in  this  thesis)  since  for  most  cases  the  convective  heat 
transfer  rate  will  be  insignificant  when  compared  to  the  absorbed  flux 
intensity.  However,  for  laminar  flow  and  for  very  low  flux  intensities, 
modification  of  Eq  (23)  appears  to  be  warranted. 

Solution  to  the  convective  heat  transfer  rate,  by  the  use  of  Eq 
(23),  combined  with  Eqs  (11)  and  (12),  the  solutions  to  the  absorbed 
flux  and  radiative  hea'c  transfer  rates,  completes  the  means  of  solving 
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the  rate  of  total  heat  transfer  (reference  Eq  (13))  through  the  exposed 
face  of  any  surface  element  in  the  model's  elemental  array.  Having  thus 
purposefully  digressed  from  this  chapter's  titled  content  in  order  to 
complete  the  modeling  of  the  heat  transfer  proper*  attentioa  is  returned 
to  modeling  the  melt  transfer. 

Modeling  the  Melt  Flow— a  Numerical  Algorithm 

With  a  means  to  calculate  the  shear  stress  x .  at  the  surface  of 

w 

the  molten  metal  (reference  Eq  (15)),  the  mass  flow  within  the  melt 
layer  is  derivable  from  Newton's  equation  (Ref  11:124) 


T 


(28) 


where  y  is  the  dynamic  viscosity  of  the  liquid  metal,  u  is  the  velocity 
of  the  assumed  laminar  flow,  and  y  is  the  normal  distance  from  the 
solid-liquid  interface  within  the  metal  plate  as  depicted  in  Fig.  4. 


■H 
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{§ 

Km 

Fig.  4.  Velocity  Profile  Within  the  Melt  Layer 
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With  an  assumed  steady  flow,  the  shear  stress  remains  essentially  constant 
throughout  the  melt  layer  while  the  dynamic  viscosity  varies  with  the 
temperature  of  the  liquid.  The  velocity  at  any  point  within  the  melt 
layer  can  be  determined  by  solving  Eq  (28)  for  u  as  a  function  of  y. 

w(y)  =  i  I  -1—dy  (29) 

u  4  uQ/) 


A  numerical  approximation  to  the  above  equation,  for  the  velocity  mid¬ 
way  in  an  i,jth  finite  element,  is 


Within  an  increment  of  time  At  the  th  element  is  transferred  down¬ 
stream  a  distance* 

(3,) 


Equivalently,  allowing  the  array  of  spatial  cells  to  remain  stationary 
and  merely  transferring  the  proper  mass  from  one  cell  to  its  adjacent 
downstream  cell,  the  "mass  transfer"  from  the  i,j th  cell  to  the  i- 1,,/th 
cell  is 


TFR.  . 
1,3 


=  *Jt* 

Ax 


Ay. 


1,3 


(32) 


where  TFR.  .  is,  for  computational  purposes,  expressed  as  thickness. 

i'jJ 

This  accountability  of  the  mass  flow  is  illustrated  in  rig.  5.  The 
scheme  portrayed  Is  consistent  with  reality  only  if  the  downstream  ele¬ 
ment  Is  also  molten  or  is  void.  Thus  to  assure  that  the  melt  has  a 
valid  cell  in  which  to  flow,  and  to  simulate  laminar  flow  conditions, 
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Fig.  5.  Graphical  Representation  qf  the  Computational  Algorithm  for 
Transfer  of  Flow  to  Adjacent  Downstream  Finite  Element 

the  model  was  programmed  to  transfer  the  flow  as  hypothetically  exem¬ 
plified  in  Fig.  6.  Note  that  the  array  of  finite  elements  representing 
the  initial  cross  section  of  the  plate  must  be  augmented  with  initially 
empty  cells  to  accommodate  the  melt  flow  and  the  subsequent  lip  formation. 

After  one  time  increment  \t  the  bottom  array  in  Fig.  6,  depending 
upon  the  surface  shear  forces  and  viscosities  within  the  elements,  might 
appear  as  the  top  array  in  Fig.  7.  Observe  that  some  of  the  cells  are 
only  partially  filled  while  others  are  "overflowing."  A  subprogram 
SQUEEZE  (see  Appendix  A)  was  devised  to  redistribute  the  liquid  among 
the  axially  adjacent  elements  in  eaci  array  column  to  eliminate  the 
voids  and/or  overflow  while  ensuring  that  mass  and  thermal  energy  were 
conserved  and  that  the  temperature  gradient  was  valid.  After  execution 
of  subroutine  SQUEF.ZE,  the  hole's  modeled  profile  might  appear  as  in  the 
bottom  array  of  Fig.  7. 

The  uppermost  element  in  each  column  is  allowed  to  range  in  thick¬ 
ness  only  between  ty/2  and  3Ay/2,  where  by  is  the  constant  axial  dimension 
of  all  embedded  ceils.  This  range  ensures  stability  in  the  numerical 
solution  to  the  partial  differential  heat  conduction  equation  without 
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Fig.  6.  A  Numerical  Finite  Element  Model  (Below)  of  a  Hypothesized 
Melt  Flow  (Above) 


undue  degradation  in  modeling  accuracy  (reference  Eq  (7)).  When,  after 
execution  of  subroutine  SQUEEZE,  an  uppermost  element  is  less  than  Ay/2 
thick,  its  energy  and  mass  must  be  combined  with  those  of  its  underlying 
element.  If  the  underlying  element  is  nonmolten,  the  resultant  merged 
element  may  be  molten  or  frozen  dependent  upen  the  combined  mass  and 
thermal  energy  involved.  A  subroutine  FREEZE  (see  Appendix  A),  actuated 
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Fig.  7.  The  Numerically  Modelled  Transfer  of  Melt  Resulting  in  Voids 
and  Excesses  Within  Certain  Cells  (See  Above)  is  Corrected 
(Below)  by  a  Subroutine  SQUEEZE 

within  the  model  for  such  exigencies,  computes  the  resultant  temperatures 
of  such  merged  elements. 

Assembl ing  the  Total  Model 

The  total  model  of  this  study  consists  of  an  iterative  sequencing 
througi  the  various  component  models  which,  as  previously  outlined, 
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approximate  the  several  phenomena  comprising  the  totai  melt- through 
process  of  laser  light  incident  upon  a  metal  plate.  The  sequential 
order  of  the  component  models,  within  each  iteration  of  the  total  model, 
starts  with  the  temperature  and  mass  of  every  finite  element  known, 
Supplementary  matrices  store  tentative  temperature  and  thickness  values 
pending  successful  completion  of  each  iteration.  Property  values  of 
the  elements  at  their  known  temperatures  are  calculated  in  subroutines 
CONDUCT,  SPHEAT,  and  VISCSTY  (the  density  of  the  metal  has  been  taken 
as  invariant  with  temperature).  Heat  transfer  rates  av,  ,  a  _  , 

and  ^radiation  at  exP°ssc*  ^ace  are  comPuted  ’-n  subroutines  FlUXX, 
FRCTION,  and  SCHMIDT.  Heat  balances  for  all  elements  over  the  maximal 
time  increment  At  (reference  £q  (7))  are  also  computed  in  the  subroutine 
SCHMID!  (titled  after  E.  Schmidt,  the  developer  of  a  numerical  and 
graphical  method  for  the  solution  to  one-dimensional  transient  heating 
(Ref  16:294)).  Solid-to-slush  and  slush-to-liquid  phase  changes  are 
searched  for  within  the  cellular  array  in  the  subroutine  XACTIME;  and, 
if  found,  the  subroutine  calculates  the  exact  time  increment  at  which 
the  affected  element  changes  to  slush  or  liquid  as  the  case  may  be  and 
reexecutes  SCHMIDT  with  the  revised  time  increment.  A  subsequent  check 
for  liquid- to-slush  or  s  1  ush-to--so'l id  phase  changes  is  performed  in  sub¬ 
routine  COOIIT,  If  such  a  temperature  drop  in  any  cell  is  detected,  the 
new  temperature  of  the  cell  is  recalculated  using  the  specific  heat  of 
the  cooler  phase  over  that  portion  of  the  time  increment  after  which  the 
phase  change  occurs.  Next,  vaporization,  as  Dreviously  treated,  is 
accounted  for  in  the  subroutine  VAPORIZ.  The  melt  flow,  also  dealt  with 
earlier  and  modeled  in  subroutines  TRNS£rRX ,  SQUEEZE,  and  FREEZE,  concludes 
the  sequence  of  modeled  effects  constituting  the  total  surface  recession 
pnenomencn. 
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V.  Results 

The  essential  result  of  this  thesis  study  is  "QUEST,"  a  model  to 
predict  melt-through  times  and  generate  supplementary  data  concerning 
the  transient  heating  and  surface  recession  dynamics  for  laser  beams 
impinging  upon  flat  metal  plates  in  an  aerodynamic  environment.  The 
model,  computer  programmed  in  FORTRAN  Extended  (a  variant  of  FORTRAN  IV), 
is  presented  in  Appendix  A.  The  version  presented  is  for  a  0.1  cm 
thick  titanium  sheet  at  sea  level  subjected  to  a  mach  0.50  wind  and 
absorbing  a  Gaussian-distributed  energy  flux  of  10  kw  over  a  10  square 
centimeter  spot.  QUEST  was  additionally  exercised  for  aluminum,  mag¬ 
nesium,  and  stainless  steel  and  for  various  thicknesses,  wind  velocities, 
absorbed  laser  powers,  and  spot  sizes.  Numerical  results  from  these 
computer  runs  are  tabulated  in  Appendix  D.  Qualitative  findings  from 
those  results  are  presented  in  the  following  chapter. 

Computer  Program  Specifications 

QUEST  consis*:  of  13  program  subroutines  and  approximately  500 
executable  statements.  Twenty- three  library  subroutines  are  employed, 
including  six  used  in  producing  graphical  plots  of  output  data.  The 
program  requires  less  chan  50,000  octal  v/ord  memory  capacity.  Although 
thus  far  run  exclusively  on  a  CDC  6600  computer,  QUEST  should  be  adaptable 
to  other  digital  computers,  with  perhaps  the  exception  of  tne  CALC0MP 
plotting  feature.  However,  the  data  that  are  plotted  nre  also  printed, 
ano  thus  this  feature  could  be  eliminated. 
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Input,  which  includes  tnach  number,. altitude,  distance  from  the 
leading  edge  of  the  plate  to  the  laser  spot,  initial  temperature  of  the 
metal,  and  the  total  rate  of  energy  absorption,  is  by  data  cards.  Other 
parameters  peculiar  to  each  metal  are  incorporated  directly  within  the 
appropriate  program  subroutines. 

Execution  time  for  the.  program  varies  with  the  input  parameters, 
‘.nicker  metals,  lower  energy  flux  intensities,  and  finer  stratification 
of  the  plate's  cellular  array  all  require  additional  computation  time. 
Total  run  times  have  ranged  from  less  than  10  to  over  500  seconds.  The 
example  presented  in  Appendix  A  compiled  in  17  seconds  and  required  an 
additional  253  seconds  for  execution.  Conversion  of  the  source  deck 
from  the  FORTRAN  IV  programming  language  to  binary  machine  language 
eliminates  the  compilation  time. 

Output,  as  programmed  and  an  example  of  which  may  be  found  in 
Appendix  A.,  consists  of  an  abbreviated  three  line  diagnostic  printout 
after  each  occurrence  of  a  phase  change  within  the  finite  element  array 
plus  an  extended  printout  of  analytic  data  at  the  conclusion  of  melt- 
through.  In  addition  to  the  melt- through  time,  the  final  printout 
provides  the  input  parameters;  selected  thermodynamic  property  values 
of  the  metal;  the  final  temperature  array;  comparative  values  of  the 
incident  (absorbed)  laser  energy,  convective  and  radiative  heat  trans¬ 
ferred,  and  energy  consumed  by  vaporization  along  the  spot-bisecting, 
wind  line;  the  remaining  total  and  non^molten  thicknesses  of  the  plate 
along  the  same  diagnostic  line;  the  temperature  history  of  the  plate's 
back  surface  opposite  the  center  of  the  laser  beam  spot;  the  position 
o*  solid-liquid  interface  at  the  center  of  the  spot  versus  time;  and 
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the  surface  recession  rate  history  at  the  center  of  the  hole.  Addi¬ 
tional  output  inform  ttion,  if  needed,  may  be  extracted  with  additional 
programming. 

Model  Validation 

The  numerical  results  from  any  complex  simulation  must  be  suspect 
until  proven  true.  Unfortunately,  if  the  results  could  be  fully  vali¬ 
dated,  there  probably  would  have  been  no  need  for  a  model  in  the  first 
place.  Such  is  the  case  with  QUEST— only  limited  means  exist  to  even 
partially  substantiate  the  results.  Despite  these  inherent  obstacles, 
great  care  was  taken  in  programming  QUEST  to  attain  what  is  believed  to 
be  mathematical  and  logical  integrity.  QUEST  was  additionally  tested 
fo"  sensitivity  and  results  were  compared  with  theoretical  values  and 
experimental  findings. 

Model  Sensitivity.  The  sensitivity  of  the  model  QUEST  to  various 
axial  dimensions  of  the  diagnostic  cellular  array  was  tested  in  order 
to  find  an  acceptable  trade-off  between  the  accuracy  of  the  results  and 
required  computer  processing  time.  The  combined  test  results  for 
aluminum,  0.05,  0.10,  and  0.20  :m  thick  at  mach  numbers  of  0.0,  0.1, 
and  0.5  and  for  10  kw  of  power  over  a  1  cm2  spot,  are  shown  in  Fig.  8. 

Deviations  of  melt-through  times  from  those  predicted  with  a 
20-tiered  array  were  less  than  ]%  for  reductions  in  the  axial  array 
dimension  n  to  10.  With  n  reduced  to  5,  however,  the  deviation  increased 
to  about  5%.  In  making  trade-offs  between  accuracy  and  computer  run 
times,  a  rough  guide  to  the  computations—and  the  computer  time  required— 
is  Eq  (7);  halving  the  number  of  tiers  will  approximately  quadruple  the 
incremental  time  in  each  forward  time  step  and  thus  reduce  by  four 
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Fig.  8.  Sensitivity  of  Melt-Through  Times  Predicted  by  QUEST  to 
Changes  in  n,  the  Number  of  Tiers  in  the  Analytic  Array 
of  Finite  Elements 


times  the  computer  processing  time  required.  The  results  presented  in 
Appendix  D  were  obtained  with  10-tiered  arrays  (exclusive  of  the  aug¬ 
mentation  required  for  mass  flow  and  lip  formation)  for  1  cm2  incident 
beam  spots  and  with  5-tiered  arrays  for  10  cm2  spots.  The  array  was 
limited  to  5  tiers  in  these  latter  cases  to  reduce  computer  run  times. 

The  computer-programmed  model  QUEST  is  dimensioned  for  a  30- 
tiered  array.  Approximately  one-third  of  these  tiers,  however,  must  be 
reserved  for  modeling  the  lip  formation,  thus  allowing  up  to  20  tiers 
for  the  initial  plate  thickness.  If  less  accuracy,  both  in  modeling 
the  hole  and  lip  profile  and  in  predicting  the  melt-through  time,  is 
permissible  and  the  array's  axial  dimension  were  accordingly  halved, 
the  central  computer  memory  required  by  QUEST  could  be  reduced  by 
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approximately  6000  octal  words.  If,  in  addition,  the  CALCOMP  plotting 
feature  were  removed  and  the  FORTRAN  programming  wcr;  teplaced  with 
binary  machine  language,  the  total  memory  core  required  by  QUEST  reduces 
to  less  than  30,000  octal  words— a  significant  reduction  in  computer 
requirements  with  only  limited  restriction  in  achievable  accuracies  of 
predicted  melt- through  times. 

Comparison  of  Predicted  Results  with  Known  Solutions.  A  complete 
absence  of  known  solutions  to  the  total  phenomenon  of  melt- through  in 
an  aerodynamic  environment  and  a  dearth  of  exact  solutions  to  even 
fictitiously  simplified  melt- through  processes  restrict  the  ability  to 
validate  the  results  of  QUEST.  The  most  satisfactory  of  the  simplified 
known  solutions  found  to  gauge  the  soundness  of  QUEST  is  that  of  melt- 
through  with  no  wind  but  with  "instantaneously-vanishing"  melt  and 
one-dimensional  (axial)  heat  conduction.  Solutions  to  melt-through 
times  under  such  simplifying  assumptions  reduce  to 

t  =  —  (33) 

melt-through  ?o(mol  wt) 

where 

M  =  the  difference  in  the  heat  content  between  the  melting 
point  and  the  initial  temperature  of  the  metal  (calories/ 
gm-mole) , 

p  =  density  (gm/cm3), 
y  =  thickness  of  plate, 

qo  =  absorbed  flux  intensity  at  mid-beam  (see  Eq  (4)),  and 
mol  wt  =  molecular  weight  of  the  metal. 

Heat  contents  for  aluminui1,  magnesium,  and  titanium  were  computed  from 
applicable  formulas  presented  by  Wicks  and  Block  in  Ref  17  (see  Appendix 
C).  The  values  for  stainless  steel  w«re  approximated  by  those  found  for 
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iron.  Tne  melt-through  times  according  to  Eq  (33)  for  the  metals,  spot 
sizes,  and  laser  powers  presented  in  Appendix  D  are  given  In  Table  I. 

For  all  but  near-hypersonic  conditions,  the  tabulated  melt-through 
times  in  Table  I  are  less  than  the  corresponding  times  modeled  by  QUEST 
and  presented  in  Appendix  D.  Times  from  QUEST  ranged  to  3QQ*  of  the 
corresponding  times  arrived  at  from  Eq  (33) — the  maximum  difference 
occurring  for  0.20  cm  thick  titanium  absorbing  1  kw  power  over  a  10  cm2 
spot  subjected  to  mach  2  airflow.  The  significant  difference  in  this 
cited  case  is  primarily  attributable  to  aerodynamic  cooling.  The  con¬ 
vective  heat  losses,  as  computed  by  QUEST,  were  over  half  (56%)  of  the 
laser  energy  absorbed  (which  alone  would  more  than  double  the  melt- 
through  time  as  gauged  by  Eq  (33).  Additional  heat  loss  by  radiation 
(5%  of  that  absorbed)  and  by  radial  heat  conduction  combined  with  the 
heating  of  the  liquified  metal  account  for  the  balance  of  the  additional 
Incident  power  required  and  the  consequential  longer  melt-through  times. 
Titanium,  as  opposed  to  the  other  metals  modeled,  exhibits  greater  con¬ 
vective  heat  losses  due  to  its  high  melting  temperature.  Its  melting 
point  exceeds  the  recovery  temperature  within  the  boundary  layer  (refer¬ 
ence  Eqs  (23)  and  (24))  at  sea  level  wind  velocities  of  up  to  mach  4.0. 

It  must  be  pointed  out  that  for  this  particular  case,  if  corrective 
factors  for  unheated  starting  lengths  and  variant  surface  temperatures 
had  been  applied  to  the  convective  cooling  as  nresently  approximated  by 
Eq  (23),  melt-through  could  possibly  have  been  precluded  altogether. 

This  would  happen  if  the  corrections  increased  the  convective  cooling 
to  a  point  where  the  absorbed  flux  and  the  combined  heat  losses  were 
equal.  Thus,  as  previously  recommended,  further  refinement  to  Eq  (23) 
appears  to  be  in  order  for  modeling  cases  involving  very  low  flux  intensities. 
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TABLE  I 

Melt-Through  Times  (in  Seconds)  of  Metal  Plates  Assuming 
Instantaneously  Vanishing  Melt,  One-Dimensional 
Heat  Conduction,  and  Ho  Wind 
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Comparison  of  Predicted  and  Experimental  Results.  An  experi¬ 
mental  result  as  reported  In  Ref  10  was  compared  with  the  predictions 
of  QUEST.  The  experiment  involved  a  16  mil  sheet  of  stainless  steel 
irrauiated  over  a  4.676  cm2  spot  by  a  9  kw  laser  beam.  Melt-through 
was  judged  from  high-speed  movie  photography  to  occur  at  0.40  seconds. 
The  thermal  history  of  the  rear  surface  was  recorded  by  means  of  a 
transducer  aligned  with  the  beam  axis  and  affixed  to  the  back  of  the 
plate. 

Corresponding  melt-through  times  predicted  from  QUEST,  assuming 
constant  absorptances  of  0.2  and  0.3  with  no  wind  (and  thus  no  melt 
removal),  were  0.51  and  0.34  seconds,  respectively.  The  rear  surface 
temperature  histories  as  experimentally  and  predictively  obtained  are 
compared  in  Fig.  9.  Also  Included  are  rear  surface  temperatures  as 
predicted  by  Torvik  (Ref  10:30)  assuming  identical  absorptances  but 
with  instantaneous  melt  removal.  As  pointed  out  by  Torvik  (Ref  10: 
31-32),  the  comparative  results  indicate  that  the  absorptance  of  stain¬ 
less  steel  must  increase  significantly  during  the  heating  to  account  for 
the  nonlinearity  in  the  experimentally  obtained  temperature  history. 

Differences  also  appear  between  the  rear  surface  temperature 
histories  as  predicted  by  Torvik  and  by  QUEST.  The  differences  are 
primarily  attributed  to  the  fact  that  QUEST  utilizes  temperature- 
dependent  property  values  as  opposed  to  constant  properties  assumed  in 
the  former  model.  In  this  case  of  stainless  steel  (AISI  304),  QUEST 
employs  a  specific  heat  varying  from  0.106  cal/gm/°K  at  300°K  to 
0.166  cal/gm/°K  at  just  below  the  melting  temperature,  while  Torvik 
(Ref  10:31)  assumed  a  constant  specific  heat  of  0.100  cal/gm/°K  over 
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Fig.  9.  Rear  Surface  Temperatures  of  a  16  mil  Stainless  Steel 
(AISI  304)  Plate  as  Determined  from  Experiment  Results 
and  as  Predicted  by  Modeling 


the  entire  temperature  range.  (Equations  used  in  QUEST  for  thermodynamic 
values  of  metals  as  a  function  of  temperature  appear  in  Appendix  C.) 

The  validity  of  QUEST  cannot  be  ascertained  with  this  single 
comparison  to  experimental  results  due  to  the  necessity  of  assuming  an 
absorptivity  of  photons  by  the  metal.  The  absorptivity  itself  is  a 
phenomenon  requiring  temperature-dependent  modeling.  Despite  being 
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limited  at  present  to  constant  absorptivlties*  however,  QUEST  does 
predict. melt-through  times  of  the  same  order  as  the  experimentally 
obtained  results.  Such  results,  coupled  with  the  comparisons  to  known 
solutions  in  the  previous  section,  are  at  least  indicative  of  the 
model's  credibility. 
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VI.  Findings 

The  limited  numerical  results  thus  far  computed  by  QUEST  (see 
Appendix  D)  have  been  evaluated  in  part  to  determine  qualitatively  the 
general  effects  of  wind  on  melt-through  and  the  influence  of  other 
parameters  on  these  effects.  These  findings  must  be  considered  as 
tentative,  pending  further  validation  of  the  model,  and  should  be 
extrapolated  to  cases  not  yet  tested  with  extreme  caution. 

Aerodynamic  Effects  on  Melt-Through  Times 

The  melt- through  times  computed  by  QUEST  for  the  metals,  thick¬ 
nesses,  wind  velocities,  energy  fluxes,  and  beam  sizes  as  tabulated  in 
Appendix  D  are  graphically  presented  in  normalized  form  in  Figs.  10 
through  13  with  mach  number  being  the  Independent  variable.  The  graphs 
reveal  no  simple  relationship  between  melt-through  time  and  mach 
number;  that  Is,  the  overall  effect  of  air  flow  on  the  melting  process 
must  in  turn  be  a  function  of  the  remaining  parameters.  For  instance, 
with  spot  size  and  energy  flux  unspecified,  one  cannot  even  predict 
whether  wind  will  retard  or  speed  melt-through.  To  wit,  0.2  cm  thick 
titanium  in  a  mach  2.0  airstream  required  from  over  twice  (for  1  kw 
Gaussian  absorption  over  a  10  cm2  spot)  to  about  one-fourth,  (for  10  kw/ 

1  cm2  spot)  the  time  to  melt-through  when  compared  to  the  same  metal 
sheet  subjected  to  the  same  conditions  in  still  air  (see  Fig.  12).  A 
further  indication  of  the  complicacy  involved  in  analyzing  the  aero¬ 
dynamic  effects  is  found  with  0.2  cm  thick  titanium  and  stainless  steel 
absorbing  1  kw/1  cm2.  In  these  two  cases,  instead  of  a  monotcnic  decrease 
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in  melt-through  time  with  increase  in  wind  velocity,  the  time  to  melt- 
through  actually  increased  as  the  mach  number  increased  from  1.0  to  2.0. 
Other  apparent  anomalies  occurred  for  0.05  cm  thick  sheets  absorbing 
1  kw/1  cm2  and  10  kw/10  cm2  where  a  slight  increase  in  melt-through 
times  was  often  noted  at  low. subsonic  velocities  followed  by  the  other¬ 
wise  characteristic  decrease  in  time  to  melt- through.  All  of  these 
apparent  anomalies  can  be  rationalized  by  reconsidering  the  two  basic 
mechanisms  by  which  wind  affects  laser-induced  melting. 

Aerodynamic  Heating/Cooling  Effects-  The  first  mechanism,  aero¬ 
dynamic  heating  or  cooling  by  convective  heat  transfer,  either  hastens 
or  retards  melt- through  by  respectively  augmenting  or  detracting  from 
the  laser-absorbed  energy.  Which  effect  occurs  is  dependent  upon  the 
relative  surface  and  aerodynamic  boundary  layer  recovery  temperatures 
and  their  respective  enthalpies  (reference  Eq  (23)).  The  dependence  of 
the  recovery  temperature  upon  the  square  of  the  free  airstream  velocity 
(Eq  (24))  plus  the  strong  dependence  of  the  magnitude  of  the  convective 
heat  transfer  on  an  additional  power  of  the  free  stream  velocity  (Eq 
(23))  forces  the  effect  of  this  first  mechanism  to  vary  significantly 
with  mach  number.  At  low  mach  numbers  aerodynamic  cooling  is  to  be 
expected  within  the  heated  spot;  but  because  of  the  relatively  slow 
free  airstream  velocity  the  magnitude  of  such  cooling  will  be  limited. 
Often,  because  of  the  limited  cooling,  the  retardation  effect  on  melt- 
through  times  was  found  to  be  masked  by  a  stronger  melt  transfer  effect. 
One  set  of  cases  where  such  masking  did  not  occur  was  for  1  kw  of  power 
absorption  over  10  cm2  spots.  In  these  cases,  for  all  metals  (see  lower, 
left  graph  in  Figs.  10-13),  aerodynamic  cooling  exceeded  all  other  effects 
for  mach  numbers  up  to  and,  in  the  cases  of  titan'um  and  stainless  steel, 
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In  excess  of  2.5.  Thereafter,  Increases  in  aerodynamic  heating  with 
increased  wind  velocities  accelerated  melt- through  to  times  approaching 
or  less  than  melt-through  times  in  still  air.  Other  «ases  where  aero¬ 
dynamic  heating  effects  are  evident  are  for  1  kw  absorption  over  1  cm2 
spots  (upper  left  graphs  in  Figs.  10-13).  The  decrease  in  melt-through 
times  observed  at  mach  numbers  greater  than  3.0  is  attributed,  for  these 
particular  cases,  nearly  exclusively  to  aerodynamic  heating.  Such  aero¬ 
dynamic  heating  was  found  possible  to  reduce  melt-through  times  to  even 
less  than  those  theoretically  possible  with  one-dimensional  (axial)  heat 
flow  and  instantaneous  melt  removal  (viz  0.1  and  0.2  cm  thick  Mg  exposed 
to  1  kw  over  10  cm2  at  mach  4.0— reference  Table  I  and  Appendix  D). 

Melt  Transfer  Effects.  Mass  flow,  or  "melt  transfer,"  the  secord 
mechanism  by  which  wind  affects  melt-through  times,  generally  hastens 
melt- through  by  removing  a  portion  of  the  melt  over  the  hottest  portion 
of  the  spot  and  thereby  eliminating  it  from  further  heating.  The  energy 
thus  "saved"  is  absorbed  by  the  underlying  elements  which  heat  and  con¬ 
sequently  melt  faster,  thus  perpetuating  the  accelerated  melting  process. 
The  effect  was  most  evident  in  the  thicker  metals  (as  can  be  seen  in  all 
but  the  lower  left  graph  in  Figs.  10-13).  This  characteristic  is  explain¬ 
able  due  to  the  fact  that  for  a  given  laser  beam  in  still  air,  thicker 
plates  either  attain  higher  front  face  and  Interior  temperatures  or  con¬ 
sume  a  greater  portion  of  the  incident  energy  in  vaporization  than  do 
thinner  sheets.  Stated  differently,  in  the  absence  of  wind,  thicker 
sheets  will  at  melt-through  exhibit  a  higher  average  heat  content  per 
unit  mass  than  will  a  thinner  sheet  exposed  to  an  Identical  laser  beam. 
Prompt  removal  of  melt  from  thicker  metals  at  or  close  to  the  liquidus 
temperature  thus  has  a  proportionally  greater  effect  in  reducing 
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NORMALIZED  MELT-THROUGH 


Fig.  10.  Dimensionless  Melt-Through  Times  (Normalized  by  Dividing 
by  the  Melt-Through  Time  at  0.0  Mach)  for  0.05,  0.10,  and 
0.20  cm  Thick  Aluminum  Plates  Absorbing  1  and  10  kw  of 
Light  over  1  and  10  cm2  Spots  at  Mach  Numbers  from  0  to  4 
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Fig.  11.  Dimensionless  Melt-Through  Times  (Normalized  by  Dividing 
by  the  Melt-Through  Time  at  0.0  Mach)  for  0.05,  0,10,  and 
0.20  cm  Thick  Magnesium  Plates  Absorbing  1  and  10  kw  of 
Light  over  1  and  10  cm2  Spots  at  Mach  Numbers  from  0  to  4 
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Fig.  12.  Dimensionless  Melt-Through  Times  (Normalized  by  Dividing 
by  the  Melt-Through  Time  at  0.0  Mach)  for  0.05,  0.10,  and 
0.20  cm  Thick  Titanium  Plates  Absorbing  1  and  10  kw  of 
Light  Over  1  and  10  cm2  Spots  at  Mach  Numbers  from  0  to  4 


NORMALIZED  MELT-THROUGH  TIME  t/t 
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Fig,  13.  Dimensionless  Melt-Through  Times  (Normalized  by  Dividing 
by  the  Melt-Through  Time  at  0.0  Mach)  for  0.05,  O.TO,  and 
0.2G  cm  Thick  Stainless  Steel  Plates  Absorbing  1  and  10  kw 
of  Light  Over  1  and  10  cm2  Spots  at  Mach  Numbers  from  0  to  4 
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melt- through  times  than  with  similar  melt  removal  from  thinner  sheets 
identically  radiated.  Another  characteristic  of  the  mass  flow  mechanism 
is  that  its  effect  quickly  increases  with  increased  wind  velocity  until 
it  approaches  a  maximum.  The  theoretical  maximum  effect  would  be  instan¬ 
taneous  melt  removal;  but,  for  the  spot  sizes  tested,  this  maximum  was 
virtually  achieved  as  early  as  sea  level  mach  0.5  (in  the  cases  of 
extremely  high  flux  intensities,  i.e.,  10  kw/1  cm2,  the  maximum  was 
approached  at  about  mach  1.0).  Further  increases  in  wind  velocities 
did  not  appreciably  hasten  melt-through  due  to  the  effect  of  mass  flow. 
The  changes  that  are  noted  in  melt- through  times  at  higher  mach -numbers 
are  attributed  instead  to  aerodynamic  heating/cooling  effects  which  are, 
as  previously  stated,  strongly  a  function  of  wind  velocity  particularly 
at  the  higher  mach  numbers. 

The  overall  effect  of  wind  on  melt- through  is  simply  the  summation 
of  the  two  individual  effects  noted  above.  But,  since  the  individual 
effects  can  vary  considerably  with  changes  in  the  metal,  the  metal 
thickness,  the  flux,  and  the  beam  or  spot  size,  these  other  parameters 
must  also  be  examined  as  to  their  influence  on  melt- through  and  the 
aerodynamic  effects. 

The  Influence  of  the  Metal  on 
Aerodynamic  Effects 

The  corresponding  melt- through  time  versus  mach  number  curves  for 
the  four  metals  in  Figs.  10-13  exhibit  definite  similarities  in  form  and 
in  some  of  the  cases  in  magnitude  as  well.  These  similarities  indicate 
that  the  metals  have  similar  influences,  qualitatively,  on  the  effects 
of  sir  flow  on  laser-induced  melting.  The  curves  of  Figs.  10-13  in  their 
normalized  form  do,  however,  discount  quantitative  differences  between 
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the  metals.  Therefore,  the  effect  of  the  metal  on  "aerostatic"  {meaning 
in  still  air)  melt-through  times  is  necessary  to  appreciate  the  magnitude 
of  those  differences  that  do  exist  in  the  influence  of  the  metal  on 
aerodynamic  effects. 

Briefly,  the  relatives  times  for  aerostatic  melt-through  for 
identical  thicknesses  and  energy  absorption,  taking  the  value  for  magnesium 
as  1.0,  were  found  to  be  approximately  1.5  for  aluminum,  3.5  for  titanium, 
and  4.5  for  stainless  steel.  Interestingly,  the  densities  of  the  same 
respective  Metals,  relative  to  magnesium,  are  1.5,  2.5,  and  4.5.  The 
correspondence  of  the  relative  densities  to  relative  melt-through  times 
indicates  that  it  is  the  difference  in  density  that  largely  accounts  for 
the  variation  between  these  four  metals  in  times  to  melt  through  in  still 
air.  Titanium,  the  only  exception  to  the  above  curiousity,  requires 
additional  time  for  malt- through  due  to  its  high  melting  point  (1988°K  cf 
923°K  for  Mg).  Stainless  steel  also  has  a  high  melting  point  (approximately 
1730°K),  but  a  low  specific  heat  (»0.166  cal/gm/°K  cf  «0.312  cal/gm/°K 
for  Mg)  keeps  the  heat  content  per  unit  mass  at  melting  comparable  to 
that  of  magnesium.  Other  properties  of  metals  that  can  influence  aero¬ 
static  melt- through  times,  but  which  in  the  case  of  the  above  metals 
apparently  "balanced  out"  in  the  net  effect,  are  the  heat  of  fusion,  the 
boiling  point,  the  heat  of  vaporize t ion,  and  conductivity.  Viscosity, 
in  addition  to  several  of  the  above  properties,  can  influence  aerodynamic 
effects  on  melt- through  times. 

Aerodynamic  effects  on  melt-through  were  observed  to  be  most  pro¬ 
nounced  on  titanium  and  stainless  steel  (see  Figs.  10-13).  Lesser  effects 
were  noted  on  magnesium  and  aluminum.  The  causes  of  these  differences 
between  the  four  metals  were  traced  to  three  properties:  the  respective 
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conductivities,  specific  heats,  and  melting  temperatures  of  the  metals. 

The  conductivity  and  specific  heat  were  found  to  Indirectly  affect  the 
mass  flow  effect  while  the  melting  point  was  determined  to  directly 
influence  aerodynamic  heating/cooling. 

The  Influence  of  Metallic  Properties  on  MeU  Flow  Effects.  In 
the  case  of  mass  flow,  poor  conductivity  forces  an  Increased  temperature 
gradient  and  consequently  higher  average  heat  content  in  any  melt  not 
removed  by  aerodynamic  frictional  forces.  As  previously  explained,  the 
faster  the  removal  of  such  heat-cumulating  material,  the  faster  the 
melt- through.  The  most  heat-cumulative  of  the  molten  metals  are  those 
with  the  poorest  conductivities  and  the  highest  specific  heats.  Liquid 
titanium,  stainless  steel,  and  magnesium  hwe  approximately  0.15,  0.25, 
and  0.75  the  conductivity  of  molten  aluminum,  and  the  corresponding 
relative  axial  temperature  gradients,  disregarding  radial  heat  conduc¬ 
tion,  will  be  approximately  proportional  to  the  reciprocals  of  the 
relative  conductivities.  Temperature  gradients  at  time  of  melt-through 
for  0.05,  0.10,  and  0.20  cm  thick  titanium  and  aluminum  as  modeled  by 
QUEST  for  1  kw  power  absorption  over  a  1  cm2  spot  (see  Fig.  14)  cleu.Ty 
demonstrate  the  significant  differences  possible  in  temperature  gradients 
between  metals.  The  products  of  the  relative  temperature  gradients  and 
the  relative  specific  heats  (see  Appendix  C  for  property  values  of  the 
metals)  should  provide  an  index  of  the  potential  amount  of  heat  acquirable 
by  unremoved  melt  M '  x  for  a  given  flux  intensity  and  plate  thickness 
(low  boiling  points  combined  with  small  heat  of  vaporization  in  the 
case  of  relatively  thick  metals  could  r  duce  the  accuracy  of  such  an 
index).  The  resultant  products  are  Ti,  6.8;  stainless  steel,  4.3;  Mg, 

2.6;  and  A1 ,  1.6. 
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Fig.  14.  Interior  Temperatures  of  0.05,  0.10,  and  0.20  cm  Thick 
Aluminum  and  Titanium  Plates  ar.  Time  of  Melt-Through 
After  Exposure  to  1  kw  Radiation  Over  a  1  cm2  Spot  as 
Modeled  by  QUEST 
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Another  property  of  metal,  viscosity,  intuitively  should  deter¬ 
mine  the  mach  number  at  which  melt  flow  achieves  a  given  percentage  of 
its  ultimate  potential  effect.  However,  despite  titanium  and  stainless 
steel  having  viscosities  at  the  melting  temperature  about  four  times 
as  great  as  magnesium  and  aluminum,  no  significant  shift  between  metals 
was  observable  from  Figs;  10-13  in  the  mach  number  at  which  melt  flew 
produced  similar  effects.  It  appears,  from  the  limited  data  available, 
that  mass  flow  in  general  approached  its  maximum  effect  between  mach 
0.5  and  1.0.  However,  it  is  anticipated  that  further  sensitivity 
testing  in  the  subsonic  regime  will  reveal  that  viscosity  does  have 
some  discernible  but  limited  influence  in  the  mach  number-mass  flow 
effect  relationship. 

The  Influence  of  Metallic  Properties  on  Aerodynamic  Cooling/ 

Heating  Effects.  As  indicated  previously,  the  melting  point  is  the 
primary  determinant  of  a  metal's  influence  on  the  convective  heat  transfer 
effects.  From  Eq  (23),  it  is  the  difference  between  the  aerodynamic 
boundary  layer  recovery  enthalpy  and  the  surface  enthalpy  (or  the 
difference  between  their  related  temperatures)  that  not  only  determines 
whether  cooling  or  heating  occurs  but  also  affects  the  magnitude  of 
whichever  effect  takes  place.  Since  the  melting  point  serves  as  a 
lower  limit  to  the  surface  temperature  at  the  spot's  center  during  the 
latter,  and  in  most  cases  the  major  portion  of  the  melt-through  process, 
this  property  is  strongly  influential  in  the  aerodynamic  heating  or 
cooling  effect.  Titanium  and  stainless  steel,  with  high  melting  points 
(approximately  1988°K  and  1730°K,  respectively),  demonstrate  their 
influence  on  aerodynamic  cooling  in  the  lower  left  graphs  of  Figs. 

12-13.  Compared  to  the  corresponding  graphs  for  aluminum  and  magnesium 


56 


GAW/MC/72-17 


(Figs.  10-11)  the  magnitude  of  the  cooling  effect  on  delaying  me«t- 
„  through  for  titanium  and  stainless  steel  is  significantly  amplified 
and  peaks  at  nearly  twice  the  mach  number.  The  elevated  melting  tem¬ 
peratures  of  titanium  and  stainless  steel  also  manifest  their  influence 
on  aerodynamic  cooling  in  the  cases  of  0.2  cm  thick  sheets  exposed  to 
1  kw  over  1  cm2  where  melt- through  times  actually  increased  from  mach 
1.0  to  mach  2.0.  No  similar  effect  was  noted  with  cooler  melting 
(933°K  and  923°K)  aluminum  and  magnesium.  Conversely,  aerodynamic 
heating  effects  are  enhanced  by  the  lower  melting  temperatures  of  alumi¬ 
num  and  magnesium.  The  enhancement  is  manifested  by  faster  melt-through 
at  reduced  free-airstream  velocities.  Such  is  evidenced  in  the  cases 
of  1  kw/1  cm2  and  10  kw/10  cm2  power  absorption  per  spot  size  (upper 
left  and  lower  right  graphs  of  Figs.  10-11)  where  appreciable  reduc¬ 
tions  In  melt-through  times  commenced  as  early  as  mach  2.0  for  these 
two  metals  while  little  or  no  reduction  occurred  in  stainless  steel  and 
titanium  until  after  mach  3.0. 

The  Influence  uf  Plate  Thickness 
on  Melt-Through  Times 

In  Ref  12  Torvik  determined  that  the  ratio  of  total  power  to 
plate  thickness  was  the  criterion  that  decided  the  accuracy  of  a  one¬ 
dimensional  (axial)  approximation  to  interior  heat  flow  as  opposed  to 
three-dimensional  reality.  That  is,  for  a  given  power,  it  is  the  thickness, 
with  the  spot  size  having  essentially  no  influence,  that  determines  the 
proportion  of  radial  heat  loss  to  the  heat  retained  directly  beneath 
the  laser  spot.  These  findings  were  based  on  instantaneous  melt  removal. 
Ho'ever,  another  influence  of  plate  thickness,  in  addition  to  its  effect 
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on  radial  heat  loss,  must  be  accounted  for  in  non-instantaneous  melt 
removal  cases.  This  additional  factor  is  the  heat  content  within  the 
molten  metal.  As  can  be  seen  from  Fig.  14,  the  heat  content  in  the 
melt  vor  a  particular  metal  in  the  aerostatic  case  approximately  quad¬ 
ruples  in  magnitude  with  a  doubling  of  thickness  (this  can  be  viewed 
as  doubling  the  base  and  nearly  doubling  the  height  of  the  approximate 
triangles  formed  from  the  graph's  axes  and  the  interior  temperatures 
within  the  molten  slab--the  area  of  which  is  proportional  to  the  added 
heat  content  of  the  metal  above  the  melting  point).  This  nonlinear  heat 
requirement,  combined  with  similarly  nonproportional  heat  requirements 
for  radial  heat  conduction  with  increases  in  thickness,  accounts  for 
the  increase  in  melt-through  time  per  unit  thickness  with  increased 
plate  thickness  as  calculated  from  the  results  tabulated  in  Appendix  D 
and  as  presented  in  Fig.  16.  The  ratio  of  aerostatic  melt-through 
times  per  cm  of  thickness  (tl/cm)/(t2/cm) ,  due  to  the  nonlinear  heating 
requirements  of  the  melt  for  different  plate  thicknesses  should  be 
approximately 


tz  /cin 


t\  /cm 


[a//  ..  +  (H  -H  .  jl 

L  melt2  '  m.p.  start 'j 

A//  ,  +  (//  -H  . 

meltj  \  m.p.  start/J 


(34) 


where  M  ,  refers  to  the  heat  acquired  per  unit  mass  of  melt  in  excess 
melt 

to  the  heat  content  of  the  same  metal  at  the  melting  point  /;  ,  and 

m.p, 

H  _  is  the  heat  content  of  the  metal  prior  to  heating.  For  the  case 
start  r  3 

of  titanium  plates  0.05  and  0.2  cm  in  thickness,  the  above  equation  becomes 


t(0.2  cm  Ti) /cm  ^  f(1150°K/2) (0,167  cal/gm°K)  +  (341  cal/gm)] 
C (0.05  cm  Ti)/cm  [ (400°K/2) (0. 167  cal/gm°K)  +  (341  cal/gm) 


=»  1.165 


(35) 
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Fig.  15.  Time  to  Melt-Through  per  Unit  Thickness  of  Stainless 
Steel,  Titanium,  Aluminum,  and  Magnesium  as  a  Function 
of  Plate  Thickness 
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where  and  ^meltz  were  approximated  from  the  temperature  pro¬ 

files  in  Fig.  14  in  conjunction  with  the  specific  heat  for  titanium 
(see  Appendix  C).  The  actual  predicted  ratio,  from  Fig.  15,  is  1.38. 
The  difference  is  attributed  to  the  additional  retardation  effect  of 
increased  radial  heat  conduction  with  increased  plate  thickness.  If 
such  is  the  case,  the  retardation  from  radial  heat  loss  as  modeled  by 
QUEST  and  from  melt  heating  are  roughly  the  same  for  the  aerostatic 
melt-through  case.  However,  it  mus'  be  remembered  that  QUEST  is  only 
a  two-dimensional,  model  and  thus  does  not  account  for  all  radial  heat 
loss;  and  that,  therefore,  radial  heat  flow  should  have  an  even  greater 
effect  than  melt  heating. 

It  was  noted  that  extrapolation  of  the  aerostatic  melt-through 
time  per  thickness  curves  in  Fig.  15  to  the  point  of  zero  thickness 
nearly  coincided  with  the  corresponding  values  computed  for  one¬ 
dimensional  heat  flow  and  instantaneous  melt  removal  from  Table  I. 

This  seems  reasonable,  but  similar  extrapolation  of  the  aerodynamic 
(for  mach  equal  to  1.0)  curves  gives  values  indicating  longer  melt- 
through  times  per  unit  thickness.  A  possible  explanation  for  this 
curiosity  would  be  a  retardation  effect  on  melt- through  due  to  aero¬ 
dynamic  cooling  as  was  noted  at  lower  mach  numbers  for  0.05  cm  thick 
plates  absorbing  1  kw  of  power  over  1  cm*  (see  Figs.  10-13). 

The  aerodynamic  curves  of  Fig.  15  also  have  less  slope,  repre¬ 
senting  the  fact  that  with  increasing  plate  thickness  the  wind  negates, 
through  melt  removal,  more  and  mere  of  the  potential  retardation  effect 
of  heat  accumulation  within  the  melt  layer.  The  remaining  slope  in 
the  aerodynamic  curves  is  mostly  attributed  to  radial  heat  losses  since 
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mach  1.0,  sea  level  air  *low  should  provide  near-instantaneous  melt 
removal  over  a  1  cm2  spot. 

The  Influence  of  Flux  levels  and 
Spot  Sizes  on  Aerodynamic  Effects 

Since  just  two  power  levels  (1  and  10  kw)  and  beam  sizes  (1  and 
10  cm2)  have  thus  far  been  tested  by  QUEST  in  the  simulation  of  laser 
melt-through,  only  limited  inferences  are  possible  as  to  the  influence 
of  each  of  these  parameters  on  the  aerodynamic  effects  on  laser-induced 
melting.  Notwithstanding,  the  similarity  in  form  and  magnitude  of  the 
corresponding  curves  within  the  upper  left  and  lower  right  graphs  of 
Figs.  10-13  (which  depict  the  aerodynamic  effects  for  1  kw/1  cm2  and 
lU  kw/1  cm2)  suggests  that  the  ratio  of  power  to  spot  size,  or  flux  inten 
sity,  may  serve  as  a  gauge  in  determining  the  effect  of  wind  on  melt- 
through.  It  is  clear,  not  only  from  comparisons  between  the  above  two 
and  the  remaining  two  graphs  within  each  of  Figs.  10-13  but  also  from 
logical  considerations,  that  flux  intensity  is  strongly,  and  possibly 
predominantly,  influential  in  determining  the  effect  of  wind  on  melt- 
through.  Low  intensities  allow  aerodynamic  cooling  or  heating  (which 
one  depends  upon  the  mach  number  and  the  melting  point  of  the  metal) 
to  strongly  manifest  their  respective  effects  of  melt-through  retar¬ 
dation  or  acceleration.  At  very  low  flux  intensities  convective  heat 
transfer  due  to  aerodynamic  cooling  (and  to  a  lesser  extent,  radiation*) 


*For  example,  for  flux  intensities  of  0.1  kw/cm2  on  0.2  cm  thick 
titanium  and  aluminum  at  mach  2.0  and  evaluated  for  conditions  at  the 
spot  center,  heat  losses  due  to  radiation  accounted  for  6.0  and  0.01% 
respectively  of  the  heat  absorbed  (the  higher  melting  point  of  titanium 
primarily  accounted  for  the  three  order  of  magnitude  difference).  Thus 
radiation  was  found  to  be  insignificant  except  for  metals  with  high 
melting  points  which  are  exposed  to  very  low  (<0.1  kw/cm2)flux  intensities. 
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could  even  counteract  laser  heating  prior  to  the  metal  rising  to  the 
solidus  temperature  such  that  melt-through  would  never  occur.  On  the 
other  hand,  increasing  the  flux  intensity  influences  the  melt  transfer 
effect  In  a  manner  quite  similar  to  that  of  reducing  the  metal's  con¬ 
ductivity.  That  is,  for  aerostatic  melt-through,  increased  flux  forces 
an  increased  axial  temperature  gradient  within  the  molten  metal.  Equa¬ 
tion  (9)  manipulated  to  the  form  of  Eq  (36),  mathematically  illustrates 
the  effect  on  the  axial  temperature  gradients  by  changes  in  either 
conductivity  or  flux  intensity.  That  is. 


T.  .  .  -  T. 


(36) 


or  equivalently. 


£  =  ££  (37) 

k  by  1  ' 

An  ultimate  consequence  of  large  temperature  gradients  within  a  suf¬ 
ficiently  thick  melt  layer  is  vaporization.  But  even  if  vaporization 
does  not  occur,  removal  by  aerodynamic  friction  of  potential  "heat¬ 
robbing"  melt  will  reduce  melt-through  times.  The  extent  of  the  reduc¬ 
tion  is  approximated  by  Eq  (34). 

Quantitative  evaluation  of  the  influence  of  flux  intensity  on 
aerodynamic  effects,  difficult  due  to  the  limited  data  points  and  the 
simultaneous  Influences  of  other  parameters,  was  attempted  by  plotting 
dimensionless  melt-through  times  for  0.2  cm  thick  titanium  and  aluminum 
at  sea  level  mach  2.0  against  logarithmically  scaled  flux  intensities 
as  shown  in  Fig,  16  (the  corresponding  times  for  stainless  steel  and 
magnesium,  although  not  plotted,  approximated  those  of  titanium 
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Fig.  16.  Melt-Through  Times  for  0.2  cm  Thick  Titanium  and  Aluminum  in 
a  Mach  2.0  Sea-Level  Airstream  Normalized  by  Dividing  by  the 
Corresponding  Mach  0.0  Melt-Through  Times  as  a  Function  of 
Power  Density 

and  aluminum  respectively).  The  solid  curves  as  deduced  from  each  of 
the  two  sets  of  data  differ  in  form— that  for  titanium  (and  stainless 
steel)  implies  a  decreasing  reduction  in  melt-through  time  with  each 
magnitude  Increase  in  flux  Intensity  while  the  one  for  aluminum  (and 
magnesium)  indicates  a  constant  reduction.  As  previously  stated,  how¬ 
ever,  with  aerodynamic  cooling  and  some  small  flux  intensity,  the 
melt-through  time  for  aluminum  must  go  to  infinity.  And,  no  matter  how 
intense  the  flux,  melt-through  will  require  some  positive  finite  time- 
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Thus,  the  curve  for  aluminum  In  Fig.  16  cannot  remain  linear,  but  must 
asymptotically  approach  both  the  abscissa  and  the  vertical  to  the 
maximum  flux  intensity  at  which  no  melt-through  occurs  for  the  given 
mach  number  and  the  other  initial  conditions.  Such  curvature  is  sug¬ 
gested  by  the  dashed  lines  of  Fig.  16,  but  computation  of  additional 
data  points  is  essential  prior  to  any  additional  postulations  in  quan¬ 
titatively  relating  flux  intensity  and  aerodynamic  effects. 

It  may  be  noted  that  the  two  data  points  for  each  metal  at  the 
1.0  kw/cm2  flux  intensity  in  Fig.  16  differed  slightly.  Further  inves¬ 
tigation  of  the  data  in  Appendix  D  will  reveal  that  for  similar  metals, 
thicknesses,  and  mach  numbers  the  1  kw/1  cm2  power-to-spot-size  ratios 
always  required  longer  to  melt-through  than  did  power-to-spot-size 
ratios  of  10  kw/1 0cm2— this  occurred  despite  identical  flux  intensities 
of  2000  joules  per  second  per  cm2  at  the  spot  center  (reference  Eq  (4)). 
The  differences,  up  to  50  percent  of  the  smaller  time,  are  ascribed  to 
lesser  radial  heat  transfer  rates  within  the  larger  spot  due  to  smaller 
radial  temperature  gradients  which  in  turn  is  a  consequence  of  the 
assumed  Gaussian  flux  distribution  described  by  Eq  (3).  Differentiation 
of  Eq  (3)  with  respect  to  the  radial  distance  r  provides  the  radial 
gradient  of  flux  intensity  as  a  function  of  the  spot  size  a;  radial 
distance;  and  flux  intensity  at  the  beam  axis  q o.  That  is, 

9r  a2  1  ' 

Easier  to  visualize  than  Eq  (38)  is  Fig,  17  which  shows  the  local  flux 
intensity  along  the  radii  of  1  and  10  cm2  spots  absorbing  1  and  10  kw  of 
power,  respectively.  The  intensities  and  their  gradients  are  identical 
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at  the  spot  centers,  but  the  magnitud-s  of  the  radial  gradient  can  other¬ 
wise  be  seen  to  be  less  for  the  larger  spot  at  any  point  within  the 
smaller  spot.  The  radial  temperature  gradients  within  the  metal  plate 
and  the  consequential  radial  heat  losses  will  approximately  correspond 
to  the  relative  radial  gradients  of  the  flux  intensities.  Thus,  larger 
radial  heat  losses  should  account  for  much  of  the  added  time  required  for 
melt-through  in  the  1  kw/1  cm2  total  power- to-spot-size  cases. 


Fig.  17.  Intensity  of  a  Gaussian  Flux  Distribution  as  a  Function  of 
Radial  Distance  from  the  Center  of  1  and  10  cm2  Spots 
Absorbing  1  and  10  kw  Power  (^0  -  2000  joules/(sec*cm2)) 

Extending  the  above  analysis  to  the  case  of  a  constant  total 
flux,  for  example  1  kw,  and  a  varying  spot  size  helps  to  explain  why 
spot  size  may  be  safely  ignored  in  gauging  the  accuracy  of  a  one¬ 
dimensional  heat  flow  approximation  for  laser  heating  (Ref  12).  Doub¬ 
ling  the  spot  area  from  1  to  2  cm2  exactly  halves  the  flux  intensity 
at  the  spot's  center  q<>  while  roughly  halving  the  radial  gradient  of 


6AW/MC/72-1 7 


the  flux, Intensity  (see  Fig.  18).  Halving  q0  will  about  double  the  time 
required  to  melt- through,  whereas  halving  the  radial  gradient  should 

^  v 

reduce  the  radial  heat  transfer  rate  by  a  factor  of  approximately  two. 
The  result  is  an  approximately  equal  radial  heat  loss  despite  the  doub¬ 
ling  of  the  spot  size.  Postulating  further,  varying  the  spot  size  for 
a  constant  laser  power  and  plate  thickness  would  have  negligible  influ¬ 
ence  in  the  total  heat  lost  through  radial  conduction;  or,  as  discovered 
by  Torvik  and  reported  by  him  in  Ref  12,  the  accuracy  in  approximating 
laser  heating  by  one-dimensional  (axial)  heat  flow  is  independent  of 
spot  size. 


Fig.  18.  Intensity  of  a  Gaussian  Flux  Distribution  as  a  Function  of 
Radial  Distance  Frun.  the  Center  of  1  and  2  cm2  Spots 
Absorbing  1  kw  of  Power 
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VI.  Cor.clusions  and  Recommendations 

Based  on  the  results  and  findings  to  date,  QUEST  appears  to  be 
a  liable  model  capable  of  closely  estimating  the  melt-through  times 
and  approximating  the  surface  recession  phenomena  of  laser-induced, 
win ’.-affected  melting  of  metal  plates  when  the  energy  absorption  rate 
is  known.  From  this  premise,  the  following  conclusions  and  recommen¬ 
dations  are  submitted. 

Conclusions 

1.  The  model,  as  comprised  in  the  computer  program  QUEST,  has 
mat  the  objective  of  this  thesis  study.  That  is,  it  stands  ready  as  a 
tool  for  further  investigation  of  laser  effects.  Substantial  data, 
ancillary  to  the  melt-through  times  as  were  primarily  reported  upon  in 
the  preceding  chapter,  are  provided  by  QUEST  to  further  assist  in  com¬ 
prehending  the  nature  of  melt-through.  The  computer  program  Is  adaptable 
to  a  wide  range  of  input  parameters  and  is  capable  of  accuracies  com¬ 
mensurate  with  the  computer  processing  time  available. 

2.  The  findings  thus  far  indicate  that  significant  utility  can 
be  obtained  from  QUEST  in  determining  aerodynamic  influences  cn  laser- 
induced  melting  that  would  otherwise  be  unpredictable. 

3.  The  effects  of  wind  on  laser-caused  melting  were  found  to  be 
by  two  mechanisms--melt  transfer  and  convective  heat  transfer.  Either 
can  predominate  at  a  given  altitude;  which  one  depends  primarily  upon 
the  flux  intensity  and  the  mach  number.  The  magnitude  of  the  aero¬ 
dynamic  effects  likewise  depends  upon  the  mach  number  and  flux  intensity, 
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but  also  varies  considerably  with  the  plate  thickness,  the  metal  and. 
the  size  of  the  laser  beam.  Further  use  of  QUEST  will  enable  a 
clearer  understanding  of  these  interrelationships  over  a  wider  range 
of  parameters. 

4.  Without  question,  the  computer  program  QUEST  can  be  improved 
.on.  A  major  restriction  (that  was  recognized  but  which  was  accepted 

since  resolution  thereof  was  beyond  the  scope  of  this  study)  is  the 
requirement  to  specify  the  flux  absorbed  by  the  metal  as  opposed  t.o 
the  flux  incident  upon  the  metal.  Incorporation  of  a  reflectivity 
model  within  QUEST  which  would  account  for  the  change  of  absorptivity 
as  a  function  of  the  metal,  its  surface  condition,  its  temperature,  and 
any  other  significant  parameters  would  enhance  the  practicality  of 
QUEST  in  the  overall  study  of  laser  effects. 

5.  Another  restriction,  similarly  accepted,  is  the  exclusion 

of  catastrophic  oxidation  effects  in  the  model.  Lesser  restrictions  are 
those  inherent  with  limited  available  data  on  the  properties  of  metals, 
particularly  at  le  elevated  temperatures  approaching,  and  in  excess  of, 
melting.  The  basic  logic  used  within  the  present  model,  however,  is 
unaffected  by  such  property  value  functions  and  therefore  will  accom¬ 
modate  revisions  as  they  may  become  available. 

Recommendations 

It  is  recommended  that: 

1.  QUEST  be  gainfully  employed  to  further  the  knowledge  of 
laser  effects  on  metals  iri  an  aerodynamic  environment. 

2.  Modeling  the  reflectivity  of  laser  light  on  melting  metals 
be  attempted  and  any  successful  results  be  incorporated  within  QUEST. 
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3.  Further  study  be  given  to  modifying  the  convective  heat 
transfer  rate  as  presently  embodied  in  QUEST  by  including  effects  of 
unheated  starting  lengths  and  variable  surface  temperatures  (reference 
Eq  (23)),  particularly  for  modeling  effects  from  low  flux  intensities 
or  with  laminar  air  flow. 

4.  The  effect  of  melt  transfer  due  to  gravity  be  compared  with 
the  effect  of  similar  melt  transfer  from  aerodynamic  friction  and  if 
the  effect  of  the  former  is  significant  that  it  be  modeled  within  QUEST. 

5.  Catastrophic  oxidation  as  a  predictable  event  within  laser- 
induced*  aerodynami c-inf luenced  melting  be  investigated  and  If  possible 
a  model  of  such  be  included  within  QUEST. 

6.  Continued  efforts  be  made  to  further  validate  QUEST  through 
experimental  results;  and,  despite  the  consistent  results  to  date  and 
the  scrutiny  to  which  the  logic  and  mathematics  of  QUEST  have  been 
repeatedly  subjected,  that  the  numerical  results  from  QUEST  continue  to 
be  treated  as  suspect  (as  should  be  the  case  with  any  computer-programmed 
model ) . 

7.  The  model  be  extended  to  other  metals  of  Interest. 

8.  Modeling  of  aerodynami cal ly  influenced,  laser-induced  melting 
of  metals  of  other  shapes  be  attempted. 
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APPENDIX  A 


QUEST,  a  Computer  Program  of  a  Numerical  Model  of  Surface  Recession 
PhenomenaTof  Metals  Subjected  to  LaserTadlaflon 
irT  an  Aerodynamic  Environment 
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PROGRAM 


5 


10 


15 


20 


25 


30 


35 


40 
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50 


55 


60 
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ouesin  coc  ssoa  ftm  V3.0-320A  opt»i  i>/i5/72 

PROGRAM  QUEST  Ti  (|(*1U?.0UTPUT.TAPE7)  - 
COff€NT***TlTANHjM«**TITANHjH*«*riTANIUM,*«TlTANlUM»«*TlTANIUM”*Tir«JiUK“‘ 
COrnON  T(30.35».TT(30.35).0Y(30,35).OCAXL(30.35).CfCR0L(30.35;. 

2  TFR(30.351.*A.'30).1«I30).NZ(30).NAT(30).N8T(30).TIX(35>  .T8X<35) . 

3  ATXrR(30>.8TXFR/30).TAU<30).TQCNVI301.7QRAO(30>.TVANlSH(30).ai3J> 

4  .O0UT130!  .QRAOI301  .SUMDYI301  .REQ13S1  ,  AC0EFF.AENrn.P.AKACH,  WELOX. 

5  APNSURE.ARHO.ARN.ATEH* .CPI  .CPSLUSH.CPV.0ELX.0ELY. DTIKT. E.MF.HV.  I. 

6  J.MIDINT.MINT.NINT.NB61 .NXTRA.N2I .PRANOTL.RHO.STMELT.TH. . TMU.TV, 

7  TISTMLT. TIME. VIS. XFR0MLE. STOP. 0ELTIME 

EQUIVALENCE  (DOT IME.XFMCL .YIN> . (XTlfC.ATXFR) . 1 TBACK.BTXFR.YMELT. 

•rvELOX) . (TFR.CAC1 .  IYSOC  IO.REQ) 

DIMENSION  ZVELQXI301 .TSACXI30) ,XFMCL(30l .XTIME130) .YINI301. 

:  7CLT(30)  .CNO130.35'  .DATAU0241  ,YSa  10(301 

0<  TA  TMLO. TKUP. TVP.HFS.HVP/ 1975. . 1968. . 3550. ,93. 9457?, 21 08. 5595/ 

DA  A  CMEMAX,CENSlTY/0.07.4,50/ 

•  CATA  >**071/0.7113/ 

CALL  PLOTS  (DATA.  1024.  71 
RE AO  903.  ALT.AMACH.XFROMLE.TO 
READ  993. Prf?.SIGMA2. RADIUS. THICK 
READ  992.KINT.NINT 

992  FORMAT (41201 

993  FORMAT  (4F20. 10) 

CALL  ATMOSIALT.ATEMP.ASlGMA.ARHO.ATHETA.AOELTA.CA.AfK.ACNO) 

ACOEFF =0 . 2* AMACH' 42 
APRSURE*A0ELTA*2I 16.22 
AVELOX«CA*AMA.CH 
ARN'ARHO*A\'ELOX«XFRO(fi.E/AMU 
AlRAr IO*ATEHP/500. 

AENiHLP* . OI836064*A7RaT 10* (0.9395396*ATRAT 10* (0 . 03243429-ATRAT 10* 

:  0.00099174211 
10  0ELX*RAD!US*2. /MINT 
DELY*THICK/N!NT 
E'0.0001 
Tf-HFS 
HV-HVP 

M10INT*(M1NT*I 1/2 

NINTT*NINT*I 

NINTTT*NINT»2 

NXTRA*35 

PRANOTL-PPPNOTL 

RHO'DENSITY 

STGP*0. 

SUMQCNVO.O 

SUMQRAO'O.O 

IML-fMLO 

TMU*TMUP 

TV«TVP 

CALL  SPH£aT((TML*TMU!/2.  .CPSLUSHI 
CALL  SPHCAT  (TML-2. ’E.CPI ) 

CALL  SPHCAT  (TV.CPV) 

CALL  SPHCAT  (T0.CP01 

DELT IME*D .5*RHO"CPO/CNDMAX 1 1 1 . / ( 1 . /DELX* *2* I .  /DEI  Y*  *21 1 

FLUX«PWR’230.9 

CALL  FLUXX( FLUX. SIGMA?) 

20  00  30  I  *  1  .MINT 
NA( t)*NINT 
NSI I l*NiNt 
N7 ( I 1«NINI 
QOUTl 1 1*0.0 
OR AO l 11*0  0 
TQCNVI 11*0.0 
TCRAOl 1 1*0  0 
T VANISH! I  1*0  0 
00  25  J»l .NiNT 
T<  1  . J) *  TO 

Tti ;  ,ji*to 

0Y( I ,J)*OCLT 
25  CONTINUE 

00  30  J-NINTt .nxTRA 
T ( ! ,J)*0  0 
TTI  I  ,J)*()  0 
0Y( I .  /'*0  0 
30  CONTINUE 
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program  oestti  page  5 

call  corcuc? 

75  Orif€*O.S*!»0»OeLY'OCLY*CPO/C.NOl  I .  I  > 

38  00  HO  1*1. MINT 

HO  r(|.NZ(I11*<0*Orif€2<«<)  *CP0*0Y<I.N7U)i)'0<ll 
iF(T!MIOfNI.NZ!NIOlNT>).LT  !TM-t>)  GO  TO  H5 
OnME*!TML-T0)/IT<MlOINT.N2IMIOtNI))-TO)*DTIME 
80  GO  TO  38 

H5  00  H5  I  *  I  .MINT 

IF! T ( I ,N2( 1 1 1 .GT . TH.-E)  NA!l)*NA(!l-l 
H5  T( 1 .MZC |l-ll*CTf| ,N2!  I l)-TO»/H.O*TO 

T  !.*«*<  TIHIOlNT.NZfMIOINT) -I  )-T01*CP0‘RHO*0eLY/O(MIDINT> 

85  GO  TO  60 

50  CONT  NUE 

00  53  I =1 .MINT 

TOCNVI I )  *QOUT  ( 1 1  *0T1F£*TQCNV!  I ) 

TQTAO! I 1*0RADI I )*0T1ME*TQRA01 I > 

90  SUCY < I i*0.0 

DO  v3  J*l .NXTRA 
53  SUM0>ll)*3UHDYII>*0Yt|.JI 
00  55  1*1 .MINT 
inN8UI.LG.Ci  GO  TO  100 
95  55  CONTINUE 

DTIME*D'LT1F£ 

CALL  FRtUON 
CALL  SCFM'DT 
CALL  XACPMG 
100  CALL  COOL  !T 

CALL  VAPOR  1 Z 

IF  (NB'MIOINTI .LG.01  GO  TO  60 
CALL  '.RNSFRX 
CALL  SQUEEZE 

105  60  00  70  1*1  .MINT 

65  00  70  J* I  .NXTRA 
70  TTl!.J)»TM  .J) 

GO  TO  50 
100  CONTINUE 

110  800  PRINT  1000.  THICK.  PWT.SICHA8.0ELX.DELY.ARN.  ALT.  AMACH.XFROM.G,  TO. 

:  TMU.  TV. HP  .HV.CPO.CPI  .CPSLUSH.CPV.CNDMAX 
1000  FCWMATHMI.'TITANIIW  *F|5.8.*  CENTIMETERS  THICK* / 

'■  1H  .‘TOTAL  BEAM  PC»ER**F9. 1 . •  KILOWATTS*/ 

:  IH  .  ‘BEAM  RADIUS  I?  SIWTAi  **FI0.5.  •  CM*/ 

115  :!H  .  *0ELXIRA0IAL  OlfCNSICN  (T  ELEMENTS) **F9. 6. •  CM*/ 

:IH  . *0ELY (  AXIAL  01 MENS I ON  OF  ELEMENTS » **F9. 6. •  CM*// 

.  IH  .‘REYNOLDS  NL«ER**Fi?.0,/ 

■  IH  .  *AL  T I  TUCE**F7. 0.  *  FEET*/ 

:  IH  ,*MACH  NUWER**F5. 3.  / 

180  :IH  ,*01STANCE  F0f*1  LEADING  E0GE*«F6.8.*  FEET*// 

:IH  .‘INITIAL  TEM°ERATURE*‘ r8.8. •  DEGREES  KELVIN*/ 

:1H  .‘FELTING  P0INT«*F|H.8,  *  DEGREES  KELVIN*/ 

IH  .’BOILING  POINT  ««F|3.8.*  DEGREES  KELVIN*// 

:IH  .‘FEAT  OF  FL7SI0N**F  I H .  3.  •  CAL0RIESIGM1/GRAM*/ 

185  IH  ,‘FEAT  OF  VAPORIZATION* *F8. 3. •  CAL0RIESIGM1/GRAM*// 

: IH  .‘SPECIFIC  HEAT  I  INITIAL  T«‘F7.6,‘  CALORIES/GRAM/OEG  K*/ 

•  IH  .‘SPECIFIC  l-EAT  I  :  T  MELT**F8.6.*  CALCftlES/GRAM/OEG  KELVIN*/ 

•IH  .‘SPECIFIC  FEAT  I  T  SLUSH**F9.6. ’  CALORIES/GRAM/OEG  K’/ 

IH  .-SPECIFIC  FEAT  i  T  UQUI0**F8.6.*  CALORIES/GRAM/OEG  K*// 

130  IH  .‘MAXIMUM  CONDUCT  I VI TY*  *F  18.6.  •  CAtORIES'SEC/CN/OCG  K*//l 

PRINT  996. 1  (T(  I  ,J1 .  1*1  .MINT) ,J<I .NXTRA) 

1998  FORMAT! IX, 87F5. 01 
PRINT  1999 

1999  FORMAT! ’IHINO  VECTOR  IS  TO  TIE  LEFT  IFOR  TEFPFRATURE  ARRAY  ABOVE 

135  l  AND  IS  UPWARD  (FOR  PRINTOUT  WHICH  FOLLOWS) */ 1 

PR'NT  5000 

8000  FORMAT!*  ELEFCNT *9X . ‘Q  IIN)*7X.*0  <CNVC>*8X,‘G  1RA0>'8X.*Q  (VPfif 
:6X.*CMS  VPRZD*7X.*CMS  LEFT*/) 

SUMQIN-0.0 

I  HO  SUMQCNV*0.0 

SUMQRAD'O  0 
SWMVPR*0.0 
00  300  1*1. MINT 
TO-Q! I )*TiME*OELX 

IH5  TQCNVI I )*TOCNV( I )*OELX 

TCRAOI I )»TORAOl I ) *CELX 
TQVPR»TVANISH<  I  l  *0ELX*RF!0*HV 

PRINT  <?C01  .  I  .  TQ,  TGCNV I  1 1  .  IQRAOi  I  i .  TQVPR,  TVANISHI  I  I  , SUWDY (  I  I 

8001  FORMAT! IG.7EI5  6/1 


ft 
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PROGRAM 

150 


155 


160 


165 


170 


175 


i  no 


195 


ISO 


IS5 


eoo 


?05 


210 


215 


220 


quest n  PACt 

SL*OIN*SUMQIN*IO 
SLMOCNV»SUMOCNV*TQCNVl  1 1 
SUMQRAD-SLMQRAO'TaRAO!  1 1 
SUHOVPR'SIMOVPR*  TQVPR 
300  CONTINUE 

HUNT  2660,  SUMO  IN 
PRINT  2620.SUM0CNV 
PRINT  2640.SUMQRAO 
PRINT  2O02.SUM0VPR 

2660  FORMAT l • | TOTAL  1NC10ENT  A8SOR0EO  E»CRGY**F25. 10*  CALORIES* » -  . 

2620  FORMAT!*  TOTAL  CALORIES  LOST  OUE  TO  C0NVECT10N«*F18. 10*  CALORIES*! 
2640  FORMAT!*  TOTAL  ENERGY  ENITTEO  OUE  TO  RAOIAT|0N»«FI8. 10*  CALORIES* l 

2002  FORMAT!*  CALORIES  EXPENDED  IN  VAPORIZATION  ( TOTAL t**FI5. 10/ I 
PRINT  2003. TITC 

2003  rORMATI*  TINE  TO  MELT  T(«0UGM**F!5. 10,*  SECONDS*"! 

DO  333  1*1. MINT 

YSCLIOf 1 1*0.0 

»7CLm * !  I -MIOINT  l  *OELX 
DO  332  J'l.NBBI 

332  YSai01|)*YSa!0ll)*0YII.J! 

IFiN8HI.EO.OI  YSOLIOl  1 1*0.0 

333  CONTINUE 
PRINT  20G4 

2004  FORMAT!*  OISTAMCE  FROM  C/L  ICMSl *5X. •  TOTAL  TH1CWCSS  <CMSl*5X. 

:  'NON-MOL  TEN  THICWCSS  ICMS)*! 

PRINT  2005. !  IXFMCL  1 1 1  .SLPOY !  1 1 .  YSOL  ID!  1 1 ! ,  1*1  .MINT! 

2005  FORMAT IF20.5.2F25. 12/ ! 

CALL  PLOT  10. 0.1.0. -31 

CALL  SCALEIXFMCL.5.00.MINT.  1 1 
CALL  SCALE (YSOL 10. 3. 00. MINT . ! ) 

SUNOY(MINT»l l^YSOL  IOIMINT ♦ I ) 

SLM0Y!M1NT*2I*YSCUD(H1NT*21 

CALL  AXIS  (0.0.0.0.36HRAOIAL  FtJSITION  FROM  BEAM  AXIS  (CHSI.-36.5  0 
•  .0  O.XFMCL  IMINI*  I  I  .XFMCL  (MINl*2! ! 

CALL  AXIS  (0.0. 0.0. I5HTHICKNESS  (CMS) . *I5.3.0.90.0.SUMDTIMINT*1I 
:  SLM)YIMINT*?!  I 

CALL  LINE!  XFMCL.  SUMOY,  MINT.  I.  *1.111 

CALL  LINE! XFMCL .YSOL 10. MINT. I.*!, II 

CALL  PLOT! IO.O.-I I .0.-31 

TBACKI 11*10 

T8ACKi2l*T!STTt.T 

XT IME « 1 1 -0.0 

XT  INC  !2I  *STICL  T 

DO  334  J*3,N1NTTT 

TPACKIJI*TIXININTTT-J*I I 

334  X!  IM  IJI*TBXININTTT-J*  1 1 
CALL  PLOT  (0.0, 1.0, -3) 

CALL  SCALE  IXTIME.5  00.NINTTT.il 
CALL  SCALE  IT8ACK. 7. 00.NINTTT.il 

CALL  AXISIO.O.O. .I9TCLAPSC0  TIME  (SECSl . - 19.5  0.0 . 0.XT IME ININTTT* I 
I  .  XT  ITC  !N!NTTT*2l  ! 

CALL  AXISIO ,0.0.0. 35HTEHPERATURE  OF  BACK  SURFACE  (OEG  K>.*35.7  0. 

90 . 0,  TBxCKININTTT  *11.  TBACKININTTT  t-?)  I 
CALL  LINE  (XT IME , TBACK .NINTT  T . 1 . *  I , I  I) 

CALL  PLOT ( 10. 0. - 1 .0. -31 
PRINT  3000 

3000  FORMAT!  •  I  TEMPERATURE  OF  8ACK  SURFACE  'DEG  Kl  VS  TIME  <5ECSi") 
PRINT  3005 

3005  FORMAT  I •  TIME  ( SECONDS* • I  OX . • TEMPERATURE  (DEGREES  KELVIN! •/ ) 
PRINT  300!  ,  ( I  XT  IME!  J)  .TBACKI  J)  I  ,J*1  .NINTH) 

3001  > CfiMAT I2C20  10/ I 
YMEL  T 1 1 1 *0  0 
XTHCl  ll-STMELT 
00  335  J*2. NINTT 
YMEL  T ! J) *OClY* ( J- 1 i 

335  XT|MEUI*T8X(NINTTT-JI 
CALL  PLOT  i0  0.1. 0.-3) 

CALL  SCALE  (Y.*ELT.5.00.NINTT.I  I 
CALL  SCALE  (XTIME.7  00.NINTT.ll 

call  AXISIO  0,0  0. 36H0 1  STANCE  FROM  ORIGINAL  SUW ACE  (CMSi ,  36.5  0. 


mi/nzm-M 


- 


-  '<■ 


m 

& 


335 


330 


335 


340 


345 


350 


355 


360 


365 


370 


OUtSf'l  PAGE  *• 

c  o.ttclknintt-:  i  .ymeltininit.?' i 

CALI  AXIS* 0.0. 0.0.  ISiHT  ITC  10  TCL'  ISECSl . *19.7. 00.90. 0.XI  IKiNiNI I 

■  * ! ! .xt ifC  inintt*?5 i 

CAU  Uf£  «rXLI.xri»t.N!NrT.|..*.|l) 

CALL  PLOTUO.O.-I.O.-3I 
call  PLOf  (0.0.I.0.-3I 
CALL  SCALE  '  XT1TC.5.00.N1NTT.  1 1 
CALL  SCALC(YTCLT.7.00.N1NIT.P 

CALL  AXISIO. 0.0.0.  (9HTITC  TO  *LT  (SECS! .- 19.5. 00. 00.0. XI  IMCININTI 
.  ♦  (  I  .XTITC(N1NIT*3: 1 

CALL  AX  15* 0.0. 0.0. ‘*901  STANCE  Of  INTERFACE  FROM  ORIGINAL  SURFACE  I 
iCHSI.^^.OO.SO.O.WCLTININITM  I .  Y>CL7(N1NTT*3*  I 
CALL  L INE (XTITC  .TTCLT  .NINTT  .  I  .♦  t  .  I II 
CALL  PLOT II0.0.-I.0.-3I 
PRINT  3003 

3003  FORMAT  I MSOL ID-L IQUIO  INTERFACE  POSITION  (CMS!  VS  TIME  ISECSl"! 

PRINT  3006 

3006  FORMAT!*  TITC  (SECONDS! • ICX. * CMS  FROM  ORIGINAL  5URFACE*! 

PRINT  3001 .  ( (XT|K(  Jl  .YTCLT  ( J)  1  .  J»  I  .NINTT ! 

ZvtLOX 1 1 !  «OEL Y/ 1 TBX IN INT 1  -STfCL  T ) 

XTlfCl  1  1*1  T0XININT)  *STMELTI/?.0 
FIN! I l*0£LY/3.0 
00  336  J*3.NINT 
YlN(Jl«J*OELY-0£lY/3.0 

ZVELOX(Jl«OCLY/lTBX<NINTT-JI-TBX(NtNTTT-JI  I 
336  XTITC(J1*(  TBXININTT-JPTBXININTTT-J!  1/3  0 
CALc  PLOT  (0.0. 1.0. -3) 

CALL  SCALE !XT  1TC.5. OO.NINT  .  1 1 
CALL  SCALE(ZVElOX.7.00.N1NT . 1 1 

CALL  AXISIO. O.O.O. I9TCLAPSE0  TlfC  ISECSl . - 19.5. 0.0. 0 .XT ITC ININT* 

■  P.XTITCiNINT«3I) 

CALL  AXISIO. 0.0.0.H4HVEL0CITY  OF  SOLIO-LIOUID  INTERFACE  (CMS/SECi. 

*4‘*.7.0,90.0.ZVELOX(NINT*1 1  .ZVELOXININW) ) 

CALL  L  INE  (XTIME.ZVELOX.NINT.I..I.IP 
CALL  PLOT (I0.0.-I.0.-3! 

CALL  PLOT  (0.0. 1.0, -31 
CALL  SCALE (YIN.5.00.NINT.  I ! 

CALL  SCALEIZVELOX.7.00.NINT. 1 1 

CALL  AXISIO. 0.0  0.36HDISTANCE  FROM  ORIGINAL  SURFACE  ICMS> . -36.5.0 . 

:0.0 . YINININT ♦ 1  I .Y!N(NINT*3II 

call  axisio. o.o.o.iwvelocity  of  sol id-l i quid  interface  icms/seo. 

■  »44.7.0.90.0.ZVELOX(NINT*I I ,ZvEL0X(NlNT*3l  l 
CALL  L1NEIYIN.ZVEL0X.NINT.I.M.IP 

CALL  PLOT (10.0. -1.0. -3) 

PRINT  3003 

3003  FORMAT (M SURFACE  RECESSION  RATE  (CMS/SEC)  VS  TIME  (SECSi  ANO  POSI 
:  T ION  OF  SOLIO-LIOUIO  INTERFACE*/! 

PRINT  3007 

3007  F0RMATI5X , *T ITt  I  SECONDS 1 '5X, *VELOC 1 TY ICMS/SEC 1  *41 . ‘CENT IMETERS* / > 

PRINT  3004 , ( (XT IME I J! .ZVELOXI Jl . YINI Jl ! .  J*l,NINTl 

3004  FORMAT (3C30  10/1 
CALL  PLOTE 

END 
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subroutine 


5 


to 


15 


20 


25 


30 


35 


*♦0 


■<5 


50 


55 


60 


65 


ATMOS 


coc  6600  tin  »3.o-2eoA  opr*t  11/15/72 


SUBROUTINE  ATMOSfZ.TM.S10KA.fiHO. TICTA.OEL TA.CA.AHU.K) 
CALLING  SEQUENCE 

CALL  ATHOS(Z.TM.e"GKA.f*O.TttTA.OELTA.CA.AMU.K» 

Z  *  GEOMETRIC  ALTITUDE  IFT) 

TM  *  MOLECULAR  SCALE  TEWCRATURE  I  DEGREES  RAMON) 

SIOMA  *  RATIO  OT  DENSITY  TO  THAT  AT  SEA  LEVEL 
PHD  -  DENSITY  Lt--SEC«'2-FT<M-4>  OR  SLUCS-FT**3 
THETA  *  RATIO  OF  lEMERATlfC  TO  THAT  AT  SEA  LEVEL 
DELTA  «  RATIO  OT  PRESSURE  TO  THAT  AT  SEA  LEVEL 
CA  *  SPEED  OF  SOUTC  (FT/SEC) 

AMU  *  VISCOSITY  COEFFICIENT  ILB-SEC-FT- '?> 


1  NORMAL. 

2  ALTITUDE  GREATER  THAN  300000. 

3  ALTITUDE  ICGAT1VE. 


FT. 


DIMENSION  HPRIMBf  II).  TfCf  1 1 1  .S!GMA8<  1 1  >  .ALMl  1 1 )  .ARAYI 1 1  .41 
EQUIVALENCE  (ARAYf I . I ) .TRRIlfif 1 1 ) . I AWAY 1 1 .2) . TM8t 1 1 ) . 

•  (ARAYf I .3) .SIOMASf 1)1. f ARAYf I .4) .ALMl I ) > 


DATA  f  f  ARAYf  I.J) 


36089.239 
82020.997 
154199.480 
173804.510 
259186.350 
295275.590 
344488. 190 
524934.380 
557742.780 
656167.980 


J- 1,41. 1*1. I 
518.666 
389.988 
389.988 
508.786 
SOB. 788 

298. 188 

293.188 

406.188 
2386.188 
■>566.188 
836. IS8 


I )/ 


l.OOOOOOOE  00 
2.9706958E-0I 
3.266575IE-02 
I .21 I7870E-0C 
5.867731  IE-04 
I .7229I56E-05 
I.7928595C-06 
9. 392I5I9E-08 
7.7658593E-I0 
5.6324877E-I0 
2.5726771C- 10 


-0.00356616  . 

0. 

0.0C164592  . 

0. 

-0.00246808  . 

0. 

0.00219456  . 
0.01097280  . 
0  00548640  . 
0.00274320  . 
0.00192024  I 


26 

17 

18 

9 

10 

1 1 
12 
!H 


OATA  Q  /  0.0 '8744 1 76  / 

X  S  /  135  72  / 

x  AMJZ  /  3.7.V998E-07  / 

X  TMZ  518.183 


IF(Z)26.I8.I7 

K*3 

GO  TO  13 

IF (Z.GT. 300000  )  K*K*l 
hPRIM-(R£/(RE«ZIW 
DO  10  t1*l .  1 1 

IFflF’RIM-H’filMBfMI  )  I  i  ,  12. 10 

CONTINJE 

M-12 

M»M-  I 

IF(ALMtM) 1 14. 15.14 

TM*TMe<M>  -ALrtfMI  •  fFRRIM-l-PPIMBlM) ) 


RE  /  2 . 085553 1 E  07  I 

PZ  /  2116.2  / 

RHOZ  r  0.0023769  / 


ATHOSOOt 
ATMOSO 02 
ATM0S003 
ATM0S004 
ATH0S005 
ATMOSO  06 
ATMOS  007 
ATHOS008 
ATM0S009 
ATMOSOIO 
ATMOSO 1 1 
ATMOSOI2 
ATMOSOI3 
ATM0S014 
ATMOSO I 5 
ATMOSO! 6 
ATMOSOI7 
ATMOS018 
ATMOSO I 9 
ATMOS  020 
ATMOS02I 
ATMOS022 
ATM0S023 
ATMOS024 
ATHOS025 
ATMOS026 
ATMOS 027 
ATMOS028 
ATMOS029 
ATMOS030 
ATM0SQ3I 
ATMO5032 
ATMOS033 
ATMOSO 34 
AYMOS035 
ATMOS036 
ATMOS037 
ATM0S038 
ATMOS039 
ATM0S040 
ATMOS04I 
ATHOS042 
ATMOS043 
AIM0SC44 

ATMOS045 
ATMOSO 46 

ATMOS047 

ATMOS048 

ATMOS049 

ATMOS050 

ATMOS05! 

ATMOS052 

ATMOS053 

ATMOS054 

ATMOS055 


SIGMA-EXPf  <  :  ,0*'0/ALMfMI  1  >  *  f  ALOGI TM3IM)  /  TM) )  1  •SIGMABfM)  ATMOS056 

00  TO  20  ATMOS057 

15  TH-TMBfH)  AiM0S056 

SIGMA-SI CMA8IM I  *EXPf- fQ*  (HFRIM-FPRIMefM! ) !  'TfC(H> )  AIMOS059 

20  RHO'RHOZ'SIGMA  ATHOS060 

TKTTA-TM/TMZ  ATMOS06I 

OELTA-SIGMA-TMTA  ATM03062 

CA-49  02177-SCRTf TMI  ATMOS063 

AKU»AMUZ*SQRT ( Th€ TA**3) 4 1  * IMZ*S> / f  TM«S> )  AIHOS064 

13  RETURN  ATMOS  065 

END  ATMOS  066 
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SUBROOIlfC  CONDUCT 


COC  6600  FIN  V3. 0*320*  OPT* I  11/15/78 


5 


10 


15 


80 


85 


30 


SUBROUTINE  CONDUCT 

COTt£NT*"TITANtUM«"TlTANIUM*"TITAN!UM"*T[TANIUh*"TlTANIUh**'mANlUM** 
CCmTN  Tl30.351.TT(30.35I.DYI30.351.C«MOa.(30.351.CNDRCIL  130.351. 
ti  TFR130.35l.NA(30>.fei30l  .NZI30)  .NAT  i3C  f  ,NBT  1 30 1 .  T1X1351 .  T8X1 351 . 

3  ATXFR130). BTXFR1301.TAUI 301. TQCNV130) .  TORADOOl .  TVANISHI 301. 0(301 

4  .O0UT1301  .QRAD130)  .sUOY(301  .REG135I  .ACOEFF.AENIH.P.AMACH.AVEI.OX. 

5  APRSURF.ARHO.ARN.ATETf’.CPI .CPSLUSH.CPV.OELX.DELY.Dl'iME.E.HF.HV. I . 

6  J. M  JOINT. MINT, NJNT.N88 1  .NXTRA.NZI  .PRANOTL  .RHO.STMELT .  T(t ,  TMU.  TV. 

•*  TISM.T. TIME. VIS. XFROMLE. STOP. DELT1ME 

DIMENSION  CNO130.35) 

EQUIVALENCE  ICND.TFR1 

DATA  CNO.CNI  .CN8.CN3.CN4/0 .075851 .  -0.00010874.  1 .3538160-07. 
:-6.894976E~!  1 . 1 .0830580-  IN/ 

NZl’NZdl 
DO  10  I =8. MINT 
10  NZi-MAXOlNZI.NZl'.n 
DO  80  1=1 .MINT 
CO  80  J*l .NZI 
m=nii.ji 

IFiTTT.GE.TMU-E)  GO  TO  15 

CNOl  l,.J)=CHQ»TTT*(CNI  *T r T  ■  (CN8* TTT *  1CN3* T  TT -CN4)  1 1 
GO  TO  80 

15  CNOl I.J1 '0.035 
80  CONTINUE 

mintt=mint-i 
00  40  1=1 .MINTT 
00  40  J=1 .NZI 

40  CNCRDLl I .Jl*8./( I ./CNOl I ,J>«! . /CNOl ! *!.JI l 
DO  30  1*1 .MINT 
NZII*NZ(1)-1 
DO  30  J-l  .NZI  I 

30  CNOAXL I  I  .  J1  =8.  /  1 1  .  /CNOl  1  .  J»  M  .  /CNOl  1  ,  JM  1  > 

RETURN 

END 


SU3ROUTINE  SPHEAT  COC  6600  ETN  V3  0-380A  OPT»l  11/15/78 

SUBROUTINE  SPHEAT  (TTT.CPIJi 

COMMENT '  #4TITAN!IW"<T1T  AN  IUM*'»T!t  AN  1UM#'#TITAN1UM***TIT  AN  1UM* "TITANIUM*** 
CCmiN  T130.35I ,  TT  1 30. 35). OY130. 351.  CNOAXU  30. 351. CNOROl  130.351 , 

8  TFR130.35)  ,NA(30>  ,T*3I  30!  ,N71 30 1  ,NAT  ( 301 ,NBT 1 30) , T IXI 35) , IBXI 35) . 

5  3  ATXFRl 301 ,BTxrR(30l . TAUI30I . TQCNVl 301 . TQRAD! 30 1 , TVANISH1 301 .0(301 

4  .QCUT 1 30 '  ,QRA(J(  30)  ,SUMDY(  301  .REQI35)  .  ACOEFF  .AENTK.P. AMACH.  AVELOX. 

5  APRScRC ,  ARt-IO.ARN.ATEMP.CPl .CPSLUSH.CPV.OELX .DElY .OT IME .E .HT ,HV . I , 

6  J.M10INT  .MINT  ,NINT  ,NB8I  .NXTRA.NZI  .PRANCTL  ,RHO .STMEL  T  .  TML  ,  TMU.  TV . 

7  TISTMLT. TIME. VIS, XFROMCE, STOP, DEI  TIME 


1° 

DATA  CO . C 1 . 08 . C 3 . C4 . C5 . CPI ' 7  1 94807E -08.8. S80977E -04.-5. 00948 1 E -7 . 
5.059365E-I0.-8.40G3R7E-I3.4  3878580-17.0, 167/ 
iriTTT  GE.TMU-E)  GO  TO  30 
in  TTT  OE  Tit -El  GO  TO  80 

10 

CPI  J>C0*TTT*(CI  •TTT<lC8<f'TMC3‘ITT*  (C4*TTT>C5)  I  n 

15 

GO  TO  100 

80 

TST=(TMU*TML)/8 

CPIJ'CO'TSTMCI  *TS!*(C?*TSTMC3*!ST*(C4*TST'C5i  Ill'MT  iTMU-Tft  1 

GO  TO  100 

30 

CPI J=CPL 

JO 

100 

CONTI  WE 

RETURN 

ENO 

71' 
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COOUT  COC  6600  FTN  V3.0-320A  OPT-i  11/15/7? 

SUBROUTltC  COOL  IT 

C<m*. Tr30.35».TT (30.35* .DYC30.35)  ,OOAXL(30.35I .CMKOLf 30.351 . 

3  TrR(30.35) .NA(30> .NBC30J .NZ<30) .NAT (30) .RBTI301 .TIXI351 ,TBXI35> . 

3  ATXFRC30)  .8TXTR(30)  .TAU(30)  .TQCNV(30)  .TQRA0I30)  .TVANtSH(30)  .0(30) 

4  .O0UK30)  .GRAOI301 ,SUK)Y(30)  .REQ(35>  .ACOEFF.AENTH.P.AMACH.AVELOX. 

5  APRSURE.ATWO.ARN.ATEW’.CPI .CPSLUSH.CPV.OELX.OtLY.DTIfC.E.W.HV.  I . 

6  J, MIOINT, MINT  .NlNT.tfiBl  .KXTRA.NZI  .PRANOTL.RHO.STMEUT.TM. . IMU.TV. 

7  TlSTM.T.TITC.VIS.XFflOMLE.STOP.OELTIHE 
DIMENSION  HfCT(30.35)  .0QT1IC(30) 

EQUIVALENCE  ( TFR.WCT  J .  (DOTUC.REQ) 

DO  10  1*1 .HINT 

IF  (OOT1T£(  I )  .NE.OTUC)  GO  TO  10 
NA( 1 )*NAT ( I ) 

*eil)*NBT(l) 

10  CONTINUE 

DO  100  1*1. MINT 
NZ1*NZ( I ) 

00  ICO  J-l.NZI 

50  IF(T([.J).G£.TT(1.J)>  GO  TO  100 

55  1F(((T(I ,J) .GE. TMU-E) .AM). (TT( I .J) . GE . TMU-E ) ) .OR. ( ( T( I .J1 .LE. TML- 

:  EI-AND.(TT(I.J).LE.TM.-E)).CR. ( ( ( T(  1 , J)  .GE. TM-E) .  AND. (T(I.J).LE. 

:  TMU-E) ) .AND. ( ITT( I ,J) .GE. TML-E) .AND. ( TT( I . J) .LE. TMU-E* ) ) *  GOTOIOO 
IF(f(l,J).LT. TML-E )  GO  to  65 
60  TTIME*DTIfC«(TMU-T(I.J))/(TT(!,J)-T!l ,JI) 

T ( I .  J)*TMU»TT1ME/RH0/CPSLUSH/DELX/0Y(  I  ,J)«»€T(  I  .J) 

IF (  NB(I).EQ.J-I)  NB(I»*J 
GO  TO  100 

65  ir(TT(l .JI.GT. TMU-E)  GO  TO  70 

TT|HE*OT|ME«(ftt,-T(|,J)!/<TT(|,J)-T(I.J>) 

T ( I ,J)«THL*TTIME/RHO/CPt/OELX/OY( I ,JI ‘LNETI I .J) 

IF  (NA(I>  .EQ.  J-l)  NA(i>J 
GO  TO  100 

70  TTIME«DTI(C*(TMU-T(I.J))/(TT(I.Jl-T(l.jn 

T ( I  .J)«TMU-TTlME/RUO/CPSLUSH/OaX/OY(  I  .J)  •LNETf  l  ,J) 

IF (  NB(I).EQ.J-I)  MB(I)*J 
IF( T( I , J) .GE. TML )  CO  TO  100 
TTIME*TTIME*l  TMU-T  ( 1 ,  J)  I  /  ( TMU-TML  ) 

T ( I .J)*TML‘TTIME/RHO/CPI/OELX'DY( I , J) *HNET ( I , J) 

IF  (NA( I )  .EO.  J-l )  NAI I >"J 
IF (  NB(D.EQ.J-I)  N8(  I ) *J 
100  CONTINUE 

TIME  -  TIME  ♦  OTIME 
STOP-STOP* I .0 

IF< (STOP. EG. 20- I  OR. (STOP. EQ. 30. I .OR. (STOP. EQ. 40. ) .OR. (STOP. FQ. 50. 

: )  .OR.  (STOP. EQ. 60. )  .OR.  (STOP. EQ. 70. )  .OR.  (STOP. EO. 80. )  1  F>0  TO  105 
IFiOTIME.EQ.OELTIME!  GO  TO  200 

105  PRINT  99S. TILE. OTIME. (NA(l).l-l, MINT). T (MIOINT ,NZ (MIOINT >) , 

:  T (MIOINT , I ! . (NB( I > , I  *  I .MINT ) .SUMDY (MIOINT ) ,  TAU(MIOINT) 

:  ,(NZ(I). 1*1. MINT) 

999  FORMAT ( 1 36H  ••»••**••••«'••* . . . 


. . .  T1ME**F 15. I2.3X, *DT IME**F 15. I2.4X, *NA( 11**3713. 

:/•  T (MIO.FRNT )**F8.3.3X, *T (MI0.8ACK) **F9. 3. 4X. *N8(II=,27I3./ 

:>  THICK*«FI4.  ie.3X,*0RAF,«*EI6. 10 ,4X.  *NZ(  1 1  **271 3) 

200  CONTINUE 
RETURN 
END 
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SUBROUTINE  FLUXX  COC  6600  FTN  V3.0-320A  OPT-I  Il/I5'72 

SUBROUTINE  FLUXX  1FLUX.S1QNA2) 

COttOK  TOO. 35)  ,TT( 30.35) .OY130.351  .CK3AXL  <30.351  .CN0RDL130.35) , 

•’  TFRS30.35)  .HA<30>  .NB130)  ,N7(30>  .NATT30)  .WK301  .TIXI35J  .TBXT351 . 

3  ATXFRI30)  .BTWRC30)  .TAU130I  .TQCW130J  .T0HAO13O)  .TVAMISHI30I  .QI30) 

5  4  .a0UT(30).GRAO(30>.SUK>YI301.flEQ(35).ACCCFF.AENT*.P.WCH  AVELOX. 

5  APRS*E.J»C.AflN.ATE»*\CPI.CPSLU94.CPV.OELX.OELY.OTINE.E  »C.h\\l. 

6  J.MIDiNT, MINT. HINT  .NEB I  .NXTRA.NZI  .PRANCTL  .TW5.5TfCLT.TH_ .  THU. TV, 

7  TlSTICT.TtfC.VIS.XFROHE.STOP.OELTIME 
00*2. ‘FLUX/13. !415926536-S1GMA2**2> 

10  MIOKIN‘MIDINT-1 

5  DO  10  W.MI0H1N 
R-DELXMMIOINT-I) 

01 1 t*Q0*EXPl-2.0* (R/S1GHA2)  “2) 

0(MINT*I-!)*0(H 
15  10  CONTINUE 

OfH101NT)*QO*EXPf-2.0‘<OELX/4.0/SIGMA2)‘-2) 

RETURN 

EK) 


SUBROUTINE  FRCTION 


COC  6600  FTN  V3.0-320A  0PT*l  11/15/72 
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SUBROUTINE  FRCTION 

COf*«N  T ( 30 . 35 )  . TT! 30.351  ,DY  130.35)  .CNOAXL  130.35)  .CfCROL  130.35) , 

2  TFR130.35) .NAI30) ,hBi30) .NZl30> .NAT130) .NBTI3CU  . T I XI35) . TSX ( 35) . 

3  ATXFR130) .BTXFR130) ,TAU(30) .TQCNVI30) .TQRA0I30) .TVANISHI30) .0130) 

4  .Q0UT1301  .QRADI30I  .S0WPI30)  .RE0I35)  .ACOEFF.AENTHLP.AMACH.AVELOX. 

5  APRSURE.ARHO.ARN.ATEMP.CPI.CPSLUSH.CPV.DELX.DELY.OTIIC.E.HF.MV.I . 

6  J.MIOINT.MINT.NINT.NBBI  .NXTRA.NZI  .PRANCTL  .RHO.STMELT.TH.TMU. TV. 

7  TISTMLT. TINE. VIS, XFROM.E. STOP. DELTIME 

OATA  E0.E!.E2.E3/i.336064E-02,9.395396E-01.3.24342fle-02.-9.9l742C 

2  -04/.E1NO.EINI .EIN2  EIN3/2.7l055aE-02. 1 .012903. -2. 01091  IE-02. 

3  5.602829E-04  /. AO. Al  ,A2.. '.3/0. 979006. 0.00651928. 0.0158536. 

4  -0.001 55636/ .B0.B1.B2.B3/0. 030633 1.1. 07642 . -0 . 1 1 3939 , 0 . 00679309/ . 

5  CP8ASE . V I SBASE . CND8ASE .  f«ASE / 0 . 2395 , 0 . 000255936.5 . dOt I B832E -05 

6  .66.377778/ 

00  100  1=1. HINT 
IFIAVELOX.GT.O. 1  GO  TO  5 
TAU1 I 1*0.0 
GO  TO  100 

5  TIH*T(  I  .NZM  i )  'C  .0036 

WH*E0»TIWMEI'TIMME2*TIM'E7) ) 

RFCTR-SORT(PRANOTL) 

10  RH»0.5*  (AENTHLP’NWI  '0 . 22‘RFCTR' ACC*’rF*AENTW.P 
RT«EIN0*RHM£INHRHME!N2*RH>EIN3) ) 

PRAW5TL-1 VI  SBASE*  SORT  <5. *RT/9. 1 *CPBA5E* (A0*RT* ! Al *RT* ( A?*RT* A3> ) ) / 
:CN08ASE/ <80*RT ■ (B1 »RT* (B2*RT *B3> ) ) ) 

RFCTR-SQRTIPRANDTL! 

IFIARN. GT. 400000. )  GO  TO  25 
15  RH*0.5* (AENTHLP*mi *0  22*RFCTR«AC0EFF*AENTHLP 
RT=EIN0*RHMEINl*RH*fEIN2‘RH«ElN3> ) 

RT£fP»RT*500. 

RMU* VI SBASE "SORT (RiEMP/900.0) 

ROENSTY-APRSURE/RTEMP/53.35 
RN'RDENSTY « AVELCX* XFROMLE/RMU 
CFX*0 . 664 /SORT !RNI 
GO  TO  50 

25  RFCTR-RFCTR** 0.6667 

RH“0.5‘  I  AENTH_P*-HH> ♦0.22,RFCTR,AC0CFF*AENTHLP 
RT*EIN0*RH‘(EIN1*RH*!EIN2*RH*EIN3) ) 

RTCMP-RT'500. 

RMU- V I S8ASC • SGR T I R TEMP / 900 . 0 1 
ROENSTY-APRSURE/RTEMP/53. 35 
RN-ROENSTY'AVELOX'XrROMLE/RMU 
CFX-0 . 370/  ( I ALOG101RN) )  *  '*2  504! 

50  TAUl  I  )*7.440937*Ar<H0*AVCt.OX**2*CFX 

RCVRYH-AENTH.P-hBASE'RFCTR’O  00001 1 1 0035* AVELOX • 

COMMENT  OOUT  IS  IN  CAI./SEC/CM*  *2 

OOUT 1  i  1  -0  244121  t76*CFX/PRAAOTL * *0  8667<ROENSTY*AVELOXMRCVRYH- 
WH'HBASF ) 

100  CONTINUE 
RE  TURN 
ENO 


rrrr 
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aBRourihc  freeze 
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SG8R0UTIJC; FREEZE-  .  .  .  . 

COMMON  TI30;35>  .TTI30. 35)  .DYI30.35)  .£>F,AXL(30. 35)  .CNDROU 30.35) . 

2  TFR(30.35>iNAI30).>B(30),NZi30)..NAT(30).WT<30)..TIX(35>.TBX(35). 
3;ATXFR(30WBWR(30>  .TAUI30)  ;TQCNV(30)  .TQfiAO<30>  .TVANiSH(30)  .0130) 

4  .0001  <30!  ,QRAO(30) , SUMOY 1 30 ) ;  REG  1 35 ) .  ACCCFF.',  AENTHLP . AMACH . AVELOX . 

5  AF>RSURE;ARf01ARN.'AT£f^;C?l'.CPSajSH;CPV,CieLX.DeLY.OTlMe.E.)r ,HV,  1 . 

6  J.MtB!NT.MlNT.NlNr;N8ei.NXIRA.N2I.;PRAN0TL.RH0.STMELT.TML.IMU.TV. 

7  TISTH.T. TIME, VIS. XFROH.E.SIOP.OELTIME 
IFIAVELOX.EG.O.O)  GO  TO  22 

NZI*NZ( I ) 

CALL  SPKAT I  (TML*T I ! .NZi-l  >  1/2,0. CPL ) 

CALL  SPFCAT ( I Tl I ,NZ1 ) ♦THUI/2. ,CP2) 

1FCT(I.N21-I).GE.TML-E)  GO  TO  10 

IFI  I TML-TI 1 .NZ1 -I ) ) *CPL*DYI l.NZI-l>.LT.II TMU-TML) «CPSLUSH*(T( 1 ,NZI 
:  1-THUJ  *CP2) *DY(  1  .NZI )  I  GO  TO  8 

T(l,NZI-l)»Ttl.NZI-ll*(ITN.-Ttl  .NZi-l) )*CPL*(TMU-TML)*CPSLUSH»IT( I 
:  .NZ1 l-TMU! *CP2)*DY( I .NZI )/ I IOY( 1 .NZI-1 )*DY( I ,NZI ) >*CPL) 

GO  TO  20 

8  T(t,NZl-l)"TML*((<TMU-TML)«CPSLUSH*!TII.NZl>-TMU)*CPe'*DY(I.NZI)-( 
:TML-T(I.NZ1-1))*CPL*DY(I.NZI-I))/(<0Y(|.NZI>*DY(1.NZI-I))»CPSLUSH> 
NAi I >*NA( I l-l 

1FITU.NZI-  li.GT.TMU)  CO  TO  9 
GO  TO  20 

9  T I 1  .NZI-I )»TMU«  (Ttl  .NZI-l l-TMU) *CPSLUSH/CP2 
IFII.EQ.MIOINT)  T1XINBI 1>  >ET( 1,1) 

IF ( 1 .EQ.MIOINT)  TBXINBI 1 ) l^TIME+DTIME 
N8( I )*N8( 1 l-l 
GO  TO  20 

10  1F( (TMU-TI 1. NZI-I ))*CPSLUSH«DYfl. NZI-I ).LT. (Til. NZI l-TMU) *CP2* 

:  DYI1.NZD)  GO  TO  12 

T  I  1  .NZI-I  >*T( I ,NZI - 1 ) ♦  1 1 (TMU-T 1 1  .NZI-I ) )  *CPSLUSH*( T  ( I  ,NZ1 )  -THU)  * 

:  CP2> )«OY( 1 ,NZ1 I / IIOYI I .NZI-l l*DY( I .NZI ) )*CPSLUSH) 

GO  TO  20 

12  T(l,NZl-l>*rMU*((T(t,NZI>-TMU>*CP2*DY(T.NZI>-lTMU-T(I.NZl-l>>* 

:  CPSLUSH«OY( I ,NZ1 -I ) ) / (OYI I ,NZI )*DY ( I .NZI-I ) )/CP2 
IFI I. EQ.MIOINT)  TIXIN6I I ) )*T( I , I ) 

IFII .EQ.MIOINT)  T8XIN8I I ) )«TIME«DT1ME 
N8I I )  *NBI  1 ) - 1 

20  DY1 I .NZI-I )*0Y< I .NZI-I >*0Yi I .NZI ) 

OYI I ,NZ1)*0.0 

Tl I .NZI )*0.0 
NZ 1 1 1 «NZ 1 1 ) - 1 

21  IF (OYI I ,NZI 1 ) ) .LE. < I .5*0ELY) )  GO  TO  22 
OYI I ,NZ 1 11*1 )*0Y( I .NZI I ) l-DELY 

Til ,NZ( 1 1 ♦ I ) *T 1 1 .NZI 1 1 ) 

OYI I ,NZ 1 1 ) ) *0£LY 

IflTd.NZIin.LT.TML)  NAI  !  )  *NAI  I ) » l 
IFI Tl I .NZI 1)1 .LT.TMU)  NBI I )«N8( 11*1 
NZ  1 1 )  *NZ  1 1 )  ♦  I 
GO  TO  21 

22  CONTINUE 
RETURN 
END 
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10 


15 


20 


25 


30 


35 


40 


45 


suaROuTiic  somor 

COmDN  T(30.35).TT(30i35).DY(30.35>  OCtAXU30.35)  .CMW0LC30.35) . 

2  TFRI30.35) ;NA(30>  .W(30l  ,NZ(30)  .NAV.(30>  .NBTI30)  JtX(35)  ,TBI((35> . 

3  ATXFR(30).8TXrR(3Bl.TAU(301.TCCMVf301.TGRAO?30).fVANlSH(301.Q(301 

4  ,OOUT(3Q)  .0RAOI30)  ,SUCY(30)  .«Ci!35) , ACQEFF-;  AENTHLP .AMACH . AVELOX . 

5  APRajRE.AI^.AW.Amff  ,CR»  .CWLUSHiCPV.mx.CCLY.DTlME.E.lCiHV.  1 . 

6  J.MID1NT,M1NT^N1NT11®RI  .NXTRA^.MZI  iPRANOTLi  WO.STTCLT.TM.  .TMj.TV. 

7  TiST«.T.Tif€;vis,xFdcn.E.STCff.0CLri«: 

OIHENSION  l#CT(30.35>' 

EQUIVALENCE  (TFR.WCT) 

OAT A  STEFAN/ 1.35506 IE- 12/ 

CALL  COCUCT 
00  100  1*1 .MINT 
NZI»NZ( I > 

00  100  J*I.NE1 

CALL  SPlCATITTl I .JI.CPIJ) 

IF(l.fC.l)  00  TO  5 

HXFRL*0.0 

CO  TO  10 

5  HXFRL*CNCROL(l-l,..il/OELX4AMINl(OY(I.J).DY(l-!.J>>MTT(l-l.J)-TT(l. 
:  J)> 

10  1F(I .NE.MINT)  GO  TO  15 
HXFRR«0.Q 
GO  TO  20 

15  HXFRR*CNOROL( I .  J) /0ELX*AHIN1 (0Y< I .J) .OY( 1*1 . J) I 4 CTTI 1*1 . J) -TT( 1 . J) 
:  ) 

20  IFIJ.NE.  DCO  TO  25 
HXrR6*0.0 
GO  TO  30 

25  HXFR8«CM)AXL(  1 . J-l  1/0CLV«ICLX« (TT(  1 .  J-l  1-TTI I . J1 1 
30  IF (J.NE.NZ1 1  GO  TO  35 

IF(T( I .J) .GE. 950.0)  GO  TO  32 

ORAO( 1 1*TT( 1 . J1 4 ‘H'STEF AN* (0. 08*0.000 I 018101 92* TT ( 1 .  Jl ) 

GO  TO  33 

32  QRAOI 1 >*TT< 1 . J) 444*STEFAN4 (0.20*0.00001 17647*TT( I . J) ) 

33  HXrRT * (Q(  1 1  *C0UT ( 1 1  -GRADt  1 )  1 40ELX 
GO  TO  40 

35  HXFRT-CNOAXL ( I .  JI/0ELY*0ELX4 ( TT(  I  .J*l  )-TT(  1  .J) ) 

cohk:nt*4  4  4To  make  model  one:  dimensional,  eliminate  hxfrr 

40  TNETI I ,J)*HXFRL*HXFRR*HXFRB*HXFRT 

T(  I , J) *TT ( I  .JHOTlfC/RHO/CPlJ'DY!  I  .J>/0£LX4WET(  1 .  J> 

1F( J.EQ.NAl  1 1 1  ATXffil  I  )»RH04CPS  J4DY(‘I  .  J)  *OELX/H«:T(  I . J) 

1F(  J.E0.N8I  I >  >  8TXFRI 1 1  *RH04CP1 J40YI  t .  Jl  ‘DELX/WCT  ( 1 .  J) 

100  CONTINUE 
RETURN 
END 


HXFRL  IN  NEXT  CARO 


SUBROUTINE  V1SCSTY  COC  6600  FTN  V3.0-320A  OP1-1  11/15/72 

SUeMOUTIlC  VISCSTY  (TXI 

CGTtCNT44  *T!  TANIUN444TI TAN!lH',44l  I  TAH1UM444T1  TANILM*44TI TANU31444?  1 TANIUM* 44 
COTMON  T(30. 351. TTI3B. 3*1. 07(30, 35). CNOAXLI30. 35)  .000.(30. 38). 

2  TFRI30. 35)  .NA(30>.»«<30), *7(301. MAT(30).NBT<30>,T1X(35I.T*((36). 

5  3  ATXFRI 30)  .8TXFRU0) . TAU(3Z  1 .  TQCW(  301 ,  TORAOI 30) ,  T VAN  1 SW 30)  ,Q(  30) 

4  ,Q0UT(30)  .QUAD!  JO)  .SUMDYI JJO)  .REQ135) , ACOEFF . AENTHLP.  AMACH, AVELOX. 

5  APRaJRErtUWD.ABN.ATOff  ,CP  i.,<mU6M.CPV,0ClX,0EtY.DTIME.E.FF  .W.  1 . 

6  J.MIOINT. MINT. MINT .Wir.VKTRA.Wri  .PRANDTl .RHO.STMELT . TML . TMU.TV. 

7  T1STH.T, TIME. VIS.WWW.E. STOP. 0CLTI1C 

10  DATA  BIG/ 1.E300/ 

1FITX.GE  THJ-El  00  TO  10 

V1S*0IO 

GO  TO  100 

10  VIS*0.064EXP( ( 1915. O-TX); 125  0>  *0.005 
15  100  CONTINUE 

RETURN 
ETC 
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SUBROUT 1NE:  SQUEEZE 

COtCN  1(30.35) .TTI30.351 ,0V (30. 351 .CNDAXLI30.3SI .CNDROLI30.35) . 

2  TFRI30.35).NA(3Q).fei30>.NZ(30>.NAT(30).NBT(30).TIX(35).TBX(35>. 

3  ATXFR(301.UTXFR(30).TAU(30) .TdCNV(30) ,TQRAD(30).TVAN1SH(30).Q(30) 

4  ,QOUT(301'.GRAO(3'M  .SLH)Y(30)  .REQI35) . ACOCF1? , AENTH.P . AMACH, AVELOX . 

5  APRSURE.ARHOi  AflN.ATEl^.CPl .CPSLUSHICPV.DELX.OELY.OTIME.E.HF.HV. I . 

6  J, MIDI NT iMI N>. MINT <NB81  ;NXTRA,NZI  ;PRANOTl,,PHO.STTELT,Ttt_.TNU.TV, 

7  T 1SIM. T .  T IME . V I S ,-XFROMLE .STOP. DELT 1MC 
IFIAVELOX.EQ.O.O)  GO  TO  300 

DO  200  1=1 .MINT 
NB0I *NB( I )*1 
NZI-NZIl) 

1F(NZ1 .LT.NBS! !  GO  TO  200 
N8PLUS2*NB( 1 ) *2 
NZIl«NZd)-l 
SUNDYLQ'O. 

DO  1  J*NBU1 ,NZI 

1  S!  HOYLQ^SUMDYLQ  >DY  I 1 ,  J) 

2  IFIOYI 1 ,N8( 1 1 )*SUMOYLQ.GT. 1 .5*DELY>  GO  TO  5 

3  Cl  4-  J-NBPLUS2.NZI 

T 1 1  ,NBS!  I* (OYI I  .NBBI  )  *  T  ( l  .N06I )  *DY(  I ,  Jl  4T( 1 , Jl) /(OYI I  ,N08I )  *OY(  1  ,  J 
:)) 

DYl I ,NB6I I»DYF1 .NBBI ) *OY( I . J) 

T  ( I  ,J)«0.0 

4  OYII.JI-O.O 
NZIll’NBB! 

NZ 1 *NE8 1 
CALL  FREEZE 

IF  (OYI 1  ,MZ(  I ) )  .GT.  1 .5*D£LY)  GO  TO  140 
GO  TO  200 

5  IF(0Y( 1 ,NB( 1 ) I-0ELYI3. 15.10 

10  T ( I ,NBPLUS21»(Tl I ,NB6I ) *OY( I ,NB8i )»T( 1 .NBPLUS2) *DY( I .NBPLUS2) )/ 

: (OYI I .NB8I 1 *OY( 1 .NBPLUS2) ! 

OYI  1  ,N0PLUS21*OYI  I  .i«81 )  *OY(  1 ,1«PLU52) 

Til .NBBI )*T( 1 ,N8t I ) ) 

OYI I .NBBI ) «OY( 1 ,NB( l ) ) -OELY 
OYI 1 ,NG( I ) )*OELY 
ie(l)-N0(|)+l 
NBBI  -N01 11*1 
NBPLUS2“NBPLUS2»I 

1FITI I ,NB( 1 ) ) .LT.THL-E)  NAI 1 l*NAt 1 )♦! 

GO  TO  3 

15  00  70  J*NB61 .NZ1 I 
REQU)*DELY-DY(  I  ,J) 

IF IREQI Jl )  20.70.40 
20  IFlJ.CO.NZIin  GO  TO  70 

Tl 1 ,J»1 l.TMJt(O) 1 1 , J»1 ) • I T I  I ,J»I l-TMUI-REOIJMITI 1 .JI-TMU) )/ (DYI I , 
:  JH  1-FVEOIJI  1 
OYIl ,J*I )*DY( 1 ,J*I l-REQ'J) 

OY!|.JI*OELY 
GO  TO  70 
40  JJ*J*1 

00  60  K'JJ.NZI 

IFlOYI 1 .K) .EQ.O. )  GO  10  60 

IFIDYll.XI.LT. REOIJI)  GO  TO  50 


GAW/MC/72- 


suhroutijc 


60 


65 


70 


75 


80 


85 


90 


95 


(00 


5*5«  -/: 


>Km&*tnmurm.  imatMinm  mtfvmes a+urn. 


S0LEC2E  MCC  8 

OYii.to-OYU.icmeafji 

I ( I .  JI*TNL'*(0YU ,  Jl*i  Tt  t .  J)-T(U)*fCQ(  J)  *  (T(  1  .Kl-TMUI 1/QELY 
OY(l.J»*OaY  ■ 

GO  TO  70 

50  T 1 1 , JI*TNU*(DY{  |  tJl *(T(  t , J1-TNU>*0Y!  I  ,K)  *  I  T(  I  .Kl-THUl 1/  I0YC  I ,  Jl» 

:  0YII.K1) 

OYCI.a»»DYCl.J»«OY(I.K» 

DYI  l.tO*0:0 
REOtJl*0ELY-0YM.Jl 
60  CONTINUE 
70  CONTINUE  , 

DO  80  J«NBB1,NZI 

IFYOYI  t  .J1.GT.0. 1  GO  TO  80 

NZil)*J-l 

IN*J 

GO  TO  85 
80  CONTINUE 
GO  TO  100  ‘ 

85  00  90  J*’N.NZI 
90  T ( I ,J)*0.0 

100  IF (0Y( I  ,NZ( 1 11  .LT.O.5*0CLYI  GO  TQ  120 
IF(0YII.NZI|11.GT.1.3*GCLYI  CO  TO  140 
GO  TO  200 

120  IF(TU.NZU1-II.GC.TMU-E>  GO  TO  125 
CALL  FREEZE 
GO  TO  200 

125  TU  ,NZI  11-1  l*THU*(0£LYM  Til. NZt  11-11 -TMU)«0Y(1.NZM>>  «(TI!.NZ!m- 
:  TMUI 1 / IOELY*OY(l ,HZ( 1111 
0Y< I ,NZ( 1 1-1 ••OCLY^OYI I .NZ(  1 1 1 
OYU.NZI  111*0.0 
T <  I  ,NZU  1 1«0.0 
NZ 1 1 1 *NZ (11-1 
GO  TO  200 

mo  NEXCESS«OY(I.NZ(I»WOELY 

IFIOYI I ,NZ( I l !-NEXC£$S*OELY.L£.0.5*t£LYl  NEXCESS»NEXCESS-I 
NZX*N7l 1 1 +NEXCE5S 
NZ  I *NZ ( 1 1 ♦ ! 

OO  130  K-NZI.NZX 
DYU.K>«0ELY 
130  TII.KI-TU.NZUll 

OY ( I .NZX)«OY( I ,NZ( 1 1 1-0ELY‘NCXCESS 
OYU.NZI  1 11-0CLY 
NZ(1I*NZX 
200  CONTINUE 
300  CONTINUE 
RETURN 
EM3 
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SUBROUTl*  TRWSFRX 


COC  6600  FTN  V3.0-320A  OPT*!  U/15'72  » 


5 


10 


15 


20 


25 


30 


35 


40 


45 


50 


55 


60 


65 


70 


SUJOJTUC  TRN9FRX 

COWCN  T130.35S  .TTI3G.35) . 0Y(30. 35)  .OCAJO.  <3C, 351  .OCROL  ( 30.35* . 

2  TfR130.35t.f*A<30t.W13Qr.NZI30:.HAT130>.HBT133l  ;TtXI35) .TBX5351 . 

3  ATWR130)  ,8TXFRI30>  .TAU(301.T0CW1301  ,TCRA01301  .TVANISH1301  ,Q(30> 

•**  .O0UT<30>JJRAO1301.9UK)ri30»  JCai3S>,AC0EFF.A£NTH.P.A*<ACK.AVELCX. 
5  «=RSlflE./«W).AflN.ArCt<P.CPl  .CPSLUSH.CPV.OeLX.GEUY.OT  jft  .  I , 

6-D.Hl0lNTiMINT,^MINT;KaB?iNXrRA;MZ!:.PRAK)TL.fMJ.STItLT.TML.THJ.TV. 

7  nCTHL'riTifCiVIS.XOOtEvSTOP'.OeLTlft  ' 

IflAVELOX.EQ.O.OI  GO ''TO- 300 
5  CONTINUE' 

DO  10  1*1 .MINT 
00:10  J“l  .NXTRA 
10  TFR(!.J)*0.0 
DO  30  1*1 .MINT 
NB8l*fe<ll*I 
NZi*NZU> 

IFlNZI.LT.f«3l>  00  Tc  30 

V1SY«0.0 

00  25  J*NBBI ,N2I 

CALL  V1SCSTY  tTU.Jl! 

15  'V1SY*VISY*DY1 I . JI/VIS 

16  TFRX*TAUt 1 1  >0T1ME*1VISY-0.5*DY1  I  .JI/V1SI 
If f IFRX.LT.OELX)  GO  TO  20 
KSHIFT*TFRXi'QCi.X 
SHIFT2*ITFRX-KSH|FT‘D£LX>/DEI.X 
SHiril*1.0-SHirT2 

lFlKSHlFT.EQ.il  GO  TO  18 

OY 1 1 -KSHIFT*!  ,1121 I -KSHIFT* I 1*1)  “DY1 1 . J)  *SH1FT1 

TFR1 I -KSHIFT*!  ,NZ(  1-KSHIFT*1  1+1 )*DY( '  .J) ‘SHIFT  1 

T 1 1 -KSHIFT* I ,NZ ( I -KSHIFT < 1)»1)*T(1.J1 

0Y(  S -KSHIFT. NZl  1-KSHIFT1*!  )*DY(  I  .J)*SHIFT2 

TFR1 | -KSHIFT ,NZ< I -KSHIFT 1*1) *DY( ! ,  J)  ‘5H1FT2 

Tl 1-KSHIFT ,NZt 1 -K5K1FTI*! )*Tt I , Jl 

DY1|.J)*0.0 

Td.Jl-0,0 

NZ  C 1 1 *NZ 1 1 ) - 1 

NZl I -K5HIFTI *NZl I -KSHIFT) ♦! 

NZl I -KSHIFT* 1 1*NZ( 1 -KSHIFT*! )*l 
GO  TO  25 

IB  0Y(l-l,NZH-|l*n*SHIFT2*OY(l.JI 
TFR1I-I  .NZl 1-11*11  *SHIFT2*OY ( 1 , J) 

Tl  1  -1  ,N71 1-1 >*l  )*TU  ,J> 

TFRII.JI ‘SHIFT l«OYII,J> 

OYll.JfSHIFTMOYII.JI 
NZ1|-II*NZ1I-U*I 
00  TO  25 

20  TFR1 1 .JI^TFRX/DCLX'DYl I ,J) 

25  CONTINUE 
30  C0NVINUC 
H!WT*M1NT-| 

00  60  1*1 .H1NNT 
NILQ0*NB(|)*1 

N2LCO«N81 !  1  *NZ! 1*1)  -M3I 1*1) 

IFl  (NZI I )  -LT.Nll.QOl  .AND.  1N2lQQ.EQ.NB1  111)  GO  TO  60 
IFlNZin.LT.N2l.QDI  NZUI-N2L00 

NZI-WI 1 1 

DO  50  J*NlLCn.NZ! 

JJ«J-N81 1 1 *NBI 1*11 

0EMDM*0YI  1  ,  Jl  -  1FRI I  ,  J1*TFR(  I  ♦  I  .  JJI 

IFlOENOM.NE.O.Q)  GO  TO  40 

T(I.JI*0.0 

Nzm*Nzin-: 

00  TO  50 

40  T 1 1  .  J!*l  (OYI I  .JI-TFRt  1 .  JM  •  Tl  1 . J! ♦TFRI  1*1 ,  JJ1  *T  11  ♦  I ,  JJ1  l/OENOM 
50  0Y(I.J)*0EN0M 
60  TONTINLC 

*HLCDVNB!H1NT )  ♦  I 
NZI«NZ1MINT) 

DO  70  J-NILCO.NZI 

70  OY1MINT  ,  Jl  *OY  1MINT.J1-TFRIHINT.J) 

300  CONTINUE 
RETURN 
END 


85 
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«3j8R0UTIlC  VAPORIZ  COC  6600  rtH  V3  O-EOA  OPT-I  ll/IS/V 

SU2W OUTIME  VWRi; 

COtON  T 1 30  . »1. TV(30. *)  .0Y<30. 351. OOMl( 30. 361. 000.130.351. 
a  TFR(3O.35l.Mt(KI.I0(3O).N7i3C>.MTI3O!.»eT(3Ol.TIX>35).rBXi35l. 

3  ATXTRI  30)  ,BFXF>t(30l  .TAU<30) .  UXNVi  901 ,  TORAOi 30> . TVANISHI 30)  .Ql  30i 
5  •<  .0OJTI30l.aRAO13O*.9lK>Yl3C>.flEai35I.AC0EFr.ArNTH.P.AMACH.AvC1.0X. 

5  V>RSi«r.A»40.A»<,ArCtf».CP|  .CFSLU9T.CPV.0EI.X.0ELY.0T(TC.E,tC.hv.I. 

6  J  ,M|OINT.K|NT.NINI.!»l  .KXTOa.WI  .PRAACTl  .RMO.STTCL  T  .  Tft. .  TMJ.  TV. 

?  I1ST«.T.TII«. VIS. XFWK.C. STOP. OELTtfC 

OO  100  1*1 .MINT 

10  5  IFiTii.ltflti)  LE  TV)  GO  TO  100 

DANISH-  I Ti | ,N?I I > 1  -  TV)  "CPV'DVI  I  ,M?l  I  1 l/HV 
iF<VANl&)-OYI  I  .NZl  I  1 1  I  iO. 20.30 
10  OYI  I.NZ(tl)«0YU.NZ<l>l-VANl6H 
T  I.NZHH-TV 

16  TVANISM!  !  fTVANISMI  I  1‘VANISH 

00  TO  100 

?0  OYt  |  , NZl I > ) *0  0 
T ( I .Hil I  I  1 »0.0 
NZl  1UN71 1 1  -  I 

;0  TVANISHI  I  I  «TVAN!9U  I  )  ♦VANISH 

CO  TO  100 

30  VAN1<X*VANISH-DYII.NZ<I>> 

TVANISHI  I  !*TVANISHTI  fOYl  I  ,NZ(  I  1  I 
OYll.NZ(l))*O.Q 

as  Tu.NZdit-o  o 

NZl I l-NZl I  I  - 1 

Ti  ,  .NZl  I  1 1  »T «  I  ,NZl  1 1)  ♦VANISH*HV/CPV/OEl.Y 

GC  to  5 
too  CONTINUE 
30  RETURN 


subroutine  xactii*;  cy  geoo  ttn  v3  o-s?oa  opdi  ims/7a 

SUBROUTItC  X.'XTIK 

COCWN  Tl30.35).TT(30,35l.0YI30.35l.Oe  JCL)30  35>  .CNOROL 1 30 . 35>  .■ 

S  TER( 30 . 351 .MAI 301 .NBI30I .NZl 30 > .NAT! 30) .NBT 1 30) . TIXI 35) . T8XI 35> , 

3  ATxrfl»30).BTWR(»).TAUI30l,TaCI*V(301.TQRADI30l  TVANISHI  .TO)  .01 30) 

5  *♦  .CJXIT (  30)  .OPAOI 301  .SUTCY 1 30)  ,REQl  35< .  ACOEFF  .ACNTH.P,  AHACM.  ‘VEI.OX. 

5  APRSUHE.AWHO.AflN.  A  TEW.  CP  I  .CPSLUSH,QPV,OCLX.tJES_Y  .OT  INf  ,E ,Hr  ,HV .  I . 

6  J.KIOINT. MINT  .MINT. NBBI  .  NX  TRA ,  NZ I  .  PRANDTl  .HHO.STMEl  T  .  TNL ,  THU.  TV . 

7  TISTMLT.TIPC. VIS. XFRCW.E . STOP. IXLTINE 
OlfCNSION  DOT INE 1 30 ) 

10  EQUIVALENCE  (DOTINF.dfqi 

5  DO  50  J - l .MINT 
COT  DC l  I  1-OTItOE 
OTIHA-OTIftfE 
OTIW-OTlfCH; 

15  NATI|)>NAI|) 

n©t  .  t  i-i«t  1 1 

10  |F(NA| I )  LE  01  GO  TO  .5 

IF  ( 1 1 1  .NAI I ) )  OC  TML  -T  >  DT IHAM  TM. -TT  1 1  .NAI  I  1  n  ♦ATXFRl  I  I 

16  IF  i  Tl  I  .NBI I » I  OC  TH)-ti  OT|NB«lTNU-TTi )  ,)OI  |)  1  ferXCRl  I  I 

?o  ao  ifiotina  or  0TIM6  on  otina  gt  dy inr i  so  to  zs 

COTINEI l)*0T!NA 
NAT  I  I >«NA< ! !-l 

:r'!  C,'  NIOINT.OR  NAI I  I  )£  NINTI  00  )0  05 
TlSTlf  '•TINIDINT.li 

?5  ST)Ctr>rifC  •otina 

?5  !C0!'H6*  GT  OTIHA  op  OTil«  OT  OTIC  GO  '0  50 

OOTDtf  !  )  -OT  if® 

IF  I  I  £0  MIOINT)  TftXIAtll  I  ■>•?!*:  •OOTlW  ■  > 

IF  I  I  EC  NIOINTl  T  IXITf)l  1  )  I  «T  I  I  ll 
30  N0f  i  I-N53' i  > -I 

50  CONTINUE 

OOOT  DC -OOT I  ME  I  I  > 

OO  60  i -a. HINT 

60  C00!  I  HE  *  AN  i  S  i  .TCEI!  1C  .00'  !HF.  ■  I  1 1 

*•  .i  iooo’ic  cc  otic  i  fo  io  gy 

onc«ooci  ic 

CAU  'iCHMIUf 
OO  Pt  r  mm 
,  Mfl 
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>C!HI  Sii3i3WI  1H30  0000000*' 
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UITC  VECTOR  IS  TO  TIC  LEFT  (fOR  TEMPERATURE  ARRAY  ABOVE  )  AX)  IS  UPWARD  (FC«  PRINTOUT  WHICH  FOLLOWS) 


ELEfCNT 

0  (IN) 

Q  (CNVC) 

Q  (RAD) 

Q  IVPR)  CMS  WRZO 

CHS  LEFT 

1 

.5503922*00 

-  .26!7T'2E-01 

. 1941892-03 

0. 

0. 

.1000002*00 

2 

. 1031692*0! 

-.435U9';E-0I 

.2657992-03 

0. 

0, 

.1000002*00 

3 

!  339052*0! 

-. 73346®. -01 

. 455447E-03 

0. 

0. 

.  1000002*00 

4 

.31 17512*0! 

-. 11641 42*00 

.9719932-03 

0. 

0. 

. 1000002*00 

5 

.5025662*01 

- . 1 767732*00 

3471582-02 

0. 

0. 

.  1 1 10962*00 

6 

.7704562*01 

- .2599462*00 

■266716E-0I 

0. 

0. 

•  16668m’  *00 

7 

. !  123242*02 

- . 3420512*00 

5960202-01 

0. 

0. 

.  I7543S  00 

8 

. ! 5S729E  *02 

- .4157992*00 

.91759762-01 

0 

0. 

.1687082*00 

9 

.20532 IE *02 

-.4806252*00 

.  14 -.5692*00 

0. 

0 

.1552432  *00 

10 

.2574352*02 

-.5319872*00 

.  I92i:v-E*D0 

0. 

0. 

.  I28986£*r  J 

1! 

.3069542*02 

-.5704792*00 

.2387660*00 

0. 

0. 

.9829112-01 

12 

. 3480552  *02 

-. 59586 3E *00 

272630E*70 

c. 

0 

.6775572-01 

13 

.3753132*02 

-.6098592*00 

.2920272*00 

0. 

0. 

4202222-0! 

IT 

.3842612*0? 

-.6128762*00 

.2916172*00 

0 

0. 

.3039762-01 

IS 

.  TT..2 :  3E  *02 

-.6061552*00 

273929C  *00 

0. 

0. 

.3(71202-01 

16 

.3480552*02 

-.5888242*00 

.2396842*00 

L 

0. 

4739122-01 

17 

.3069542*0? 

-.5601022*00 

.  193I06E*00 

0 

0 

.7627142-01 

18 

.2574352*02 

-.516850E*00 

.1343712*00 

0 

0 

.9995872-01 

19 

.20532 IE *02 

-.4580I6E*00 

.7557902-0! 

0 

0. 

.1000002*00 

20 

.1557292*02 

- . 3891302*00 

. 3679892-01 

0 

0. 

.1000002*00 

21 

. 1 123242*0? 

-.3149392*00 

.  162948C  01 

0. 

0. 

I00000E*00 

22 

.7704562*01 

-.24O03OE*QO 

6380172-02 

0. 

0. 

. 1000002*00 

23 

.5025662*01 

-  I72978E  *00 

.2488762-02 

0. 

0. 

.1000002*00 

24 

.31 17512*01 

- . 1 16333E  *00 

. 97 1 44  IE -03 

0 

u 

.1000002*00 

25 

1839052*01 

-.7334422-01 

. 455447E  *03 

0. 

o 

1000002*00 

26 

.1031692*01 

-.4354732-01 

265800E-03 

0 

u 

■100000E.00 

27 

550392E  >00 

-.2617582-01 

. I94I89C-03 

0 

0 

1000002*00 

OA 
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TOTAL  INCIDENT  ABSORBED  ENEROV*  429; 1900637067  CALORIES 
TOTAL  CALORIES  LOST  ac  TO  CONVECTION*  -8.9640259259  CALORIES 
TOTAL .EfCROY-EHITTEO-OOE  TO  RADIATION*  2 1602791 5860  CALORIES 
CALORIES'EXMNDED  IN  VAPORIZATION  (TGTAL>«  0.0000000000 


TIM  'id  MELT  THROUGH* 

DISTANCE  TROHC/L  (CMS) 
-2.60000 

-2.40000 

•2.20000 
-s. 00000 
-1.80000 
-I .60000 
-t.HOOOO 
-1.20000 
-1.00000 
-.80000 
-.60000 
-.40000 
-.20000 
0.00000 
.20000 
.40000 
.60000 
.80000 
1 .00000 
1.20000 

1.40000 
1 .60000 
1 .80000 
2.00000 
2.20000 

2.40000 
2.60000 


.4027470967  SECONDS 

TOTAL  THICKNESS  (CMS) 
. 100000000000 

.100000000000 

.100000000000 

. lOOOOOOOOOOO 

.111096293699 

. 166682630527 

. 175485097892 

. 168707588854 

. 155242690182 

.128984950525 

.096291138188 

.067755664368 

.042022169830 

.030397570922 

.031712024201 

.047391247646 

.076271423712 

.059958709493 

. 100000000000 

.100000000000 

. 1 00000000000 

. 100000000000 

. 100000000000 

. 100000000000 

. 100000000000 

.100000000000 

. 100000000000 


NON-MOLTEN  THICKNESS 
.  1000.00000000 

.100000000000 

.tOQOOOQOOOGO 

.  1000.70000000 

.lit  096293659 

. 140000000000 

. 140000000000 

.133333333333 

. 120000000000 

.093333333333 

.066666666667 

.033333333333 

.0)3333333333 

O.O.'.OOOOOOOOOO 

. GO s66CG66667 

.02666666666? 

.060000000000 

. 093333333333 

. 100000000000 

. lOOOCOOOOOOQ 

. 100000000000 

. 100000000000 

. 1000000000C0 

. 100000000000 

. 100000000000 

.100000000000 

. 100000000000 


6AW/HC'/72r17 


TOTOATUBE  or  BACK  SUTACE  10EG.K)  VS  II*  (SECS) 

TIME  (SECOOS)  TDTERATuRE  IDE  OWES  KELVIN) 


1. 

.30000000006*03 

.  £3£6469 1 68E  •  00 

.16400705516*04 

.£5413405006*00 

.17468738516*04 

.£7090670496*00 

.18097586776*04 

■£940S46£94C*00 

.  I956£4575IE*04 

.£974031  £006 *00 

.  IB9O7£6403E*O4 

■  300n»7060E*p0 

.  19)70766896 *04 

.3£0l£3STSE«uC 

.19396166636*04 

.3304014*1066*00 

.  19545184476*04 

-3403305£I7E*00 

. I 9664788976*04 

•34974|9££0C*00 

, 197800531 IE*94 

. 36690038686  *00 

.19747010096*04 

■  367959004£E*00 

.  1 973008 184E *04 

•  3769947£8lE*00 

.  I97S0564?9E*04 

.39605465736  *00 

.  19754486966*04 

■  395£5I0094E*00 

.19779816516*04 

. 408747Q9676  *  00 

. 19868198896*04 

SOL  ID-U  QUID  INTERFACE  POSITION  I  CMS  I  VS  lift  (SECS) 
CrtS  FROM  ORIGINAL  SURFACE 


II*  (SECONDS) 
.83864091666*00 

.£5413405006*00 

•  £705067049E*00 

.£040545894E*66 


0. 

•  6666666667E-08 
.  13333333336-0! 

.£ooo6oooooe-oi 


.£97403i£O86*OO 
.30003970606*00 
.  3£0 1  £357856  *00 
.  330461 47G€£*00 
•  3403305£I7C*00 
.34974196806*00 
.3509oo3a6e:*oo 
. 36795900466  *00 
.  3769947891  E*00 
.  386054657 7t  *00 
. 39565 1 0094C  *00 
.40£7470967E*00 


.£6666666676-01 
. 33333333336-Oi 
.400000G300E-01 

■  4666666667E-0! 
. 533333333 3E-0 I 
•6000000000E-C1 

•  66666G6667E-0I 
. 7333333333E-0 I 
•0OOOOOOOOOE-OI 

•  63S6666S67E-0I 

■  9333333333E-0I 
. 10000000006*00 


SURFACE  RECESSION  RATE 

I  CMS/SEC  1  VS  II* 

TIHISECOWSI 

VELOCITV(CMSTSEC) 

.843391907 4C*00 

.31089041806*00 

. £6838070 I4C *00 

. 407£0309066*00 

. 87760061 78E*00 

.  46464651 49E*  00 

.£91 16068516*00 

.58790863766*00 

•303l6l4|34E*00 

.50703091976*00 

•3I440I64£I£*CO 

.5906I37373C*00 

. i£53O£5844€*00 

.64363157806*00 

.  33540599686  <00 

.676EB4 15986 <00 

. 3430363£I9E*00 

•  7CB3607605C*00 

•  3543£l  I54'*E«00 

■7£79£404l8£*00 

■  3634£S69MC*00 

.73594747676*00 

■  37£4766WcE*00 

.73781813006*00 

.  36158469876*00 

73564091976*00 

.39053883336*00 

.  7249851 8696*00 

.  39699905306  *  00 

.08335835766*00 

I  SECS*  AM)  POSITION  OF  SaiO-LIOUIO  INTERFACE 
CENTI*TERS 
.  3;'T)333333E-0£ 

. lOOOQOOOOOE-OI 
.<6666666676-01 
•8333333333C-01 
.30000000000-01 
■  3666666667E  -Cl 
43333333336-0! 

.50000000006-01 

.56666665676-0! 

6333333333E-0I 

7000000000C-0! 

.  76566666676-01 
8333333333C -01 
SOOOOOOOOOf-OI 
96666666676.  01 
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THICKNESS  (CMS) 
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Fig.  A-l.  A  Cross-Sectional  View  at  Time  of  Melt-Through  of  a  0.01 
cm  Thick  Titanium  Sheet  Subjected  to  a  Mach  0.5  Sea-Level 
Wind  and  an  Absorbed  Gauss ian-Distributed  Flux  of  10  kw 
Over  a  10  cm2  Spot. 
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TEMPERATURE  OF  BACK  SURFACE  ( DEG  K)  *tOl 
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Fig.  A-2.  The  Temperature  History  at  the  Back  Surface  (Opposite  the 
Laser  Beam  Axis)  of  a  0.1  cm  Thick  Titanium  Sheet  in  a 
Mach  0.5  Sea-Level  Airstream  Absorbing  10  kw  of  Laser 
Power  Over  a  10  cm2  Spot 
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TIME  TO  MELT  (SECS) 
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Fig.  A-3.  The  Solid  Surface  Recession  History  at  the  Center  of  the 
Laser  Beam  Spot  of  a  0.10  cm  Thick  Titanium  Sheet  Absorb¬ 
ing  10  kw  of  Laser  Power,  Gaussian-Distributed,  Over  10 
cm*  (Time  vs.  Axial  Distance) 
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Fig.  A-4.  The  Solid  Surface  Recession  History  at  the  Center  of  the 
Laser  Beam  Spot  of  a  0.10  cm  Thick  Titanium  Sheet  Absorb 
ing  10  kw  of  Laser  Power,  Gausslan-Distributed,  Over  10 
cm2  (Axial  Distance  vs.  Time) 


The  Velocity  History  of  the  Nadir  of  the  Liquid-Solid 
Interface  in  a  Laser-Formed  Pit  in  a  0.1  cm  Thick 
Titanium  Sheet  Absorbing  10  kw  of  Gaussian-Distributed 
Power  Over  a  10  cm2  Spot 
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of  the  Liquid-Solid  Interface  In  a  L»:ar-Formed  Pit  in  a 
0.1  cm  Thick  Titanium  Sheet  Absorbing  1C  kw  of  Gaussian- 
Distributed  Power  Over  10  cm2 
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APPENDIX  B 
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SUBROUTINE  ATWOSir.TM.SlCHA.FNO, THETA.OtXTA.CA.ANU.KI 


PACE  l 


..□El 


s  i  qua*  s  i  GtiAS  m*.  *exp(  -<  o*  i  K-Hpn  i  mb  <  mi  i 


ShO»«HOZ*SIUNA 

'CTA-TM/TM? 

TCI  ?  "SIGMA*  WT A 

CA-^9  *■  >  77«$ORTlTM) 

_ a , 

<  TVCTA*  *3M  I  (  TM7*SI  / «  TH»Si  > 
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b^0UTi»c  ruuxx  <rLux.siCMA?» 
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SJMUTIfC  nVTItti 


'MX  ) 
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SUBROUTINE  SCHMOT 


HO 


HNET  ( I  .  Jl  *HXFRL  <HXFRR<HXFRB<HXFRT 
TC  |  , JI*TTI  I .  Jl  <OT|ME/RHO/CPIJ  OYI  I  .J1/DCLX‘HNET!  I  .  Jl 


1  , 

J.EQ.NAI 1 1 

Li 

NO 


J.EQ.N61 1 1  [.IESljI  0rxFR(  1  l*RHO*CPIJ,OY! !  .  Jl  ’DELX/HNET 1 1  .Jl 


IOO 


CONTINUE 
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SUBROUTINE  TPNSfRx 


PA.GC  1 


SU6R0UHIC  IRNSfPX 


PACE  2 
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subroutine  xactime 
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APPENDIX  C 

Thermophysical  and  Thermodynamic  Properties  of  Aluminum,  Magiiesium, 
Titanium,  Stainless  Steel  (AISI  304),  and  Air  as  Used  in  QUEST 

Two  sources  were  primarily  used  in  obtaining  the  thermodynamic 
and  thermophysical  properties  of  the  four  metals  studied.  The  first 
(Ref  17),  from  which  specific  heats,  heats  of  fusion,  and  heats  of 
vaporization  were  taken,  is  Thermodynamic  Properties  of  65  Elements— 
Their  Oxides,  Halides,  Carbides,  and  Nitrides  compiled  by  Charles  E. 

Wicks  and  Frank  E.  Slock  and  published  in  1963  as  Bureau  of  Mines 
Bulletin  605.  This  was  considered  as  the  most  complete  and  authori¬ 
tative  source  available  for  the  calculation  of  heat  content  as  a  function 
of  temperature  from  room  temperature  (30Q°K)  through  the  boiling  point. 
The  second  source  (Ref  18),  used  to  provide  conductivities, .^densities, 
viscosities,  and  emissivities,  is  Recommended  Values  of  the  Thermo¬ 
physical.  Properties -of  Eight  Alloys,  Major  Constituent*--  ^nd,  Their 
Oxides  assembled  under  the  direction  of  Y.  S,  Touloukian  of  Purdue 
University.  This  was  the  most  recent  (1966),  the  most  thorough,  and 
in  most  cases  the  only,  source  that  could  be  found  which  presented  the 
listed  properties  as  a  function  of  temperature  over  the  critical  ranges 
of  interest.  However,  extrapolation  of  this  data  from  within  the  liquid 
state  to  the  boiling  temperature  was  neced-ary  to  cover  any  temperature, 
perhaps  unlikely  but  theoretically  possible,  that  might  be  encountered 
by  QUEST. 

Standard  atmospheric  properties  were  used  for  free  stream  air 
conditions.  Although  only  sea  level  conditions  were  used  in  the  cases 
tested  thus  far,  the  subroutine  ATMOS  (a  library  subprogram  obtained 
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i  '  ' 

i 

from  the  Computer  Center,  Aeronautical  Systems  Division,  Wright-Patterson 
AFB,  Ohio)  possesses  the  capability  to  compute  applicable  properties  of 
the  atmosphere  according  to  the  1959  ARDC  Atmospheric  Tables  to  altitudes 
of  300,000  feet.  Property  values  of  the  air  for  boundary  layer  calcula¬ 
tions,  including  conductivity,  specific  heat,  enthalpy,  and  viscosity, 
v/ere  obtained  from  Refs  19.,  20,  and  21. 

The  property  value  of  the  metals  as  used  within'' QUEST  will 
warrant  review  as  any  future  investigations  and  experiments  which  pro¬ 
vide  a  better  basis  of  estimation  . become  available.  The  present 
estimates,  however,  obtained  from  the  limited  sources  cited  and  while 
certainly  in  error  to  some  extent,  are  deemed  sufficient  to  provide 
meaningful  trends  and  comparisons  between  the  modeled  metals. 

The  property  values  as  used  in  QUEST  for  aluminum,  magnesium, 
titanium,  stainless  steel,  and  air  follow. 


4 


125 


GAW/MC/72-17 


ALUMINUM 

Density:  2.698  gm/cm3 

Melting  Point  (M.P.):  933°K 

Boiling  Point  (B.P. ):  2740°K 

Hteat  of  Fusion:  95.185  cal/gm 

Heat  of  Vaporization:  2516.667  cal/gm 

Conductivity 

T  < M. P.  k  =  [0.5325466  +  2*(-l. 839836x10““ 

+  2’(-2.360679x10"7).)3  cal/(sec-cm*°K) 

A/,?.  <T  < B.P.  k  -  [0.21749  +  Q.;P00717(2’ -  1000°K)]  cal/(sec*cm* 

Specific  Heat 

cp  =  [0.1831  +  2’(1 . 097048x10  "’l ]  cal/(gm°K) 

Cv  =  0.259  cal/(gm°K) 

Viscosity 

y  =  0.07612433  +  r[-l .460141X10"“  v  2?(l.24557><10"7  +  ff{-5. 404153X10 
+  21[1.181837^10~14  +  Tf-l .035147*10" 1 °)]} )]  poise 

Hemispherical  Total  Emissivity 

C  =  0.022  +  o.oooo6(:r) 

NOTE:  T  is  temperature  expressed  in  °K. 


T  <  M.  P. 
T  >  M.  P. 
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MAGNESIUM 


Density:  1.738  gm/cm3 

Melting.  Point  (Af.P. ):  923°K 

Boiling  Point  (S.P. ):  1393°K 


Heat  of  Fusion:  88.845  cal/gm 

Heat  of  Vaporization:  1295.656  cal/gm 

Conductivity 

T  <M.P. 

k  =  0.37762  +  ^C-O. 0000478)  cal/ (sec- cm- °K) 

M.P .  <?<B.P. 

k  =  0.1673  cal/(sec-cm-°K) 

Specific  Heat 

T  < M. P.  C  =  0.2S25S  +  !T(6.1697927X10'5)+  S’"2 (3208. 2922)  cal/ (gm»°K) 
P 


T  > M. P.  Cp  =  0.32  cal/ (gin- °K) 

Viscosity 

U  =  1 .460654^10" 1  +  T{-2 .911014x10"“  +  2’[2.0358xi0 
+  n-4. 880739x10" U)]}  poise 

Hemispherical  Total  Emissivi.ty 

e  =  0.08  +  0.000178  < ’T) 


NOTE:  T  is  temperature  expressed  in  °K. 
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TITANIUM 


Density: 

Melting  Point  (M.P.): 
Boiling  Point  (s.p. ): 
Heat  of  Fusion: 

Heat  of  Vaporization: 


4.50  gm/cm3 
1988.0°K 
3550. 0°K 
93.946  cal/gm 
2108.56  cal/gm 


Conductivity 


T  <  M.  P.  k  0i075251  +  T{-0. 00010874  +  i{l. 353216  10-7 

+  T.f-0. 294976  xlO-1 1  +  2’(1.0ll3p52xl0~1 ")  ]}) 
♦  cal/ (sec* cm* °K) 


M.P,  <T  <B.P.  k  =  0.035 


Specific  Heat 

T  < M.P.  Cp  =  0.071942  +  r[0. 0002821  +  2’(-5. 009421x10 
+  T{5.0S938S*10~10  +  r[-2.400367xi0"n 
+  y(4. 327858*10" 1 '}]})]  cal/(gm*°K) 

T>M.P.  Cp  =  0.167  cal/(gm*°K) 

Viscosity 

*  nnc  ^  -  Z*)/125] 

p  =  0.005  +  0.06(3^  1  poise 

Hemispherical  Total  Emissivity 

:T<950oK  C  =  0.08  +  1 .8181822’ 

T  >  950°K  c  =  0.20  +  1.1764  721 


NOTE:  T  is  temperature  expressed  in  °K 
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StAINLESS  STEEL  (AISI  304} 


density: 

Melting  Point  [M.P. ): 
Boiling  Point  [B.P. ): 
Heat  of  Fusion: 

Heat  of  Vaporization: 


7.9  gm/cm3 
1730.0°K 
3000. 0°K 
$6,252  cal/gm 
1515.211  cal/gm 


Conductivity 

T  < M. P.  k  =  .0.0155816  *  2'[l.  039741x10'“  +  f(-l .579824x10 

+  H1.653638X10-10  +  ?[-8.215264xl0"1“ 

+  r(i,. 565207X1 0"17)].})]  cal/(sec*cm*°K) 

M.P.  <T  < 3.P.  k  =  O.OS258  cal/ (sec  ♦  cm  •'’K) 


Specific  Heat 

T  < M. P.  Cp  =  0.06496569  +  Hi. 696444X10"“  +  ^[-1 . 119582xi0'7 
+  T[2 .,769178X10" ")]}  cal/(gm*°K) 

T  > M. P.  Cp  =  0.1791  cal/(gm-°K) 

Viscosity 

M  =  2.448646  +  2’[-4.22152xl0"3  +  2’(2.?87352xi0-6  <-  T{-1 .067079xl0"3 
+  3,[1.911987xi0-13  +  jT(-1 . 369543X10"1 7)  ] })]  poise 

Hemispherical  Total  Emissivity 

200°K  <T  <  700°K  £  =  0.095  +  0.0002(2’) 

2’>700°K  e  =  0,165  +  0.0001C21) 


NOTE:  T  is  temperature  expressed  in  °K. 
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AIRo 

t 

Standard  atmospheric  properties  at  sea  level: 
f  =  2S8.15°K 

Speed  of  sound  =  340.29205  m/sec  =  1,116.89  ft/sec 
p  -  10.1325  nt/cm2  =  14.696  Ib'/in2 
p  =  0.001225  gm/cm3  =  3.7457  slug/ft3 

Conductivity 

ko  =  5.88218832  cal/(sec‘cm,sK) 

To  =  500°R 

k  =  ko  [0.0506331  +  1. 07642 (T/T0)  -  0. 113939 (T/2V) 2 
+  0.00679309(2’/2’o)3]  cal/ (sec ‘.cm- °K) 

Specific  Heat 

Cpo  =  0.2395  cal/(gm*°K), 

To,  =  500°K 

Cp  =  Cpo[0.979  +  0. 00651928  (T/To)  +  0.0158S36(2’/2,0)2 
-  0. 00155636 (SyS’o)'3]  cal/(gra’°K) 

Enthalpy 

h0  -  66.377778  cal/gm 
To  =  500°R 

h  *  66 . 37778 [0 .01836064  +  0.9395396  (S’/?,,) 

+  0. 03243428 (f/Jo) 2  -  9.91742X10”4 (SP/fo)3]  cal/gm 

Viscosity 

Po  =  2.55936xl0‘"4  gm/(sec*cm) 

To  -  9Q0°R 

V  -  Uo  V(T/To)  gm/ (sec* cm) 


mwm-w 
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APPENDIX  D 


Numerical  Results  from  QUEST  for  0.05,  0.10,  and  0.20  cm  Thick  Aluminum 
Magnesium,  Titanium,  and  Stainless  Steel  Irradiated,; by  1.0  and  10.0  fa 
Power  Over  Spots  of  1.0  and  10.0  cm2  at  Selected 
Mach  Numbers  Between  0.0  and  4.0 
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TABLE  D-I 

Melt-Through  limes  (in  Seconds)  as  Predicted  by  QUEST  for  Aluminum 


Spot  Size  (cm2) 

r  " 

Mach  Number 

(at  sea  level) 

1 _ : _ 

1.0 

m’. 

10.0 

Plate  Thickness  (cm) 

Plate 

Thickness  (cm) 

B9 

0.10 

mm 

0.05 

0.10 

0.20 

" 

1.0  kw:  Absorbed  Power  - 

0.0 

0.03535 

0.20194 

0.52967 

0.82320 

1 .86313s 

4.41704 

0.1, 

0.08540 

0.19978 

- 

- 

- 

- 

0.5 

0.08477 

0.19316 

0.47043  , 

- 

- 

4.92163 

1.0 

0.08444 

0.19148 

0.46699 

- 

2.19655 

5.32153 

1.5 

- 

- 

- 

2.24050 

- 

2.0 

0.08377 

0.18977 

0.46090 

- 

2.17886 

5.24222 

3.0 

i 

- 

1.77671 

- 

4.0 

0.07939 

0.17682 

0.41828 

- 

1.31212 

2.83614 

10.0  kw  Absorbed  Power 


0.0 

0.00934 

0.02190 

0.05554 

0.07519 

0,15718 

0.34103 

0.1 

- 

- 

- 

0.07533 

0.15717 

- 

0.25 

- 

- 

- 

0.07540 

- 

- 

0.5 

0.00917 

0.01997 

0.04173 

0.07531 

0.153111 

0.31521 

2.0 

0.00846 

0.01678 

0.03313 

0.07448 

0.15085 

0.30828 

4.0 

0.00791 

0.01582 

0,03112 

0.07121 

0.14386 

0.29477 

¥ 


* 


* 


t  '' 


) 
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TABLE  D-II 


Melt-Through  Times  (in  Seconds)  as  Predicted  by  QUEST  for  Magnesium 


' 

_ 

Spot  Size  (cm2) 

Mach  dumber 
(at  sea  level) 

1.0 

10.0 

Plate  Thickness  (cm) 

Plate  Thickness  (cm) 

0.05 

0.10 

■B 

0.05 

0.10 

0.20 

1 .0  kw  Absorbed  Power 


0.0 

0.05846 

0.13485 

0.34751 

0.55793 

1 .22602 

2.89378 

0.1 

0.05850 

0.13349 

0.32574 

- 

- 

- 

0.5 

0.05794 

'0.12797 

- 

0.60350 

1.33226 

3.16028 

1.0 

0.05750 

0.12579 

0.29493 

- 

1.40712 

3.37728 

1.5- 

- 

- 

- 

0.64631  f 

1.43397 

- 

2.0 

0.05682 

0.12434 

0.28970 

- 

1.39618 

3.30165 

3.0 

- 

- 

- 

- 

1.16482 

2.60833 

4.0 

0.05400 

0.11659 

- 

- 

0.87862 

1.87628 

10.0  kw  Absorbed  Power 


0.0 

0.00727 

0.02108 

0.06527 

0.05332 

0.11212 

0.24451 

0.1 

- 

- 

- 

0.05341 

0.11210 

- 

0.5 

0.00694 

0.01540 

0.03087 

0.05332 

0.10934 

2.0 

0.00626 

0.01232 

0.02439 

0.05238 

0.10516 

4.0 

0.00586 

0.01156 

0.02203 

0.05005 

0.10071 

0.20484 
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TABLE  D-III 


Melt-Through  Times  (in  Seconds)  as  Predicted  by  QUEST  for  Titanium 


* 

Spot  Size  (cm2) 

Mach  Number 
(at  sea  level) 

1.0 

I- 

10.0 

Plate  Thickness  (cm) 

Plate  Thickness  (cm) 

- 

■MB 

■  mm 

0.10 

0.20 

0.05 

0.10 

0.20 

1 .0  kw  Absorbed  Power 

0.0 

0.20763 

0.46671 

1.15179 

1;  99589 

4.22383 

’  9.4575 

0.1- 

0.20713 

0.44930 

0.99254 

- 

4.40521 

9.7064 

0,5 

0.20257 

0.42399 

0.91670 

2.36293 

5.01379 

11.2717 

1 .00 

0.201312 

0.421019 

0.903722 

2.73299 

5.88005 

13.7029 

2.00 

0.20065 

0.41977 

0.91325 

3.51087 

7.81713 

20.3557 

4.00 

0.19664 

0.41074 

0.89954 

3.23944 

6.98538 

- 

6.00 

- 

( _ _ 

- 

- 

3.16585 

- 

10.0  kw  Absorbed  Power 


0.0 

0.03201 

0.10184 

0.30469 

0.19696 

0.41948 

0.92913 

0.1 

- 

- 

- 

0.19754 

0.41346 

- 

0.25 

- 

T 

- 

0.19677 

- 

- 

0.5 

0.02577 

0.05213 

0.10236 

0.19542 

0.39231 

0.77309 

2.0 

0.02141 

0.04211 

0.08358 

0.19304 

0.38146 

0.76627 

4.0 

0.02075 

;  0.04023 

- 

0.18873 

0,37820 

0.76285 
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TABLE  D-IV 

Melt-Through  Times  (in  Seconds)  as  Predicted  by  QUEST  for  Stainless  Steel 


"  ’  '  ' 

-  . 

1 

Spot  Size  (cm2) 

V 

Mach  Number 
(at  sea  level) 

1.0 

10.0 

Plate,  thickness  (cm) 

Plate  Thickness  (cm) 

■ 

0.10 

0;20 

0.05 

BBS 

mm 

’  1 

1.6  kw  Absorbed  Power 

0.0 

0.2601.8 

0.58389 

' 

1.43667 

2.50173 

5.28696 

11.7950 

0.1 

0.26056 

.  0-57556 

1.30990 

- 

- 

- 

0.5 

0.257121 

0,54425 

■‘-18365 

2.90293 

6.16154 

13.7963 

1.0 

0.25493 

0.53390 

1,14413 

3.27428 

7.02708 

16.2477 

2.0 

0.25276 

0.53106 

.  1.14560 

3.91964 

8.60810 

21,2044 

3.0 

- 

- 

- 

- 

8.64554. 

21.2096 

4.0 

i 

0.24664  , 

0.51299 

:  1.1120 

3.21920 

6; 81 835 

.15.3470 

10.0  kw  Absorbed  Power 

■ 

0.0 

0.03760 

0-11804 

0.36298 

0.24689 

0.52481 

1,15785 

0.1 

- 

- 

- 

0.24769 

0.52489 

- 

0.5 

0.03494 

0.07705 

0.14304 

j 

0.24729 

0.50430  | 

0.99877 

1.0 

- 

- 

- 

- 

- 

0.96467 

2.0 

0.03023 

0.06090 

0.10138 

| 

0.24407 

0.47639 

0.95862 

3.0 

- 

- 

- 

- 

0.95426 

4.0 

0.02709 

0.05197 

0.09482 

0.23533 

0.467544 

0.94398 
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